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Chapter  1 


INTRODUCTION 

An  important  part  of  the  analysis  of  signals  is  the  determination 
of  energy  as  a  function  of  frequency.  The  usual  approach  to  this  task 
is  to  calculate  the  classical  Fourier  Transform  of  the  signal  and  to 
investigate  the  square  of  its  magnitude.  This  approach  has  been  used 
for  years  to  analyze  all  kinds  of  signals.  Yet,  for  most  of  those  who 
do  energy  spectral  analysis,  something  is  lacking.  It  is  natural  to 
think  of  the  frequency  distribution  of  energy  as  a  changing  function  of 
time.  Our  perception  of  speech,  communications  signals,  physiological 
signals,  and,  in  fact,  all  continuous  information-bearing  signals  prompts 
this  interpretation.  There  is  a  process  taking  place  within  the  signal; 
energy  within  the  signal  is  sensed  to  be  distributing  itself  over  dif¬ 
ferent  frequencies  as  time  goes  on.  The  Fourier  Transform,  however, 
cannot  provide  the  information  for  such  a  time-frequency  energy  distri¬ 
bution.  By  its  very  definition,  the  Fourier  Transformation  uses  the 
entire  signal  and  permits  analysis  of  only  the  frequency  distribution  of 
energy  of  the  signal  as  a  whole.  To  solve  this  problem,  many  who  do 
spectral  analysis  have  taken  a  "piecewise"  approach;  the  signal  is 
broken  up  into  contiguous  pieces  and  each  piece  is  separately  Fourier 
Transformed.  The  resulting  family  of  Fourier  Transformations  is  then 
treated  as  if  it  were  the  basis  for  a  joint  time-frequency  energy 
distribution.  Despite  this  practice,  a  theoretical  justification  for 


this  approach  cannot  be  found  in  the  literature.  On  the  other  hand, 
there  are  some  interesting  approaches  which  do  have  some  theoretical 
justification  but  they  suffer  from  serious  difficulties  of  interpre¬ 
tation  and  implementation. 

The  main  contribution  of  this  work  consists  of  showing  that 

(1)  There  ^s  a  theoretical  justification  for  treating  the 
piecewise  Fourier  Transform  approach  as  a  valid  basis  for 
joint  time-frequency  energy  analysis,  and 

(2)  there  is  a  general  family  of  "windowed"  Fourier  Series 
expansions  which  provides  a  useful  framework  for  applying 
"weighting"  windows  to  time-frequency  energy  spectral 
analysis. 

In  addition,  this  work  provides  a  sound  justification  for  using  windows 
with  the  Finite  Fourier  Transformation  along  with  some  of  their  limi¬ 
tations.  A  new  perspective  on  signals  is  also  given  through  the  concept 
of  a  signal  vector  moving  through  a  vector  space  whose  basis  vectors 
provide  measurements  of  signal  energy  in  time  and  frequency. 

In  order  to  standardize  notation  and  clarify  assumptions,  the 
following  explanations  are  given.  All  signals  discussed  are  assumed  to 
be  complex  and  square  lntegrable.  That  is,  the  signal  x(t)  is  assumed 
to  have  finite  energy  E  given  by 


E 


dt  <  •  . 


(1-1) 


Observations  are  assumed  to  be  made  of  the  signal  beginning  at  a  finite 


start  time,  t.,  and  ending  at  a  finite  stop  time,  t_ .  Outside  of  the 
A  B 

time  Interval  ,  the  signal  Is  assumed  unknown.  Furthermore, 

x(t)  Is  assumed  to  be  "piecewise  smooth,"  meaning  the  signal  x(t)  and 
its  first  derivative  dx(t)/dt  are  both  piecewise  continuous.  This  Is  not 
an  overly  restrictive  requirement;  most  real-world,  physical  signals  are 
piecewise  smooth.  The  Fourier  Transform  of  the  signal  x(t)  is  defined 


*  |  x(t)  e-^*1*  dt 


(1-2) 


where  the  radian  frequency,  id,  is  related  to  the  cyclical  frequency, 
f,  by  the  expression 


ai  -  2irf  . 


Alternatively  the  Fourier  Transform  is  frequently  written  as 


(1-3) 


X(f)  -  |  x(t)  e~J2irft  dt. 


(1-4) 


Since  the  right  sides  of  (1-2)  and  (1-4)  are  identical,  X(f)  and  X(u>) 
represent  Identical  quantities.  Furthermore,  the  transform  in  (1-4) 
should  be  denoted  by  X(2irf)  in  order  to  be  consistent  with  X(u>) 
defined  in  (1-2).  However  following  common  practice,  the  notation  X(u>) 
and  X(f)  will  be  used  for  convenience.  Strictly  speaking,  X(u>)  and 
X(f)  are  two  different  functions,  the  difference  being  a  scale  factor 


to.. 


in  the  Independent  variable.  The  term  "energy  spectrum"  will  denote 
the  squared  magnitude  of  either  X(f)  or  X(w) ,  depending  upon  the 
context.  Although  use  of  the  Inverse  Fourier  Transform  will  not  be 
made  In  the  discussion,  It. Is  defined  for  a  signal  x(t)  as 


-  |  X(f)  eJ2lrft  df 


(1-5) 


or  by 


X(o>)  eJuJt  dw  . 


(1-6) 


An  Important  result  Is  Parseval's  Formula  given  by 


|x(t)|2  dt 


•srf  lw«>l: 


(1-7) 


or  by 


|x(t)|2  dt 


|x(f) 


(1-8) 


This  relationship  Is  the  basis  for  the  Interpretation  of  |x(u>)|2  or 
I  X(f)|2  as  the  frequency  distribution  of  energy  or  as  the  energy 
"density"  In  frequency  of  the  signal  x(t).  Since  the  Dirac  Delta 
function  Is  often  used  In  this  context.  It  Is  defined  here  by  Its 
familiar  Integral  property.  Specifically,  the  Dirac  Delta  function  of 
time,  6(t),  is  defined  by 


x(t)  6(t-t  )  dt  ■  x(t  )  . 


(1-9) 


Whan  Che  limits  of  Intagratlon  srs  finits,  thsrs  results 


r  x<to> 


<  t 

o 


< 


H 


|  x(t)5(t-to)  dt 

*A 


*  x(to+)  8  eo  "  eA 
51  x(to_)  5  fco  "  *B 


I  0  :  t  <  t.  or  t  >  t_, 

y.  *  o  A  o  B 


(1-10) 


where  It  Is  assumed  that  t.  <  t_  and  the  notations  t  +  and  t  denote 

A  B  o  o 

the  limit  at  tQ  taken  from  the  right  and  left,  respectively.  A  counter 
part  to  the  Dirac  Delta  function  used  with  discrete  functions  is  the 
Kronecker  Delta,  6  ,  defined  by 


(1-11) 


Both  of  these  functions  will  appear  in  the  discussions.  Finally,  it  is 
useful  to  define  the  rectangle  function,  Rect C^) »  as 


Rect(|) 


(1-12) 


The  Fourier  Transform  of  this  function  is  often  known  as  the  "sine" 
function  of  cyclical  frequency  f.  In  particular,  its  definition  is 


With  this  as  background,  a  brief  outline  of  the  dissertation  is 
given  next.  Chapter  2  discusses  the  approaches  that  have  been  proposed 
in  the  past  for  time-variant  energy  spectral  analysis.  The  original 
approach  was  proposed  by  Gabor  [1]  in  1946.  Although  many  concepts 
were  proposed  over  the  years,  it  seems  that  only  one  took  hold,  the 
"piecewise"  Fourier  Transform  approach.  Five  significant  questions  are 
posed  in  Chapter  2  which  reveal  the  inadequate  theoretical  development 
of  this  popular  approach.  Chapter  3  presents  the  proposed  vector  space 
interpretation  of  time-variant  energy  spectral  analysis.  A  clear, 
useful  concept  is  presented  which  nicely  lends  itself  to  a  theoretical 
justification  of  the  "piecewise"  Fourier  Transform  approach.  In 
addition,  a  family  of  windows  is  shown  to  be  useful  in  providing  "sharp' 
spectral  representlons  of  signals.  Finally  the  five  questions  from  the 
end  of  Chapter  2  are  answered.  Chapter  4  considers  problems  arising 
in  practical  applications  of  the  approach  proposed  in  Chapter  3.  A 
summary  and  recommendation  for  future  work  is  given  in  Chapter  5. 
Lastly,  the  appendix  lists  the  computer  programs  used  in  Chapters  3 


Chapter  2 


PREVIOUS  APPROACHES  TO  TIME-VARIANT  SPECTRAL  ANALYSIS 

In  this  chapter  the  previous  approaches  to  the  problem  of  time- 
variant  energy  spectral  analysis  are  described.  By  time-variant  energy 
spectral  analysis  is  meant  the  process  of  describing  the  energy  of  a 
given  time  function  x(t)  by  a  two  dimensional  distribution  over  time  and 
frequency.  From  this  point  of  view,  the  ordinary  concept  of  the  fre¬ 
quency  spectrum  of  a  signal,  expressed  in  terms  of  its  Fourier  Transform, 
is  not  adequate  to  describe  the  time-varying  behavior  of  signal  energy. 

An  example  is  the  acoustic  signal  generated  by  a  musician  on  his  instru¬ 
ment.  Here  the  signal  consists  of  a  time  sequence  of  musical  "notes," 
each  "note"  being  en  acoustical  vibration  of  prescribed  Intensity  and 
frequency  and  lasting  only  a  short  time.  An  observer  interested  in  what 
is  happening  within  this  signal  gathers  very  little  information  from  an 
ordinary  Fourier  frequency  spectrum  which  is  based  on  the  entire  history 
of  the  signal.  This  inadequacy  of  the  ordinary  Fourier  spectrum  motivated 
researchers  to  look  for  ways  to  construct  an  "energy  spectrum"  which 
changes  as  the  signal  changes. 

The  concepts  presented  in  this  chapter  will  be  shown  to  have  serious 
difficulties  of  interpretation.  The  main  emphasis  of  this  chapter  will 
be  to  show  these  difficulties  and  to  set  the  stage  for  the  proposed 
approach  which  is  discussed  in  Chapter  3.  The  work  of  Gabor  [1],  who 
developed  a  representation  of  a  signal  in  terms  of  s  fundamental  building 


block  in  the  tine-frequency  plane  called  the  logon.  Is  first 
discussed.  The  significant  further  work,  based  upon  Gabor,  done  by 
Bastiaans  [ 2  ]  and  Lerner  1 3 ]  is  also  pointed  out .  Fallacies  associated 
with  Gabor's  concepts  are  shown  and  the  difficulties  with  implementing 
his  representation  are  explained.  The  works  of  Page  [4],  Levin  [5],  and 
Rlhaczek  [6]  to  establish  a  two  dimensional  "instantaneous"  energy 
spectrum  are  then  discussed.  Problems  that  arise  with  the  interpretation 
and  Implementation  of  these  concepts  are  also  shown.  Finally,  the  pop¬ 
ular  Discrete  Fourier  Transform  approach  to  time-variant  spectral 
analysis  is  discussed  and  the  recently-developed  Wigner  Distribution  [7, 
8,  9,  10]  approach  is  presented. 


2.1  Gabor's  Logon  Concept 

In  1946  Gabor  published  a  landmark  paper,  "Theory  of  Coununication" 
[1]  in  which  he  sought  to  represent  time  signals  in  two  dimensions,  with 
time  and  frequency  being  coordinates.  Be  pointed  out  that  our  own 
physical  perception  of  sound  corresponds  to  a  two  dimensional  function 
of  time  and  frequency  even  though,  strictly  speaking,  the  acoustical 
vibration  detected  by  the  human  ear  is  a  function  of  time  alone.  Gabor 
defined  an  "elementary  signal"  as 


and  its  Fourier  Transform,  F  (fi  t  ,  f  ),  as 

6  O  © 


Vfi  V  fo>  ■  \  V':  V  V 


-J2irft 


dt 


£  .V/Af-f  >*  + 


e) 


(2-1-1) 


(2-1-2) 


(2-1-3) 


This  elementary  signal  Is  just  s  complex  dsoidal  function  of  frequency 
fQ  end  phase  0  which  Is  modulated  by  a  Gaussian  function  centered  at 


time  tQ.  (See  Figure  2-1- la  and  Figure  2-1-lb.) 


FIG 


S- la.  REAL  PART  OF  ^  / 

titfi  *.,£)]  tf  .  _2L_  ’ 

■f  VTir 


Fig.  r-i^  real  part  of  Ye(f) 


The  spread  of  each  function  Is  defined  by  Its  standard  deviation.  The 


property  of  "reciprocal  spreading"  is  noted,  specifically. 


1  1 


t  2ir  of  * 


(2-1-4) 


where  ot  and  are  the  standard  deviations  of  their  respective  Gaussian 
envelope  functions  in  (2-1-1)  and  (2-1-3) . 

Gabor  chose  this  "elementary  signal"  as  a  foundation  on  which 
to  construct  an  approximation  for  an  arbitrary  physical  signal.  His 
choice  was  evidently  Influenced  by  the  fact  that  this  function  has  a 
unique  property.  Of  all  well-behaved,  physically -meaningful  functions, 
this  is  the  only  function  for  which  the  product  of  "effective  duration" 
At,  and  "effective  bandwidth,"  Af,  which  are  defined  below,  takes  on 
its  minimum  value  of  %.  These  "effective"  terms  are  essentially 


2-4 


measures  of  the  mean  square  deviation  from  time  and  frequency  center- 
points,  respectively.  Specifically,  for  an  arbitrary  complex  signal 
x(t)  with  Fourier  Transform  X(f)  and  finite  total  energy  E  given  by 

E  -  j"  |x(t)|2dt,  (2-1-5) 

—O© 

the  effective  duration  At  is  defined  by 
At  -  ("y-  j"  (t-t)2  |x(t)|2  dt 

— OO  “ 

while  the  effective  bandwidth  Af  is  given  by 


(2-1-6) 


Af  - 


r 


(f-?)2 


|X(f)|2  df 


where 


(2-1-7) 


t 


t 


|x(t)|2  dt 


and 


f 


|X(f)|2  df . 


(2-1-8) 


(2-1-9) 


The  quantities  t  and  f  are  the  assumed  time  and  frequency  "center- 
points,"  respectively.  Using  the  Schwarz  inequality,  Gabor  showed 
that  the  product  AtAf  is  always  larger  than  or  equal  to  %  for  arbitrary 
functions.  However,  only  when  the  "elementary  signal"  is  substituted 


for  x(t)  in  (2-1-6)  and  (2-1-7)  is  it  found  that  t  ■  t  ,  f  ■  f  ,  and 

o  o 

AtAf  ■  h  exactly.  This  property  evidently  suggested  that  the  elementary 

signal  would  make  a  good  basis  for  building  a  two-dimensional  time- 

frequency  analysis.  Gabor's  development  built  upon  this  fact  as 

follows.  First,  the  elementary  signal,  ip  (t;  t  ,  f  ),  was  called  the 

e  oo 

(0,  0)th  element  and  the  generalized  elementary  signal,  ip^  (t ;  t^,  f^) , 
was  called  the  (k,  £)th  element.  (Notation  will  be  simplified  here  by 
writing  (^(t)  instead  of  ^  (t;  t^,  f^).)  It  is  seen  that  for  ^}(t), 
t  ■  t^»  f  ■  f^,  and  AtAf  *  as  before.  Next,  Gabor  defined  a  two- 
dimensional  grid  called  the  logon  grid  by  partitioning  the  time-frequency 
plane  into  rectangular  cells,  the  cell  centered  at  (t^,  f£)  being 
associated  with  the  (k,  A)**1  elementary  signal.  Each  cell  has  sides 
At  and  Af,  found  from  (2-1-6)  and  (2-1-7)  with  ^  (t)  substituted  for 

x(t) .  Specifically, 

At  -  /T~  _1_  (2-1-10) 

/  2  a 

and 

The  area  of  each  cell  is  therefore  AtAf  *  The  (k,  £)th  cell,  along 
with  its  associated  weight,  a^£,  is  called  the  (k,  A)th  logon  and  it 
corresponds  to  the  weighted  elementary  signal  given  by 


*kl*kl(t>  *  *  *a2<t'tk)2.3(2*<tt  +  s)  ‘ 


(2-1-12) 


FIG.  2-e  THE  LOGON  GR\D 

In  essence  then,  the  logon  grid  lends  itself  to  approximating  an  arbitrary 
signal  x(t)  by  the  approximation  &(t),  where 

*<t)  -  l  l  ak^kJ.(t)  .  (2-1-13) 

k*-«°  5.“-°° 

It  is  important  to  note  that  i|>  (t)  is  not  truncated  in  time  and  frequency 

as  is  its  associated  logon  at  t  *  t^  and  f  *  f  . 

It  would  be  expected  that,  if  fc(t)  is  a  "good"  approximation  in 
some  sense,  then  this  grid  of  logons  should  give  a  good  indication  of 
where  significant  amounts  of  energy  lie  in  the  time-frequency  plane  as 
the  signal  unfolds  in  time.  It  will  be  shown  later  that  this  is  not  the 
case;  in  fact  the  approximation  of  a  signal  by  the  sum  of  elementary 
signals  in  (2-1-13)  actually  obscures  the  contribution  made  by  the  signal 
in  a  finite  length  subinterval  to  the  time-variant  distribution  of  signal 
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energy.  Before  the  fallacies  of  the  logon  concept  are  shown,  the  problem 
of  making  x(t)  a  "good"  approximation  is  first  discussed. 


2.1.1.  Completeness  of  the  Set  of  Elementary  Signals 

The  function  &(t)  should  meet  certain  conditions  to  be  considered 
a  good  approximation.  The  usual  criterion  is  that  the  integrated 
squared  representation  error,  e,  be  minimized.  By  definition. 


e 


x(t)  -  *(t)|2  dt  . 


(2-1-14) 


Substituting  the  expression  for  R(t)  from  (2-1-13),  it  is  found  that 


e  -  1  |x(t)|2  dt  -  2 RE 

J 

“  00 
i  i  j  **ct>*M(t)dt 

—CO 

k*-®  x—00 

♦  i  i  I  !  vl|  <w(t>C  <*>*.  (2-i-i5) 

k«—  £— «  m—  n— «  1. 

where  the  order  of  integration  and  summation  has  been  interchanged 
assuming  the  series  to  be  uniformly  convergent.  In  order  to  minimize 
e  it  is  necessary  to  choose  the  appropriately.  This  means  that  e 
is  a  function  of  the  complex  variable  a^  and  optimization  theory  can 
be  applied  to  minimize  e  [{a,^}]  .  It  is  found  immediately  that, 
because  the  set  of  ^(t)  functions  are  not  orthogonal,  serious  diffi¬ 
culties  arise.  The  minimization  process  results  in  a  large  number  of 
simultaneous  linear  equations  for  the  a^ .  In  general,  this  system  of 
equations  cannot  be  readily  solved. 


Interestingly  enough,  Gabor  proposed  an  approach  for  determining 
a  set  of  coefficients  a^  that  does  not  explicitly  account  for  the  rep¬ 
resentation  error.  He  determined  a  set  of  by  using  a  technique  of 
successive  approximations  over  finite  intervals  of  x(t).  The  approach 
starts  by  considering  a  representation  of  x(t)  about  the  point  in  an 
interval  of  width  At.  Within  this  interval,  it  was  assumed,  as  a  first 
approximation,  that  there  were  no  contributions  from  neighboring 
elementary  signals  outside  this  interval..  In  effect,  the  approximation 
of  (2-1-13)  was  expressed  as 


*<*>■  l  l  W<c) 

k«-°° 

j6  V  ?  -a2(t-t)22,f 

-  •  I  I V e  e  1 

<x>  £«— oo 

-  f  r2(t^>2  i 


and  the  first  approximation,  (t) ,  was  defined  as 


<«>  ■  <je  «'a2<‘V  f  • 

lm-a> 


~-<t<tN+  (N-0,±  1,  ±2,...). 


(2-1-16) 


(2-1-17) 


(2-1-18) 


(2-1-19) 


It  is  noted  that  the  elements  in  the  sum  of  (2-1-18)  for  which  k  4  N 
are  ignored  in  the  first  approximation  x^  (t)  in  accordance  with  Gabor's 
assumption.  Next,  recognizing  that  he  needed  to  deal  with  only  the 
function  (t)  e+°2(t_tN)2  when  t.  -  <  t  <  +  ~~  ,  Gabor  made 


>v. 


the  approximation 


*Ct)  .  tx  (t)  «+a2(t'CN)2 


*19  tL  v  'J2’f‘t- 


CN  -  *t<c<\+  *r  ■ 


(2-1-20) 


Therefore,  an  approximate  Fourier  Series  representation  for 

x(t)  e+a^t"tll^2  was  obtained  in  the  specified  interval  centered  at  t  . 

N 

The  constants  a^  could  then  be  found  by  the  usual  procedure  for 
determining  Fourier  coefficients.  At  this  point,  Gabor  had  an  exact 

|  r«xj  .  \2 

representation  for  (t)  e  N  over  the  specified  interval  under 

the  assumption  of  no  overlap  from  neighboring  elementary  signals. 
However,  it  is  realized  that  when  the  coefficients  a  ^  are  obtained  in 
this  way,  the  representation  for  x(t)  would  in  fact  have  error  due  to 
overlap  from  neighboring  elementary  signals.  To  deal  with  this  problem 
Gabor  presented  a  heuristic  argument  for  iteratively  processing  the 
error  so  as  to  converge  on  a  set  of  coefficients  a^£  which  would  yield 
a  close  approximation.  However,  he  did  not  justify  in  any  detail  the 
convergence  of  the  iterative  process  nor  did  he  discuss  its  represen¬ 
tation  error. 

It  was  not  until  1980  that  the  open  literature  provided  a  key  to 
the  situation.  Martin  Baatiaans  [ 2] ,  using  a  set  of  functions  which 
were  found  to  be  "bi-orthogonal"  to  the  Gabor  elementary  signals,  was 
able  to  find  a  set  of  a^  for  Gabor's  expansion.  But  he,  too,  failed 
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to  discuss  the  representation  error  for  his  a^-  In  the  following 
development  it  is  shown  that  Bastiaans '  choice  of  gives  the  best  set 
of  coefficients  in  the  sense  of  minimizing  e  [Ta^f],  and  in  addition, 
establishes  the  completeness  of  Gabor's  set  of  elementary  signals. 
Completeness  is  important  because  it  guarantees  that  the  representation 
error  converges  to  zero  for  a  sufficiently  large  number  of  basis  func¬ 
tions  in  the  expansion.  It  is  completeness  that  allows  us  to  control 
the  representation  error  in  any  given  case  and  to  achieve  whatever 
accuracy  is  desired  for  the  approximation. 

Bastiaans  first  changed  Gabor's  notation,  defining  a  function 
8^4 (t)  in  place  of  </'k£(t).  Specifically, 

8^(0  "  <  a  )J*  e  7  (t_kT)2  ejT~  *  e~jlTk*  (2-1-21) 


is  used  instead  of  ik.(t),  where 


i  -a(t-t.  )2  j2irf.t  j0 

<»k£(t)  ■  e  k  eJ  l  eJ 


(2-1-22) 


To  make  the  transition  in  notation,  it  is  noted  that  Bastiaans 
substitutes 


kT  for 


(2-1-23) 


In  addition,  ha  sets  T  ■  *o  and  multiplies  his  Gaussian  modulated 

2  k 

clsoidal  function  by  a  factor  of  (— ) .  This  is  done  so  that 


r  i*kt 


(t>r  dt  - 1. 


(2-1-24) 


Therefore  each  elementary  signal  has  unit  energy. 

Bastlaan's  significant  findings  are: 

A.  There  exists  a  discrete  set  of  functions,  which  are 

"bi-orthogonal"  to  the  g^£(t)  and  exhibit  a  closure  relation  with 
gv.(t).  Specifically, 


f*<a(t)We>  Jt  ’  !»k  s„t  > 


(2-1-25) 


where  3^  is  the  Kronecker  delta  function  defined  by 


r  1  ,  m-k 

0  ,  m  i  k  , 


(2-1-26) 


I  I  gfct  (tj^)  Yj^  (tj)  ■  ^(tj^-tj) 
!(■««  Zm— ® 


(2-1-27) 


where  5(t)  is  the  Dirac  delta  function  defined  by 


x(t)  5(t-t  )dt  -  x(t  ). 

o  o 


(2-1-28) 


f  f 


B.  The  YkJl(t)  are  the  Gaussian -modulated  clsoldal  functions 
given  by 


(bk  (J-)?  a  ("kT)i,fV,UJ(-l)*  e-<»  + W2  . 


Yki(t) 


o'  NK  '  e 
o 


m>  (  y  -  k  -  »s) 


(2-1-29) 


where  the  constant  K  -  1.85407468...  is  the  complete  elliptic  integral 


for  the  modulus 


/r  • 

C.  The  Gabor  coefficients  a^  are  given  by 


r 


v  ■  }  x(t)yki  (t)  dt  » 


(2-1-30) 


where  x(t)  is  the  signal  to  be  expanded  into  the  discrete  set  of 

In  his  paper,  Bastiaans  does  not  explain  in  what  sense  this 
choice  of  coefficients  results  in  a  "good"  representation.  Consequently, 
this  set  of  a^  is  now  examined  to  determine  the  representation  error 
achieved. 

First,  g^(t)  is  written  in  place  of  \pkA(t)  in  the  expression 
for  the  error  function  e  [{a^}]  in  (2-1-15)  to  obtain 

e  “  j  1 2  dt  “  2RE  I  f  x*(t)gkA(t)dt} 

-OS  k ,  -<*> 

*im  \i *1  j <2-i*3i) 

k,l,m,n--«*  -• 


Again  it  la  noted  that  the  g^^(t)  are  not  orthogonal  aa  was  also  the 
case  for  the  ^(t).  Bastlaans'  expression  for  the  a^  In  (2-1-30)  Is 
now  substituted  Into  e  [  {a^}]  to  obtain 


e  ■  r i  *  dt 

-  2KE  {J  l  Jj 

x(t2)Ykl(t2)x*(tl)gk£(tl)  dtl  dt2J 

k,  4— •  -<*> 

J 

+HII  | | |"x^t2)Ykl^t2^x*^tl^Ymn(tl^8kl(t3)8nn(t3)dtldt2dt 

k,  l  ,m,n“-®  -• 


(2-1-32) 


Using  Bastlaans'  expression  In  (2-1-27),  there  results 


e  [{aia>]  "  f  |x(t)|2dt  -  2  RE  •If-  2>6(t1"t2^x*^ti^dt  dt2^ 

-0  —OB 

+  111  X(t2*  ^  8ki(t3)Yki<t2^][^  ^  8mn(t3)Ymn(tl|] 

k,  £■-«  tn,n»-«* 

i(tx)  dt1dt2dt3  (2-1-33) 

■  j" |x(t)|2  dt  -  2RE  { j"x(t2)x*(t2)dt2> 

—  OD  —OO 

♦  [f  u,>  6 (t--t,)6 (t,-t. )x*(t. )dt.dt«dt 


(2-1-34) 


A 


(2-1-35) 


-  j  lx(t)|2  dt  -  2RE  {  J  | x(t) j 2  dt}  +  j  | x(t) [ 2  dt 


0! 


(2-1-36) 


This  shows  that  the  representation  error  Is  zero  for  Bastlaans*  choice 
of  {a^}  provided  that  an  Infinite  number  of  functions  g^Ct)  are  used 
in  the  expansion.  Since  e  [{a^}]  is  by  definition  nonnegative,  a 
true  minimum  for  the  representation  error  has  been  obtained.  In 
addition,  this  shows  that  the  set  of  functions  g^(t)  are  complete 
because 


e  C<V>]  •  o 


(2-1-37) 


for  an  infinite  number  of  coefficients  in  the  doubly-indexed  set  of 
{ak£}* 

It  should  also  be  pointed  out  that  if  x(t)  is  known  only  over 
a  finite  representation  Interval  [ t^, tg] ,  the  Gabor  elementary  func¬ 
tions  are  still  complete  in  terms  of  the  Bastlaans  development  over 
this  finite  Interval.  To  see  this,  a^  is  modified  by  defining  the 
coefficient  a^,  where 


x(t)Yw(t)  dt. 


(2-1-38) 


The  representation  error  is  also  modified  by  defining  c'Ua^JJ, 


(2-1-39) 


1  T*  J  "  Jt  l*(t)  -  II  ^  g^COl2  dt. 

A  k,£— 

Substituting  (2-1-38)  into  (2-1-39),  it  follows  that 

e'Ka^J]  -  j‘B|x(t)|2dt 

-2RE  {ll  j|  x(t1)Yk4(t1)x  (t2)8ltA(t2)dtidt2} 
k,i—  CA 

*Ull  HI'  X<tl)l'k*(tl)l'<t2)V(t2)8kl(t3)C(t3)dtldt2dt 
k,t,*,n— •  A 

(2-1-40) 

By  inserting  (2-1-27)  into  this  expression,  there  results 

t'  -  0  ,  (2-1-41) 

for  the  infinite  set  of  coefficients  {a^}. 

This  discussion  shows  that  the  Gabor  representation  is  complete 
over  the  finite  Interval  as  well  as  the  infinite  interval,  with  the  only 
requirement  that 

J* |x(t)|2  dt  <  •  •  (2-1-42) 

Physically  real  signals  are  therefore  representable  in  terms  of  the 
Gabor  elementary  signals.  It  should  be  noted  that  this  representation 


is  valid  only  for  Bastiaans'  choice  of  the  parameters  in  (2-1-23). 


Specifically, 


-  kT  (2-1-43) 

and 

ft  -  Y  (2-1-44) 

are  noted.  These  conditions  give  g,^(t)  as  defined  in  (2-1-21). 

2.1.2.  The  Fallacy  of  the  Logon  Concept 

One  of  the  benefits  expected  from  the  Gabor  representation  for 
x(t)  was  that  the  rectangular  cell,  or  logon,  associated  with  its 
elementary  signal,  would  tell  where  significant  amounts  of  energy  are 
distributed  in  the  time-frequency  plane.  This  concept  is  now 
investigated. 

It  is  Important  to  note  that  the  representation  of  a  signal's 
energy  spectral  density,  as  the  signal  changes  in  time  within  the 
observation  interval,  is  of  interest.  One  attempt  at  analyzing  this  is 
to  compare  the  energy  spectrum  corresponding  to  one  subinterval  during 
the  observation  with  the  energy  spectra  of  neighboring  subintervals. 
Hopefully  this  will  give  a  picture  of  how  the  signal's  energy  re¬ 
distributes  itself  along  the  frequency  axis  as  variations  unfold  within 
the  time  domain  signal.  Gabor's  logon  grid  is  expected  to  provide  a 
convenient  description  of  this  energy  distribution  over  time  and 
frequency. 
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In  order  to  study  time-variant  energy,  a  finite  length  obser¬ 


vation  Interval  is  defined  to  be  [  t^,  t0].  As  shown  previously,  the 
Gabor  elementary  signals  are  complete  over  this  Interval.  Within  this 
Interval,  it  is  desired  to  analyze  the  progression  of  energy  spectra 
for  a  number  of  subintervals;  for  example,  the  three  subintervals  shown 
in  Figure  2-3. 

Re  C  '*003 


FIG.  2-3  SUBINTERVALS  WITHIN  THE  OBSERVATION 
The  total  energy  within  the  subinterval  centered  at  t  ■  NT  is 


dt. 


Using  Parseval's  Energy  Theorem, 


T 

NT  +  ! 


T 

NT-~ 


|x(t)|2  dt 


-M 


T 

NT  +  j 


x(t)  e'J2irftdt|2  df. 


T 

NT- j 


(2-1-45) 


(2-1-46) 


where  the  quantity 


T 


,<t) 


j  T 
NT  - 


(2-1-47) 


is  recognized  as  an  energy  spectral  density  and  will  hencefc  th  be 


denoted  by  |X(NT,  f,  T)|2.  Specifically, 


X(NT  ,  f,  T) 


T 

r+2 


x(t)  SiUtt  dt. 


j  x 

NT- y 


(2-1-48) 


It  Is  hoped  that,  whatever  the  representation  for  x(t) ,  the  resulting 
form  for  |x(NT,  f,  T)|2  should  clearly  show  the  distribution  of  energy 
In  frequency  due  solely  to  the  sighal's  variations  within  the  sub¬ 
interval  centered  at  NT.  For  example,  if  the  subinterval  under  con¬ 
sideration  contained  a  pure  sinisoid  of  frequency  10  Hertz,  we  would 
hope  that  the  expression  for  |X(NT,  f,  T)|2  would  have  a  significant 
magnitude  primarily  near  f  ■  10  Hertz. 

Now  |x(NT,  f,  T) | 2  is  evaluated  for  the  case  where  x(t)  is 
represented  by  x(t)  over  [tA>  tg] .  In  particular 


*(t)  "  J’J’  Vi  gk£(t)’ 

k, £■-« 


(2-1-49) 


where 


B  * 
*(t)Yv?(t)  dt, 

** 


(2-1-50) 


and 


,  \  J-  (t-kT)2 
■  (d>  e  • 


j*k£ 


(2-1-51) 


The  function  g^(t)  Is  recognized  as  the  Gabor  elementary  signal  from 
(2-1-21)  and  y.ff(t)  is  as  stated  in  (2-1-29).  It  is  recalled  from 


Section  2.1.1.  that  these  functions  are  complete  over  [tA,  t_]. 

A  B 


Substituting  £(t)  for  x(t)  in  (2-1-48),  there  restilts 

T 

NT  +  f 

X(NT,  f,  T)  -  x(t)  e~J2lTft  dt  (2-1-52) 


NT  +  y 

'll  |  8,^(0  e"j2irft  dt,  (2-1-53) 

k,t— «  NT  -  y 


where  Interchange  of  summation  and  integration  was  assumed  permissible. 
Similar  to  the  definition  of  X(NT,  f,  T) ,  the  function  G^(NT,  f,  T)  is 
defined  as 


GkA(NT,  f,  T) 


e-j2irft 


dt  , 


(2-1-54) 


which  is  recognized  as  an  Integral  over  a  finite  section  of  a  Gauss ian- 
enveloped  function.  Inspection  of  (2-1-54)  reveals  that  Gk^(NT,  f,  T) 
becomes  smaller  and  smaller  as  the  quantity  |k-N|  becomes  larger  and 
larger  because,  for  |k|»N,  the  integration  limits  lie  very  far  out  on 
the  "tail"  of  a  Gausslan-enveloped  cl sol dal  function.  Using  (2-1-54) 
in  (2-1-53),  it  follows  that 


X(NT,  f,  T)  -  ll  a££  GW(NT,  f,  T) ,  (2-1-55) 

k, £--« 

wh*te  a^'  is  the  weighting  associated  with  the  logon  centered  at  t  »  kT 

l 

and  f  •  —  .  It  will  now  be  shown  that  the  energy  spectral  density 


|X(NT,  f,  T) J 2  is  not  clearly  identified  by  the  logons  along  the  line 
t  *  NT  as  had  been  hoped.  These  logons  are  shown  in  Figure  2-4. 
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FIG.  2- A  LOGONS  ALONG  THE  LINE  t*  NT 

WITH  THE  ASSOCIATED  WEIGHTS 

Using  (2-1-55) ,  the  energy  spectral  density  becomes 


|X(NT,  f ,  T)  | 2  -  |  H^l  GkJl(NT,  f,  T)|2  (2-1-56) 

k,i— ® 

-  I  l  \\i\Z  |Gk£(NT,  f,  T)  | 2 
It,  i**-00 

*1111  <v(NT,  f,  I)  G^(MT,  f,  T). 

k,i,m,n«-«  (2-1-57) 

(kj*m  or  iAn) 

It  is  important  to  notice  that  G^^OTT,  f,  T) ,  as  defined  in  (2-1-54),  is 
not  related  to  the  signal  x(t)  whereas  a^  is  dependent  upon  the  signal 
throughout  the  entire  observation  interval  as  indicated  in  (2-1-50) . 


Hopefully,  |Gk£(NT,  f ,  T) |  will  be  small  for  k  +  N  so  that  the  major 
contributions  to  |x(NT,  f,  T)|2  can  be  identified  as  just  the  logons 
along  the  line  t  *  NT  in  Figure  2-4.  But  this  is  not  the  case  in 
general  as  shown  next. 

Substituting  the  form  for  g^£(t)  from  (2-1-51)  into  (2-1-54) 
there  results 


Gk£(NT,  f,  T)  -  <%)k  e"jm  j 


m+2  -^t-kT)2  4l.t 

7rZ  '  3  T  -j2irf t 

e  e  e  dt 


T 

NT- j 


(2-1-58) 


(l)h  ej*k£e-j2*fkT  J 


(N-k)T  +  f  _1  2 

e  °  e_j2Tr(f  “  T)ydy,  (2-1-59) 


(N-k)T--| 


where  the  change  of  variable  y  *  t  -  kT  was  made.  For  simplicity,  let 


M 


f  -  Y  wliere  M  is  an  integer.  Specifically, 


(N-^T+f  .  2  .2. 


G,  a  (NT,  _ 
ki  *  T 


M  .  .2;  jirk£  "ay  -jT(M-Oy 

,  T)  "  U4  e  J  t  e  e  dy< 


(N-k)T-y 


(2-1-60) 


M 


Setting  f  «  -  means  that  attention  is  being  focused  on  the  function 


|  X(NT ,  f,  T) | 2 1 
f 


-  |X(NT ,  T)|2 

M  1 

T 


(2-1-61) 


I  11  aUGk£(NT*f’T>l2 


(2-1-62) 


which  is  the  energy  density  at  the  specific  frequency  —  due  to  the  signal 

within  the  time  interval  of  width  T  centered  at  NT.  By  substituting 

(2-1-60)  into  (2-1-62) ,  it  is  noticed  that  all  the  logons  along  the  line 

t  ■  NT  for  which  £  i*  M  will  influence  the  calculation  of  energy  spectral 
M 

density  at  f  ■  -  and,  in  addition,  the  logons  away  from  the  line  t  -  NT 
(for  which  k  4  N)  will  also  contribute  to  the  spectral  density.  The 
latter  contribution  is  more  troublesome  because  it  indicates  that  logons 
along  the  lines  t  ■  (N  ±  1)T,  (N  ±  2)T,...,  etc.  "leak"  their  associated 
weights  aN"+1  ^  ,  aN"+2  ^...etc.  into  the  measurement  of  spectral  density 
of  the  (N,  M)th  cell  in  the  time-frequency  plane.  To  study  the  magni¬ 
tude  of  this  "leakage,"  consider  Figure  2-5  in  which  several  logons 

M 

along  the  line  f  -  —  are  shaded.  The  cell  of  interest  is  the  (N,M) 
cell  as  shown. 


(N,  CELL  IS  CELL  OF  INTEREST 
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FIG.  2-5  NEIGHBORING  LOGONS  ALONG  i »  $ 

M 

For  the  cells  along  the  line  f  *  —  ,  the  frequency  index  £*  M  while  the 
time  index  k  »  N,  Nil,  N±2 .  etc.  With  £  ■  M,  (2-1-60)  becomes 


G,  w(NT ,  5,  T)  -  (kh  e  jlTkM  f 


(N-k)T  +  y 

kM  2  y2  dy 


(2-1-63) 


which  Is  a  real  quantity  because  e 
error  function  erf(t),  defined  as 


(-1)  .  Using  the  Gaussian 


er£(t)  ■  m  \  •*’ /2  iz' 


and  making  the  change  of  variable  z  * 


(2-1-64) 


/E~ 

v  o  y,  (2-1- 


63)  can  be  written 


—  f —  X 

GkM(NT»  f*  T)  “  (2<J)U  j 


1  izl lO*  /t[(m)i+  21 


-Si 2  . 

e  dz 


rjr  t 

(N-k)  T  -  2"J 


(2-1-65) 


-  (2o)U  (-!)*“{ erf  [/  a  (N-k  +  Jj)T] 


-  erf  [/  o  (N  -  k  -  *s)TJ). 


(2-1-66) 


Recalling  that  Bastlaans  set  T  ■  Jo  ,  it  follows  that 


GkM(NT’  f*  T)  “  (20)j£  <“l>k!krf  (N  -  k  +  *5)1 


-  erf  [*^2tT  (N  -  k  -  %)]). 


(2-1-67) 


L,  kM  M 

Table  2-1  presents  the  values  of  (2o)^  (-1)  GW<(NT,  T)  in  actual 
magnitude  and  in  decibels  relative  to  the  value  at  k  »  N.  (Values  of 
erf(t)  were  found  in  Feller  [13],  page  167,  and  the  decibels  column  was 
calculated  by  the  expression  10  log  ‘ 


TABLE  B-t 

MAGNITUDES  FOR  FACTORS  ALONG  THE  UNE  -f»  *• 

(2<r)*(-0kM  GkH(NT,  $,T) 

,  ACTUM.  US  RfUTlVB  TO 

R  MACHtTUOt  k«  H  VMAJt 


N  0.7893  0.0 


N 

0.7893 

N  t  1 

O.IOSS 

N±£ 

8.642*\05’ 

N*3 

<  io"6 

<  ST8.9 


The  values  In  the  table  show  that  there  are  three  significant 
M 

— ,  T)  functions.  These  functions  are  those  for  which  k  -  N  ± 
and  k  -  N  and  they  multiply  their  associated  a^  coefficients  in 
(2-1-62).  These  associated  a^'  coefficients  correspond  to  the  logons 
shown  doubly  shaded  in  Figure  2-5.  Depending  upon  the  coefficients 
a^,  the  other  four  shaded  logons  may  also  contribute  significantly. 
This  is  a  serious  situation  because  it  runs  counter  to  our  expectation 
that  only  logons  along  the  line  t  ■  NT  should  represent  the  energy 
spectral  density  of  the  time  Interval  centered  at  NT.  Herein  lies  the 


fallacy  of  the  Gabor  logon:  what  Gabor  offered  as  a  conceptually  con¬ 
venient  tool  Is  not  convenient  for  analyzing  and  isolating  energy 
densities  from  Individual  subintervals.  The  logon  grid  misleads  one 
Into  thinking  that  there  are  distinct,  continuous  blocks  in  the  time- 
frequency  plane  where  the  signal  energy  is  concentrated.  Not  only  are 
the  Individual  energy  densities  not  concentrated  in  their  respective 
logons,  but  they  are  also  "smeared  out”  by  the  expression  In  (2-1-62) 
involving  the  logon  grid  weights  a^.  At  the  root  of  this  problem  is 
the  fact  that  the  logon  grid  is  just  not  a  legitimate  description  of 
a  signal  being  approximated  by  a  sum  of  Gausslan-enveloped  cisoldal 
basis  functions.  The  concept  of  the  logon  grid  is,  by  its  very  defini¬ 
tion,  a  two-dimensional  "map"  of  contiguous  blocks,  each  block  having 
a  weight  a^  indicating  the  strength  of  a  corresponding  basis  function 
g^(t).  It  attempts  to  localize  the  influence  of  individual  basis 
functions  to  individual  cells  of  finite  width  in  time  and  frequency. 
However,  because  the  basis  functions  are  not  strictly  limited  to  their 
individual  cells,  there  may  be  severe  leakage  effects  caused  by  large 
weighting  coefficients  a^',  thus  invalidating  the  Interpretation  of 
where  signal  energy  is  truly  concentrated.  The  logon  grid  is  therefore 
not  a  reliable  tool  for  determining  signal  energy  nor  for  representing 
the  approximation 

09 

*(t)  ■  £  I  ®k£^* 
k,£— • 


(2-1-68) 


2.1.3.  Lerner1  s  Elementary  Functions 


In  his  publication  of  1961,  Robert  Lemer  [3]  generalized  the 
Gabor  approach  to  the  logon  grid  by  extending  the  signal  representation 
to  more  general  elementary  signals.  He  allowed  the  elementary  signal 
to  have  the  general  form  v.  0(t),  where 


vki(t)  -  v(t  -  kT)e 


j4r*t 

T 


(2-1-69) 


v(t)  being  a  "convenient"  finite  energy  function  whose  energy  is  con¬ 
centrated  near  t  -  0  and  whose  energy  spectrum  |v(f)|2  is  concentrated 
near  f  ■  0.  This  is  quite  similar  to  Gabor's  elementary  signal  with 
its  Gaussian  envelope  function  multiplying  a  clsoidal  pulse.  Lerner 
approximated  a  given  signal  x(t)  with  the  function  x(t) ,  where 


*(t>-  n  \i\tM 

k,f«-« 


(2-1-70) 


Although  he  gave  arguments  to  support  his  claim  that  the  set  of  v  0(t) 


is  complete  for  practically  all  well-behaved  v(t)  with  small  time- 
bandwidth  product,  Lerner  did  not  establish  that  his  grid  of  a^'  weights 
provides  any  better  analysis  of  energy  than  Gabor's.  If,  once  again, 


T  T 

attention  is  focused  on  x(t)  in  the  subinterval  [NT  -  y  ,  NT  +  —] 


and  vk)^(t)  is  substituted  for  g^^t)  in  (2-1-53),  the  Fourier  Transform 
of  the  signal  in  that  subinterval  is 


X(NT,  f,  T)  -  H  a^  j  VkA(t) 


(2-1-71) 


k  ,4—  NT--| 


Ut  Vk£(NT*  f»  T)  h*  defll»ed»  *8  was  G.  (NT,  f,  T)  in  (2-1-54);  that  is 


T 

VkA(NT,  f,  T)  -  J  vkfc(t)  e"j2irft  dt. 

NT--| 


It  follows  that  |X(NT,  f,  T)j2  may  be  expressed  as 


|X(NT,  f,  I)|2»|[I  a^  Vk£(NT,  f,  T)|2 
k,4*-« 

“  I  I  Ia,j|2  |Vk£  (NT,  f,  T)  | 2 

k9  £■— «o 

oo 

*  HII  »kt<»T,  £,  T)  (HT,  f,  I). 

k,A,m,n«-«  . 

(ki*m  or  4i*m) 


(2-1-72) 


(2-1-73) 


Since  Lerner  did  not  consider  signal  energy  in  his  work,  he  failed  to 
require  the  energy  spectrum  of  a  finite  length  subinterval,  here  given 
by  |X(NT,  f,  T) | 2 ,  to  have  a  clear  relationship  to  the  signal  in  the 
corresponding  subinterval.  When  this  is  desirable,  the  vk£(t)  should  be 
chosen  to  yield  this  property,  yet  Lerner  neglected  this  issue. 

It  will  be  shown  in  Chapter  3  that  this  property  is  achieved  by 
defining  the  v^(t)  to  be  time-limited.  It  will  be  emphasized  there  that 
this  property  clarifies  the  relationship  between  the  individual  sub¬ 
interval  and  its  energy  spectrum.  It  is  sufficient  to  say  here  that 


Lerner  did  not  go  Into  the  Issue  of  time-variant  energy  spectra,  nor 
did  he  indicate  concern  for  the  possible  misconceptions  arising  from  a 
grid-type  of  representation  in  the  time- frequency  plane. 


!  2.2  Page’s  Instantaneous  Power  Spectrum 

In  1952,  C.  H.  Page  [4]  attempted  to  express  the  intuitive  concept 
of  a  time-variant  energy  spectrum.  He  approached  the  definition  of  an 
energy  density  distributed  in  time  andj  frequency  by  considering  the 
S  properties  which  that  energy  density  should  satisfy.  The  difficiency 

i 

of  the  ordinary  Fourier  energy  spectrum,  determined  by  the  entire  history 
of  the  signal  from  t*-°°  to  t«  00 ,  was  very  evident  to  him.  He  pointed 
out  that,  when  a  man  and  woman  sing  alternate  verses  of  a  song,  our 

I 

Intuition  predicts  that  the  energy  spectrum  should  change  with  frequency 
as  time  goes  on.  Yet  the  ordinary  Fourier  energy  spectrum  "averages  out" 
the  energy  spectra  over  consecutive,  finite-length  time  intervals. 

With  this  in  mind,  it  was  postulated  that  there  must  be  some  two- 
dimensional  function  of  time  and  frequency  which  becomes  progressively 
more  similar  to  the  ordinary  Fourier  spectrum  as  more  and  more  of  the 
signal's  history  is  included  in  the  observation  Interval. 

Page  defined  an  energy  density  p(t,f)  in  the  time-frequency  plane 
such  that  the  total  energy  in  the  signal  from  t*-*°  to  t*T  is 

T 

I  |  p(t,f)  dfdt.  (2-2-1) 


The  rate  of  change  of  total  energy  with  time  T  is  the  Instantaneous 
power  (as  a  function  of  T)  and  is  given  by 


■—  J  j  p(t,f)  dfdt  -  j  p(T,f)df 


(2-2-2) 


which  shows  that  p(t,f)  can  be  viewed  as  an  Instantaneous  power  spectrum 


at  the  time  t»T.  To  represent  p(T,f)  in  terms  of  the  given  signal  x(t). 


it  is  noted  by  Parseval's  Theorem  that  the  total  energy  up  to  time  T  is 


1 

|  |x(t)|2  dt  -  j*  |j"x(t)  e"j2lTftdt|2  df. 


(2-2-3) 


Since  (2-2-1)  expresses  the  same  total  energy  from  t**-00  to  t*T,  then 


j  j"  p(t,f)  dfdt  -  J*  || 


,f)  dfdt  •  \  if  x(t)  e"j2irftdt|2  df.  (2-2-4) 


Differentiating  both  sides  of  this  equation  with  respect  to  T  produces 


the  result 


T 

f“ p(T,f )  df  -  j~  | j  x(t)  e”^2lrftdt | 2  df.  (2-2-5) 


With  the  total  energy  being  finite,  the  integral  on  the  right  hand  side 


exists.  Furthermore,  the  partial  derivative  of  the  integrand  on  the 


right  side  is  seen  to  be  a  continuous  function  of  T  and  f.  Specifically 


T  T 

|  f  x(t)  e‘j2irftdt|2  -  2RE  {x(T)e_J2irfTj  x*(t)e+J2fff  ’ tdt> . 
*  I  — «0 


(2-2-6) 


m 


r*. 


Therefore  the  partial  derivative  may  be  taken  within  the  integral  on  the 
right  side  of  (2-2-5) ,  giving 


_3 

3T 


“  r  3T"  ,|Tx(t)  e’j2irftdtl2  df‘  (2-2-7) 


Using  (2-2-5) ,  it  follows  that 


T 

f  p(T,f)  df  -  j"  g|-  ||  x(t)  e~j2irftdt|2  df  (2-2-8) 


which  Page  considered  a  sufficient  condition  to  determine  p(T,f) .  In 
particular, 

p(T,f)  -  a|-  l|  *CO  e"J2irftdt|2 

-40 

-  2RE  (x(T)  e"J2lTfT  |  x*(t)  e+j2irf tdt}.  (2-2-9) 

_00 

Strictly  speaking,  (2-2-8)  is  an  equality  only  between  integrals  and  not, 
in  general,  between  integrands.  Nevertheless,  Page  treated  (2-2-9)  as  a 
useful  expression  for  p(T,f). 

An  Interesting  modification  of  this  concept  was  presented  by  Morris 


Levin  [5]  in  1964.  Levin  defined  the  "running  Fourier  Transforms" 


<«)  -  ]  x(t)  e 


(2-2-10) 


*j+<f)  “  |*  x(t)  e”J2lTftdt. 


Just  as  Page  Identified  the  expression 


(2-2-11) 


—  if 

3T  'J 


x(t)  e'j2irftdt|2 


(2-2-12) 


with  the  instantaneous  pover  spectrum  of  the  signal  in  the  interval 
(-«>,T),  Levin  identified  the  expression 


if  *c«>  .-J2'ftdtp 


(2-2-13) 


vlth  the  instantaneous  power  spectrum  of  the  signal  in  the  Interval 
(T,»).  Using  (2-2-10)  and  (2-2-11),  he  defined 


sf-  I  V(f> 


(2-2-14) 


Jr  lx_+(f)  1 2 


(2-2-15) 


as  the  Instantaneous  power  spectra  for  the  intervals  (-®,T)  and  (T,«), 
respectively.  Past  and  future  values  of  the  signal  x(t)  can  then  be 
treated  symmetrically  by  defining  the  new  instantaneous  power  spectrum 
p'(T,f)  as 

P  (T»f)  -  h  [jj|-  I^"(f)|2  -  gf-  |XT+(f)  I2 ]  .  (2-2-16) 

This  gives  the  mathematically  pleasing  result  that 
p "(T,f)  -  RE  (x(T)e~^27rfT  j"x*(t)  e+j52irftdt) 

.OO 

-  RE  (x(T)e"J2lTfTX*(f)}t  (2-2-17) 

where  X(f)  is  the  ordinary  Fourier  Transform  of  x(t)  taken  from  t--®  to 
t®+®.  Levin  then  showed  that  p'(T,f)  had  the  following  desirable 


properties: 


f  p'(,» 


f)  df  -  |x(T) 1 2  , 


(2-2-18) 


j”  p'(T,f)  dT  -  |x(f) | 2 


(2-2-19) 


and 


foo  ^oo  roo 

p'(T.f)  dT  df  -  |x(t)p 


dt. 


(2-2-20) 
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This  last  property  suggests  that  p'(T,f)  may  be  used  as  an  energy  density. 
All  of  these  properties,  while  mathematically  pleasing,  still  do  not 
justify  the  assumption  by  Page  that  the  Integrands  on  both  sides  of 
(2-2-8)  are  equal. 

In  order  to  put  this  instantaneous  spectrum  concept  on  a  firmer 
footing,  August  Rlhaczek  (6]  developed  a  different  approach  to  p'(T,f). 

He  considered  a  complex  signal  within  a  time  interval  AT  centered  at 


t*tQ.  Based  on  this  time  Interval,  the  total  complex  energy  E within  a 

frequency  band  AB  centered  at  f  was  defined  as 

o 


l 


V-f 


x(t)  y  (t)  dt 


(2-2-21) 


*  AT 
o  2 


where 


V1  {RMt  [fs2] x(f) 


(2-2-22) 


M-C: 


lf-f0!  <  f 

I  f-f 0 1  >  ^ 


(2-2-23) 


It  Is  noted  that  y(t)  is  the  Inverse  Fourier  Transform  of  only  that  part 
of  the  signal's  transform  within  the  frequency  band  of  interest.  Rlhaczek 
allowed  the  band  AB  to  be  chosen  so  small  that  X(f)  did  not  change  over 
it.  Specifically, 


2-34 


y(t)  -  x(fo)  '5* 


-l  (Rect  fc|2  : 
AB 


■  X(fo)  e+^irfot  AB  sine  (tAB), 


(2-2-24) 


Furthermore,  AB  was  allowed  to  become  the  Infinitesimally  small  quantity 
AB  so  that 


AB  sine  (TAB)  -  AB. 


(2-2-25) 


Using  this  in  (2-2-24) ,  there  results. 


y(t)  -  X(fQ)  AB  e+J2irf0t  . 


(2-2-26) 


Substituting  this  expression  for  y(t)  into  (2-2-21)  gives 

.  AT 
t  +  ~ 
o  2 

E  -  h  X*(f  )  AB  f  x(t)  e~i2vfol  dt. 

C  °  'AT 

o"  2 


(2-2-27) 


Rlhaczek  assumed  next  that  AT  is  so  small  that  x(t)  does  not  change  its 


value  within  it,  so  that 

«.♦¥ 

E  ■  *S  x(t  )X*(f  )  AB  f 
C  0  °  i  AT 

Co"  2 


-j2irf  t 
e  J  o  dt. 


(2-2-28) 


4  x  (t  )  X*(f  )  AB  [AT  e“^2irf0t0  ainc  (f  at)  ] ,  (2-2-29) 
o  o 


Letting  AT  become  the  infinitesimally  small  quantity  AT,  there  results 


AT  sine  (fAT)  -  6T.  (2-2-30)  ^ 

M 

Using  this  in  (2-2-29) ,  the  final  result  is 

1  -  %x(t  )  X*(f  )  e~^irfoto  5B  fix.  (2-2-31) 

C  O  O 

Inspecting  this  expression,  the  energy  density  at  the  point  (t  ,f  )  is 

o  o 

defined  to  be 

x(tQ)  X*(fQ)  e_j2irfoto  ,  (2-2-32) 


thi  product  5B6T  being  the  elementary  cell  area  in  the  time-frequency 
plane  and  the  factor  of  \  being  due  to  the  complex  notation.  £t  is 
interesting  that  this  energy  density  is  the  complex  version  of  Levin's 
expression  for  the  energy  density  (or  instantaneous  power  spectrum) 
p'(t,f)  in  (2-2-17). 

Although  these  attempts  to  define  a  useful,  intuitively  pleasing 

energy  density  over  the  time-frequency  plane  may  have  established  some 

interesting  properties,  they  fall  far  short  of  providing  a  useful  tool 

for  time-variant  energy  spectral  analysis.  The  reason  for  this  is  that 

no  satisfactory  explanation  can  be  given  for  the  fact  that  the  "energy 

density"  in  (2-2-17)  and  (2-2-32)  can  become  negative.  Rihaczek  alluded 

to  this  fact  when  he  stated  that  (2-2-32)  is  not  really  a  measure  of  the 

energy  distribution  at  (t  , f  ).  Readmitted  that  positive  values  of  his 

o  o 

energy  density  may  be  offset  by  negative  values  at  neighboring  points, 
so  that  very  little  energy  may  actually  be  concentrated  where  the 
energy  density  has  large  values.  Consequently,  it  was  concluded  that 
the  energy  density  must  be  Integrated  over  a  time-frequency  cell  chosen 


large  enough  to  prevent  significant  interaction  or  leakage  from  adjacent 
cells.  This  would  hopefully  yield  a  positive  energy  for  each  cell.  There 
is  a  similarity  here  to  the  problem  of  neighboring  cell  "leakage"  on  the 
Gabor  logon  grid.  Unfortunately,  knowledge  of  the  proper  cell  size 
requires  fine-grained  knowledge  of  the  signal's  characteristics  in  time 
and  frequency  which  are  unknown  in  many  real  problems.  Furthermore,  the 
requirement  to  furnish  the  proper  cell  size  for  integration  of  the  energy 
density  demands  some  of  the  very  information  that  one  would  use  time- 
variant  energy  spectral  analysis  to  determine! 

Another  major  difficulty  lies  in  the  Implementation  of  (2-2-17)  and 
(2-2-32).  Those  expressions  require  the  signal's  Fourier  Transform, 

X(f) ,  to  be  accurately  known.  However,  X(f)  relies  upon  knowledge  of  the 
signal  x(t)  over  the  infinite  time  interval  (-»,»).  When  x(t)  is  very 
long  it  is  not  easy  to  calculate  X(f)  with  the  precision  desired  for 
fine-grained  analysis  of  the  energy  density  of  (2-2-17)  or  (2-2-32) . 


2.3  The  Wigner  Distribution 

A  recent  discovery  [7,  8,  9]  of  some  importance  is  that  a  certain 
integral  expression  related  to  the  Ambiguity  function  has  interesting 
properties  when  used  for  time-variant  spectral  analysis.  This  integral 
expression  is  called  the  Wigner  distribution  after  the  man  who  intro¬ 
duced  it  in  the  context  of  quantum  mechanics.  The  cross  Wigner  distri¬ 
bution  of  two  signals  x(t)  and  y(t)  is  defined  by 


(t,io) 


x(t  +  h  y*(t 


-  h 


(2-3-1) 


where  u>  -  2irf  is  the  radian  frequency.  The  auto-Wigner  distribution  of 
the  signal  x(t)  is  then 


Wx  (t,u>)  -  |°°  x(t  +  j)  x*(t-  -j)  e‘j“T  dt  (2-3-2) 

which  is  similar  to  the  Ambiguity  function  Ax(t,uj)  given  by 

Ax(t,ai)  “  |~x(t  +  j)  x  (t  -  |)  e  dr.  (2-3-3) 

—<x> 

It  is  noted  that  the  roles  of  t  and  jr  are  interchanged  in  (2-3-2)  and 
(2-3-3)  but  the  interchange  only  affects  the  arguments  of  the  signal. 
It  Is  interesting  that  (2-3-2)  and  (2-3-3)  are  in  general  different 
representations.  However,  when  the  signal  is  real  and  symmetric  (even 
or  odd),  that  is, 

x(t)  -  ±  x (— t)  (2-3-4) 


then 


(t,w)  -  ±  2Ax(2t,  2oi)  .  (2-3-5) 

Some  of  the  advantages  of  using  the  auto  Wigner  distribution  (WD)  for 
signal  representation  instead  of  the  ambiguity  function  (AF)  are  the 
following. 

1.  The  WD  of  any  (real  or  complex)  signal  is  real  while  the  AF  is  in 
general  complex. 

2.  The  WD  of  a  real  signal  is  an  even  function  of  frequency  while  the 
AF  in  general  is  not. 
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3.  The  effect  of  a  time  shift  t  on  the  signal  is  to  shift  the  time 

argument  of  the  WD  by  t  while  the  effect  on  the  AF  is  to  add  a  linear 

o 

term*  -ajt  ,  to  the  phase  angle. 

0 

4.  The  effect  of  a  frequency  shift  u  on  the  signal  (or  equivalently 

o 

multiplication  by  e^a>°t)  is  to  shift  the  frequency  argument  of  the  WD 

by  u  while  the  effect  on  the  AF  is  to  add  a  linear  term,  hi  t,  to 
o  o 

the  phase  angle. 

5.  The  properties  of  a  joint  time-frequency  density  function  given  in 
(2-2-18) ,  (2-2-19)  and  (2-2-20)  are  satisfied  by  the  WD  but  not  in 
general  by  the  AF.  However,  for  real  signals,  the  AF  satisfies 
similar  relationships.  Specifically,  given  a  signal  with  total 
energy  E, 


00 

Wx  (t,a>)  dt  du  -  E,  (2-3-6) 

—00 
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2tt 


j  Wx  (V  u)  d“ 


(2-3-7) 


and 


j  W^  (t,u)Q)  dt  *  | X(ioo)  J 2 .  (2-3-8) 

••00 

On  the  other  hand,  only  for  a  real  signal  with  total  energy  E  there 
results 


00 

Ax(t,io)  dt  do»  •  2E, 


JL 

2tt 


(2-3-9) 


Ax<Vw)  d“  “  x  ("2^» 


(2-3-10) 
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and 
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A  (t,o)  )  dt 
x  *  o 


(2-3-11) 


This  shows  Chat  Che  HD  is  Che  more  convenient  choice  when  density-function 
properties  are  desirable. 

6.  The  HD  of  a  signal  tlmelimited  to  a  certain  Interval  will  itself  be 
tlmellmlted  to  the  same  Interval  whereas  the  AF  does  not  in  general 
exhibit  this  property.  Specifically,  If 


x(t)  -  0  ,  t  <  t  and  t  >  t.  , 

8  D 


(2-3-12) 


then 


H  (t.w)  ■  0  ,  t  <  t  and  t  >  t.  .  (2-3-13) 

X  SO 

7.  The  HD  of  a  signal  bandllmlted  to  a  certain  frequency  Interval  will 
itself  be  bandllmlted  to  the  same  Interval  whereas  the  AF  does  not 
in  general  exhibit  this  property.  Specifically,  if 

X(w)  ■  0  ,  at  <  u  and  u  >  w.  ,  (2-3-14) 


sa 


ri 


8.  The  instantaneous  frequency  n^(t)  of  a  complex  signal  can  be  found  by 
determining  the  average  frequency  of  the  WD.  Specifically,  the 
average  frequency  of  the  WD  is 


f  “  V* ■ 

f  Wx(t,u)) 


u>)  doj 


IH  {-j£  in  x(t)>  -  &x(t) 


(2-3-16) 


u>)  du 


This  property  is  not  exhibited  by  the  AF. 

The  preceding  properties  are  very  useful  in  describing  signal  character¬ 
istics  which  vary  in  time  and  frequency.  However,  a  crucial  point  is 
reached  when  considering  the  Wigner  distribution  as  an  energy  density. 

The  problem  is  again  that  this  "density  function"  is  not  in  general 
positive.  This  poses  great  difficulties  of  Interpretation.  For  example, 
it  may  be  asked,  "What  does  a  negative  value  for  energy  density  physically 
mean?"  or  "What  is  the  significance  of  a  negative  value  for  the  inte¬ 
gration  of  the  density  over  a  finite  region  in  the  time-frequency  plane?" 
Evidently  those  who  introduced  the  Wigner  distribution  as  a  tool  for  time- 
variant  spectral  analysis  have  not  found  answers  to  these  questions 
[7,  8,  9,  10].  It  might  be  suggested  that  the  absolute  value  or  the 
squared  value  of  the  Wigner  distribution  be  used.  However,  a  very 
significant  insight  gained  by  these  researchers  is  that  the  imposition  of 
a  positivity  requirement  on  any  tlme-fequency  energy  density  function 
prevents  it  from  satisfying  most  of  the  properties  listed  above. 


"On  the  other  hand  it  was  shown  that  positivity  of  such  a  repre¬ 
sentation  excludes  it  from  having  properties  that  are  desirable 
for  the  extraction  of  Instantaneous  information  like  the  instan¬ 
taneous  power,  Instantaneous  frequency,  and  the  finite  support 
property.  If  such  instantaneous  characteristics  should  be 
obtainable  from  such  a  time-frequency  representation,  then  we 
have  to  accept  negative  values  for  it.  This  seems  to  be  the  only 
way  to  accommodate  Heisenberg's  uncertainty  relation." [9] 

The  meaning  of  the  Heisenberg  uncertainty  relation  in  this  context  is  that 
an  arbitrarily  sharp  frequency  representation  is  Impossible  when  analyzing 
arbitrarily  short  time  intervals.  It  is  Interesting  that  the  authors  of 
the  paper  quoted  above  felt  they  were  forced  to  accept  negative  values 
for  an  energy  density.  It  will  be  shown  in  Chapter  3  that  a  different 
approach  to  time-variant  energy  avoids  this  problem. 

Another  noteworthy  fact  pointed  out  by  these  same  researchers  is 
that  the  Wigner  distribution  belongs  to  a  class  of  related  time-frequency 
representations  among  which  are  the  ambiguity  function  and  the  repre¬ 
sentation  of  Page  and  Rlhaczek.  The  relationship  between  them  is  nicely 
developed  but  here,  again,  the  interpretation  of  negative  energy  density 
is  the  main  stumbling  block  to  satisfactory  application. 

2 . A  Conventional  Discrete  Fourier  Transform  Techniques 
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N-l  . 

W£i>  '  l  1  ”  • 


where  4-0,  1,  ...»  N-l 


and  f  -  £ 
£  T 


(2-4-2) 


(2-4-3) 

(2-4-4) 


Since  the  discovery  of  computer  algorithms  for  the  rapid  calculation  of 
(2-4-2) ,  the  popularity  of  the  DFT  for  all  kinds  of  spectral  analysis 
has  grown  enormously.  In  a  typical  application,  the  signal  x(t)  Is 
first  sampled  at  an  appropriate  rate.  The  magnitude  squared  of  XjjpT(f£) 
is  then  determined  as  an  approximation  to  the  energy  spectral  density 
|XFr(f)|2  evaluated  at  the  specific  frequencies  f  -  f^.  This  result  is 
then  interpreted  as  the  sampled  version  of  the  energy  spectral  density 
for  the  signal  in  the  finite-length  observation  interval  [o,T].  However 
in  applications  to  time-variant  energy  spectral  analysis,  the  procedure 
commonly  used  is  to  partition  the  entire  observation  interval  into 
subintervals,  as  was  done  in  Figure  2-3.  For  each  subinterval  centered 
at  t^  and  having  width  T,  the  time-indexed  version  of  (2-4-1)  is  approxi¬ 
mated.  This  time- indexed  version,  XpT<f ;  t^) ,  is  defined  as 


T 

tk+  2 


XFT(f ;  tk)  -  j  x(t)  e"J2*ftdt, 


(2-4-5) 


t  -- 
*k  2 


and  it  is  approximated  by  the  time-indexed  version  of  (2-4-2), 
t^),  given  by 


n-k  +  j-  1 


XDFrCf£5  tfc)  I  H  x(tn>  e  i  4  n  , 


(2-4-6) 
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*3n| 


N  being  an  even  integer. 

The  magnitude  squared  1 2  ls  then  treflted  as  the  energy 

spectrum  for  the  subinterval  at  t^  and  the  family  of  energy  spectra  for 
a  set  of  t^  is  considered  in  common  usage  to  be  the  time-variant  energy 
spectra,  as  sketched  in  Figure  2-6. 
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RGr.  2-6  THE  FfcMtLY  OF  ENERGY  SPECTRK  f 

Although  it  ls  usual  to  interpret  this  two  dimensional  plot  as  an  energy 
density  over  time  and  frequency,  there  is  no  justification  in  the  open 
literature  for  this  interpretation.  In  particular,  there  are  no 
satisfying  answers  to  the  following  questions: 

(1)  What  is  the  relationship  between  the  energy  spectrum  of  the 
entire  signal  and  the  energy  spectra  of  individual  subintervals? 


(2)  What  interpretation  should  be  given  to  the  values  of 

I  XDFT<f  &  *  at  arbltrary  frequencies  f£? 

(3)  At  what  specific  frequencies,  f^,  should  t^)!2  be 

used  to  provide  physical  interpretation  of  a  "time-variant  energy 
spectrum? 

(4)  What  are  the  effects  of  averaging  several  sets  of  |XjjpT(f n;tk^ 
from  different  time  subintervals  t^?  That  is,  what  significance  can  be 
given  to  the  expression 

Ki 

l  IWV  VI2  ’  (2-4-10) 
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(5)  What  interpretation  should  be  given  to  the  family  of  energy 

spectra,  t^l2,  when  the  signal  x(t)  is  multiplied  by  a 

"weighting  window"  before  it  is  transformed?  (It  is  common  to  do  this 
in  order  to  attempt  better  resolution  of  sinusoids.) 

This  last  question  Introduces  the  idea  of  "window  weighting"  with  the 
window  function  w(t).  The  windowed  DFT,  (f£) ,  is  defined  as 

Wf«>  -  "f  »<V  «<'„>  e'J2rftt"  • 

n«0 

l  -  0,  1,...,  N-l, 

where  w(t)  is  usually  tapered  to  zero  at  the  edge  points  t-tQ  and  t*tN_^. 
(The  indices  l_ and  n  are  given  in  (2-4-7),  (2-4-8),  and  (2-4-9).)  A 
detailed  review  of  many  window  functions  now  in  use  is  found  in  a  paper 
by  Frederic  Harris  [12]. 


Another  common  practice  Is  to  transform  a  sequence  of  overlapping 


partitions  of  the  signal  x(t)  as  shown  in  Figure  2-7. 
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FI6.  2-7  OVER  LAP  PI  NS  PARTITIONS 

Transforms  of  overlapped  partitions  are  sometimes  used  with  the  window 
weighting  technique  to  compensate  for  the  window  function's  small  values 
near  the  edge  of  the  partition.  For  example,  if  the  signal  contains 
short  burst 8  of  sinusoids,  the  regular  non-overlap  processing  with  window 
weighting  might  de-emphasize  a  burst  which  occurs  at  a  partition  boundary. 

All  of  these  techniques  have  the  advantage  that  they  are  easily 
Implemented  on  a  computer.  However,  the  groundwork  for  interpreting 
results  is  for  the  most  part  lacking.  The  five  questions  posed  above 
point  out  the  lack  of  any  theoretical  framework  and  justification  for 
the  application  of  the  DFT  to  time-variant  energy  spectral  analysis.  The 
proposed  approach  presented  in  Chapter  3  will  be  shown  to  establish  this 
theoretical  framework  and  in  addition  answer  the  questions  posed  here. 


Chapter  3 

A  NEW  APPROACH  TO  TIME-VARIANT  ENERGY  SPECTRAL  ANALYSIS 


This  chapter  presents  a  new  concept  for  describing  the  time- 
frequency  energy  distribution  of  a  signal.  Up  until  this  point  what  was 
meant  by  a  time-frequency  energy  distribution  had  not  been  adequately 
explained.  For  our  purposes,  the  term  time-variant  energy  distribution 
will  be  taken  to  denote  a  discrete  set  of  energy  values  indexed  by 
specific  frequencies  and  specific  finite  time  intervals.  This  set  of 
values  can  be  visualized  as  the  sampled  version  of  a  two-dimensional 
nonnegative  function  of  time  and  frequency.  Hence,  the  term  time- 
frequency  energy  distribution  will  be  used  synonymously.  Negative  values 
of  energy  will  not  arise  and  the  term  energy  density  will  not  be  used. 

This  rules  out,  therefore,  the  occurrence  of  negative  energy  densities 
as  encountered  in  Chapter  2.  Secondly,  the  term  frequency,  or  more 
precisely  spectral  frequency,  will  be  used  in  conjunction  with  "shaped" 
sinusoidal  functions  of  time  having  finite  duration.  This  contrasts  with 
the  traditional  usage  which  involves  infinitely  long  sinusoids  and  it  has 
the  advantage  of  a  clearer  interpretation  when  finite-length  signal 
observations  are  being  analyzed.  Hence,  the  determination  of  a  sinusoidal 
signal's  spectral  frequency  at  a  certain  time  will  require  the  definition 
of  that  signal  over  a  finite,  nonzero-length  time  interval  which  encloses 
the  time  in  question.  It  is  clear,  therefore,  that  the  notions  of  time 
and  frequency  in  the  context  of  energy  spectral  analysis  are  not  separable. 


The  meaning  of  the  term  time-frequency  energy  distribution  as  it  will  be 


developed  here  rests  upon  the  physical  relationship  between  signal 
energy,  spectral  frequency,  and  an  associated  time  interval  as  different 
time  intervals  come  under  analysis.  This  three-fold  relati  .ship  will 
be  gradually  clarified  and  interpreted  through  the  use  of  linear  vector 
space  concepts. 

3.1.  The  Vector  Space  Description  of  Time-Variant  Enerzv  Spectral 


Analysis 

One  of  the  most  desirable  qualities  sought  in  a  time-frequency 
energy  distribution  is  the  clear  indication  of  progressive  changes  of 
the  frequency  distribution  of  energy  from  one  moment  to  the  next  within 
a  signal's  observation  interval.  It  is  therefore  natural  to  break  up  an 
observed  signal  into  contiguous  subintervals,  each  of  which  is  analyzed 
sequentially.  For  our  purposes,  the  signal  observation  interval  will 
be  denoted  by  [tA,t^]  while  each  subinterval  will  be  defined  as  T 
seconds  long  and  centered  at  a  particular  time  t.  as  shown  in  Figure  3-1 


c - START  IN6 

index 

FIRM.  INDEX 

FIG.  3-1 

SIGNAL  OBSERVATION 

INTERVAL 

As  Indicated  in  the  figure, 

Vi  -  *k  - T' 

where  the  indices  are  integers  given  by 

k  “  Ks*  Ks  +  1 . Kf’  (3-1-2> 

It  is  assumed  that  there  are  N  subintervals,  that  is, 

Kf  -  Kg  -  N  +  1.  (3-1-3) 

T  T 

The  portion  of  the  signal  contained  in  the  interval  [t^-  t^+  j  ]  will 

be  referred  to  as  the  kth  segment.  It  is  important  to  note  that  the 
signal  is  assumed  to  be  unknown  outside  of  [t.,t„]  and  all  attention 
will  be  given  to  what  is  going  on  within  this  total  observation  interval. 

The  proposed  approach  consists  in  representing  the  signal  in  each 
subinterval  by  a  point  in  a  linear  vector  space.  The  location  of  this 
point  and  its  motion  in  time  will  automatically  indicate  associations 
with  various  types  of  sinusoidal  functions  of  specific  frequencies  and 
a  convenient  description  of  time-variant  signal  energy  will  emerge. 
Specifically,  the  entire  observed  signal  x(t)  is  approximated  by  the 
function  £(t)  where 

Kf  M 

*<t>  ■  [  I  "b  v<e) 

k-K  t— M 


(3-1-4) 


and  where  the  doubly— indexed  orthonormal  basis  functions  ^^(t)  are  given 

by 

/w(t-tk)  ejV  ;  |t-tk!  <| 

0  5  lt_tJ  y  \  >  (3-1-5) 

w(t)  being  known  as  the  "window".  Also  the  radian  frequencies  are  given 
by 


where  Q  is  referred  to  as  the  "frequency  spacing  factor,"  a  non-zero 
integer  uniquely  associated  with  the  window  w(t).  Inspection  of  (3-1-4) 
and  (3-1-5)  shows  that  the  signal  in  each  segment  is  being  approximated 
by  a  discrete  set  of  basis  functions,  that  these  basis  functions  are 
non-overlapping  functions  of  time,  and  that  each  basis  function  is  a 
"shaped"  or  "windowed"  complex  cisoidal  function,  e^<i)£t,  characterized 
by  its  particular  radian  frequency 

It  is  quite  significant  that  the  square  root  of  the  window 
function,  /w(t) ,  is  the  shaping  tool.  An  illustration  of  one  such 
window  is  given  in  Figure  3-2  along  with  the  real  part  of  a  typical 
basis  function  Here  the  window  function  is  the  so-called 


"Bartlett"  window,  given  by 


4  r\ 


0  ;  1 1 1  >  T/2, 


with  Q“2  as  its  associated  frequency  spacing  factor. 


(3-1-7) 


F\G.  3-1  THE  BARTLETT  WINDOW,  ITS  SQUARE  ROOT, 
AND  A  T\P\CAL  BASIS  FUNCT\ON 

Since  the  basis  functions  defined  in  (3-1-5)  are  non-overlapping  and  time 
limited  to  intervals  of  width  T,  the  window  function  w(t)  will  be 
assumed  to  be  time-limited  to  [-T/2,  T/2]  without  loss  of  generality. 
Strictly  speaking,  the  window  function  by  itself  need  not  be  time-limited 
but  the  only  portion  of  it  of  any  consequence  is  the  portion  in  the 
Interval  [-T/2,  T/2].  This  will  be  seen  shortly.  The  assumption  of  a 
time-limited  window  is  consistent  with  the  basic  motivation  to  represent 
the  signal's  frequency  distribution  of  energy  as  a  separate  entity  for 
each  time  subinterval.  It  is  interesting  to  note  that  Gabor's  choice  of 
basis  functions  in  (2-1-1)  does  not  have  this  property  which  is  one  of 
the  underlying  reasons  why  the  logon  concept  is  not  useful. 

Introducing  the  Fourier  Transform  of  the  time-limited  window 
function  as  W^(u),  where 


T/2 

V  (n\\  m  f  e  ^ 


(3-1-8) 


it  is  required  for  the  orthonormality  of  the  basis  functions  that 


WT(0)  -  1 
and  that 


(3-1-10) 


(3-1-11) 


n  -  ±1,  ±2,...  .  (3-1-12) 

Of  primary  importance  here  is  the  requirement  in  (3-1-11)  which  states 
that  WT(u>)  must  have  an  infinite  number  of  equally  spaced  zeroes,  the 
spacing  being  Although  not  ruling  out  complex  window  functions,  in 

which  case  the  procedure  here  would  call  for  replacing  w (t)  by  |w  (t) | 
in  (3-1-8)  and  w(t  -  t^)  by  jw(t  -  t^)  I  in  (3-1-5),  it  is  much 
easier  in  the  computation  of  coefficients  to  be  discussed  later  to  use 
real  window  functions.  Hence  w(t)  will  be  assumed  to  be  real  here  while 
the  signal  x(t)  will  be  allowed  to  be  complex  in  general. 

It  is  interesting  to  note  in  passing  that  the  Bartlett  window  in 
(3-1-7)  satisfies  (3-1-10)  and  (3-1-11).  In  particular,  for  the 


Bartlett  window. 


Simplifying  the  integral  in  the  k  -  m  case,  it  is  found  that 


T 

\+2 


T/2 


f  v(t-t.  )  e'^VV'dt  -  e'^VV^  J 

1  T  *•  -1 


t  -- 
Ck  2 


.  .  -j (w  -w,)u. 
w(u)e  J  n  l  du. 


T/2 


(3-1-18) 


where 


u  *  t  -  t, 


(3-1-19) 


and 


du  ■  dt. 


(3-1-20) 


The  right  side  of  (3-1-18)  is  just  WT(w)  evaluated  at  the  frequency 


(w  -ui.)  and  multiplied  by  a  constant  complex  exponential.  Specifically 
n  i 


(3-1-18)  is 


rJ((vVtk  w(a)) 

i  n  x, 


e  W„ 


^  (n-i) 


(3-1-21) 


Using  (3-1-10)  and  (3-1-11),  it  follows  that 


e  T  *  W„ 


(n-i) 


In  * 


(3-1-22) 


so  that  (3-1-17)  becomes 


Sln  ;  k  '  * 


11  ’  < 


0  ;  k  i<  m 


(3-1-23) 


i 
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J-10 


w(t)  “  Rect 


T-8 


0  a(t) 


(3-1-24) 


(T-8) 
2 


a(t-x)  dx. 


(3-1-25) 


-llzR 


In  addition. 


w(o)  -  | 


(T-8) 
2 


a(x)  dx 


(3-1-26) 


S2=SL 

2 


1  T 

which  evaluates  to  when  B  <  The  Fourier  Transform  of  w(t) , 


which  la  identical  to  WT(w)  in  this  case,  is  the  product  of  the  trans 


form  of  Rect 


T-B 


WT(to)  -  W(a>)  -  A(w)  • 


and  the  transform  of  a(t) .  Specifically, 


2  sinl 


<D 


(3-1-27) 


where 


A(u) 


[  .<t>  e-Jut 


dt. 


(3-1-28) 


Clearly  W^(w)  has  zeroes  at  the  uniformly  spaced  frequencies 


u  ■  n 


2w_ 

T-B 


f 


(3-1-29) 


n  -  ±1,  ±2 


(3-1-30) 


as  well  as  at  the  zeroes  of  A(co)  as  sketched  in  Figure  3-5. 


FIS-.  3-S-  COMPOSITION  OF  w(c4  AS  A  PRODUCT 
OF  SOLID  AND  DASHED  LINES 

It  follows  from  (3-1-27)  and  from  the  (j^g)  area  property  of  a(t)  given 
in  Figure  3-4  that 


WT(0)  -  1 


(3-1-31) 


in  satisfaction  of  the  orthonormality  requirement  in  (3-1-10) .  However 
in  order  that  (3-1-29)  satisfy  the  orthonormality  requirement  in 
(3-1-11) ,  S  must  be  chosen  so  that 


(3-1-32) 


8  -  T  (^i)  , 


(3-1-33) 


where  Q  is  the  integer  "frequency  spacing  factor"  used  previously. 

Thus  the  parameter  S,  which  controls  the  time  duration  of  a(t)  and  the 
area  under  its  curve,  can  be  appropriately  chosen  to  ensure  the  ortho¬ 
normality  of  the  basis  set.  The  main  benefit  therefore  from  using 


able  parameter  6.  Specifically,  these  time  domain  requirements  are 


8/2 

|  a(t)  dt  -  ^  (3-1-34) 

-8/2 


and 


6  “  T(^)  .  (3-1-35) 

In  actual  practice,  once  values  for  the  segment  width  T  and  the 

frequency  spacing  Q  have  been  given,  8  can  then  be  determined  from 

(3-1-35)  and  a  suitable  symmetric  positive  function  a(t)  can  be  chosen 

to  have  area  and  duration  8.  The  specific  shape  of  a(t)  is  a  use- 
r— p 

ful  design  tool  and  will  be  discussed  later. 

As  an  example  of  a  member  of  the  class,  when  Q»l,  (3-1-35) 
requires  8  to  vanish.  Taking  the  limit  of  both  sides  of  (3-1-34)  as 
8  -*■  0  gives  the  requirement  that 


8/2 

lim  f  a(t)  dt  -  i 

8  -»■  0  j  T 

-6/2 


(3-1-36) 


which  can  be  satisfied  by  setting 


a(t)  -  ±  6(t), 


which  is  the  Dirac  Delta  function  with  area  —  . 


(3-1-37) 


Substituting  (3-1-37)  into  (3-1-24),  there  results 


w(t)  -  i  Rect  I 


> 


(3-1-38) 


the  so-called  "Rectangular"  window.  It  is  interesting  to  note  that  by 
setting  Q-l,  the  smallest  possible  "frequency  spacing"  of  the  basis 
function  in  (3-1-5)  is  achieved.  This  spacing  is  ±~,  i^-,  ...etc. 

As  another  example,  when  Q-2,  (3—1—35)  requires  that  B  - 
Examining  (3-1-34)  for  this  case,  it  is  found  that 


T/4 

* 

a(t)  dt 

T/4 


2 

T  * 


(3-1-39) 


There  are  many  symmetric  positive  functions  a(t)  that  will  satisfy  this 
requirement.  One  such  function  is 


a(t)  -  —  Rect 


(3-1-40) 


sketched  in  Figure  3-6 . 
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FIG.  3-6  THE  FUNCTION  a(t)  CHOSEN  IN  (3-1-40) 

When  (3-1-40)  is  substituted  into  (3-1-25)  there  results 
T/4 

w(t)  -  J  a(t  -  t)  dr  (3-1-41) 

-T/4 


Fifi.  3-10  the  “raised  cosine"  window  or  (  3-i-4<V) 


3-17 

This  window  is  known  as  the  "Raised  Cosine"  or  the  "Hanning"  window.  It 
is  interesting  to  note  that  when  these  windows  are  used,  Q  ■  2  and  the 
frequency  spacing  of  the  basis  functions  in  (3-1-5)  becomes  ±^-  ,  , 


The  frequency  spacing  factor  is  an  interesting  quantity.  By 
inspection  of  (3-1-5)  and  (3-1-6) ,  it  is  evident  that  the  basis  functions 
are  shaped  cisoldal  functions  with  radian  frequency  oi^.  The  expression 
for  a)^  in  (3-1-6)  indicates  that  only  integer  multiples  of  the 
fundamental  radian  frequency  are  allowable  frequencies  for  the  shaped 
cisoldal  functions.  Thus  the  larger  the  value  for  Q  the  larger  the 
fundamental  frequency  and  the  larger  the  difference  in  frequency  between 
each  basis  function. 

Now  that  the  orthonormality  of  the  basis  functions  in  (3-1-5)  has 

been  established  and  a  means  for  their  selection  discussed,  it  is 

appropriate  to  consider  the  goodness  of  the  approximation  in  (3-1-4). 

Clearly  the  quality  of  the  approximation  depends  upon  the  choice  of 

coefficients  c,  „  and  one  of  the  most  common  criteria  of  selection  is 
ki. 

the  integrated  minimum  mean  squared  error  criterion.  This  criterion 
requires  that  the  c,  .  be  chosen  in  such  a  way  that  the  integrated  mean 
squared  error. 


e 


j x ( t )  -  R(t) | 2  dt 


(3-1-50) 


be  minimized.  Substituting  (3-1-4)  into  (3-1-50),  there  results 


3-18 


r  B  f 

a 

l*(t)  -  I 

JtA  k-K 

l  Ck£  *ki(t)l 

£*-M  KJt 

(3-1-51) 


Kf  M 


~  2RM  l  l 


k-K  £— M  JtA 

s  A 


x  (t)  *k£  (t)  dt 


Kf  M 


+  l  l  I  ck5>  I  2 » 

1 _ ir  n  1/ 


k-K  £»- M 
s 


(3-1-52) 


where  use  was  made  of  the  orthonoraality  of  the  set  to  derive 

the  last  term  on  the  right  side  of  (3-1-52).  Since  the  signal  is 
assumed  to  be  complex  in  general,  the  expansion  coefficients  in  the 
approximation  must  also  be  complex  in  general.  Thus, 


ck£  “  ®k£  +  j  ’ 


(3-1-53) 


where  g^  is  the  real  part  and  h^  the  imaginary  part.  Making  this 
substitution  in  (3-1-52)  it  is  found  that 


|x(t) I2  dt 


r  Kf  M  fc i 

-  2RE(  l  l  (gk£  +  j 

I  k-K  l— M  't. 

^  s  / 

Kf  M 

+  l  l  +  h?t) 

k-K  l— M  *  k 


x  (t)  *M(t)  dt 


(3-1-54) 


-  2  Jk,  i-H  •»  “1  l  X  (t>^<t>  dt] 


Kf  M  r  rtB  *  l 

+  2  Jk  { Jt 


Kf  M  Kf  M 

I  I  gjjA  +  I  I  h£ 

k-K  £--M  k-K  b*-M 

8  S 


(3-1-55) 


It  is  seen  now  that  e  is  actually  a  function  of  the  two  sets  of  variables 
(®k£^  and  80  t*ie  minimization  of  a  two  dimensional  function 

£(^®k£^  .{^J)  must  be  performed.  There  are  two  criteria  for  a  true 
minimum  at  some  point  •  gQ  and  h^  *  hQ.  Specifically, 


(3-1-56) 
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5ST72 


<*o’ho> 


(*o'ho) 


3gU3\l 


<Vh„> 


(3-1-57) 


where  all  partial  derivatives  are  evaluated  at  the  point  (g  ,h  ) .  From 

o  o 

0-1-55)  it  is  seen  that,  for  a  specific  k  and  £., 
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and 


9  e 
3gk£ 

92e 

3glJ 

3e 

3hk£ 
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Imposing  the  conditions  in  (3-1-56),  it  follows  that 


*k£ 


-  RE  i 

‘VV  *■ 


B  x*(t)  t|»ka(t)  dt 
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L>Z 


and  that 


-IM< 


<Bo’ho> 


B  x*(t)  dt 


or,  equivalently,  that 


"ki. 


gk£ 


(g^»h  ) 


i  bki. 
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(grt.h  > 
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(3-1-58) 

(3-1-59) 

(3-1-60) 

(3-1-61) 

(3-1-62) 

(3-1-63) 

(3-1-64) 


(3-1-65) 


-  RE 


X  (t)ik  .(t)  dt  V  -  j  IM« 


x  dt 
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REi  [  B  x (t)  «*4(t)  dt  l  +  j  IM 


(3-1-66) 


x(t)^k£(t)  dt 


(3-1-67) 


x(t)  ^k£(t)  dt. 


(3-1-68) 


In  addition,  when  (3-1-59),  (3-1-61),  and  (3-1-62)  are  substituted  into 
condition  #2  in  (3-1-57),  it  is  found  that  the  expression  given  in 
(3-1-68)  is  a  true  minimum.  This  expression  is  commonly  known  as  the 
Fourier  coefficient  of  x(t)  relative  to  the  orthonormal  set  {^k£(t)>. 
Taking  c^  to  be  defined  in  this  way,  that  is, 


c  ■ 

ki 


*(t)  ll»k£  (t)  dt. 


(3-1-69) 


it  follows  from  (3-1-52)  that  the  minimum  error  is 


fB 

|x(t)P 
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k-K  A--M 


(3-1-70) 


(3-1-71) 


The  optimum  coefficient  in  (3-1-69)  can  also  be  written  as 


'ki 
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fck+7 

t  -  - 
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x(t)  /w(t-t^) 


-J2$  U 


dt 


(3-1-72) 


which  is  quite  significant  in  that  the  coefficient  for  a  given  segment 
centered  at  t  ■  t^  is  completely  independent  of  the  signal  outside  that 
kth  segment.  In  terms  of  the  approximation  *(t)  in  (3-1-4),  this  means 
that  the  £(t)  function  can  be  viewed  as  the  sum  of  individual  segment  approxi 
mations.  Defining  the  approximation  of  the  k*^  individual  segment  to  be 
x^t),  where 

\(t>  ■  I  Mcki  *U  (t)-  <3-1-73) 


k  -  Kgl  +  1 . Kf  ,  (3-1-74) 

it  follows  that  the  total  approximation  is 
Kf 

x(t)  -  l  *.(t).  (3-1-75) 

k-K  K 
s 

The  fact  that  each  x^t)  Is  solely  under  the  influence  of  the  signal 
x(t)  in  its  own  subinterval  provides  the  framework  within  which  the 
vector  space  approach  is  now  developed. 

It  is  recalled  that  the  proposed  approach  concerned  the  repre¬ 
sentation  of  the  signal  x(t)  in  each  subinterval  by  a  point  in  a  linear 
vector  space.  This  vector  space  is  now  specified  as  the  complex  vector 


space  of  piecewlse-continuous  square  lntegrable  functions  defined  over 
a  particular  time  subinterval.  The  basis  vectors  of  this  space  are 
the  complex  functions  ^^(t)  with  k  constant  so  that  the  factor 

/w(t-t^)  in  (3-1-5)  is  fixed  for  the  vector  space.  All  the  usual  vector 
space  properties  are  valid  here  as  well  as  the  inner  product.  Speci¬ 
fically,  for  three  vectors  £k(t),  £k(t),  and  zk(t)  3Pace  having 

Pk£.’  yki’  and  Zk£  38  their  corresponding  Fourier  coefficients  relative 
to  the  set  { ^ C t: > }  ,  the  following  conditions  are  satisfied. 

(1)  Equality  Relationships 
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*k(t) 


ykl  "  Zkl’  Pk2 


yk2  +  zk2  *  * ' 


(3-1-76) 


(2)  Sum  Relationships 
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(3-1-77) 


(3)  Multiplication  by  Complex  Scalar  Relationships.  Given  the  complex 
number  a. 


£k(t) 
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(3-1-78) 


(4)  Difference  Relationships 


(5)  Zero  Vector  or  Origin  Definition 


yk(t)  "  0  ^  ykl  *  °’  yk2  0 . 


(3-1-80) 


(6)  Inner  Product  Definition 


<  yk(t)*  2k(t)  > 
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(7)  Norm  or  Length  Definition 


II  yk(t)  ||  -  <yk(t),  yk(t)> 


(3-1-84) 
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Also  the  norm  ||  y^(t)  “  z^(t)  ||  is  called  the  distance  between  y^(t) 
and  z.  (t).  It  is  important  to  note  that  the  symbol  y.  (t)  represents 


tne  tc  segment  approximation  of  the  signal  y(t)  defined  over  [t  ,t  ] 

A  B 


and  that 


•c 


B  y(t)  *k *(t)  dt 


(3-1-87) 


represents  the  corresponding  Fourier  coefficients. 

Using  the  previously-discussed  k**1  segment  approximation  x^Xt) 
if  follows  from  (3-1-84)  and  (3-1-73)  that  the  norm  squared  of  the 
vector  x^t)  is 
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a  very  important  result.  This  states  that  the  energy  in  the  kth 
approximation,  ^(t),  as  seen  in  the  right  side  of  (3-1-88),  is  just 
the  sum  of  its  squared  components  in  the  vector  space,  or  equivalently, 
the  distance  squared  from  the  origin  to  the  vector  space  point  associ¬ 
ated  with  x^t).  This  geometric  interpretation  arises  from  the  fact 
that  the  basis  vectors,  the  ^(t),  are  orthonormal.  Now,  using 
(3-1-75),  the  energy  in  the  total  approximation,  St( t) ,  is  expressed 


(3-1-91) 


~f 

-  I 

k«K 

s 


M 


l 

M 


(3-1-92) 


(3-1-93) 
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Thus  the  energy  in  the  total  approximation  is  the  sum  of  the  squared 
magnitudes  of  all  the  individual  coefficients  c^- 

These  squared  magnitude  values,  the  |ckJJ2  are  quite  significant 
in  themselves.  They  represent  the  energy  in  each  individual  building 
block  of  the  representation,  namely,  each  weighted  basis  function 
cki  ^kfc^’  This  ^act  *8  made  evident  by  considering  the  energy  in  an 
arbitrary  building  block,  that  is. 


f B  Ki  *nt(t)|2  dt  ■  |ckJ2  LB  !*M(t>|2  dt '  |ckt|2  • 

CA  *A 

(3-1-95) 

But  more  than  this,  these  |  ckJl  |  2  values  represent  a  distribution  of 
energy  in  time  and  frequency;  in  time  by  virtue  of  the  association 
with  the  kth  segment  of  the  signal,  and  in  frequency  by  virtue  of  the 
association  with  the  2th  shaped  cisoidal  function,  /w(t-t,  )  e^w2C. 


This  is  strikingly  similar  to  the  family  of  energy  spectra  produced 
by  the  conventional  Discrete  Fourier  Transform  application  discussed  in 
Section  2.4.  Here,  however,  a  theoretical  -justification  has  been  given 
for  interpreting  these  values  as  an  energy  distribution.  In  addition, 
the  relative  contributions  of  each  shaped  cisoidal  building  block  to 
the  energy  of  the  total  approximation  can  now  be  easily  compared. 

In  terms  of  the  vector  space  concept,  the  distribution  of 
energy  in  time  and  frequency  is  seen  in  the  movement  of  the  kC  segment 
approximation  of  vector  £^(0  from  one  location  to  another  as  the  index 
k  changes.  This  movement  produces  two  important  pieces  of  information: 

(1)  The  squared  distance  from  the  origin  to  the  point  associated 
with  £^(0  may  change  from  the  kt^1  time  Interval  to  the  (k+l)tl1  time 
interval,  where  the  window  function  w(t-t^)  in  the  representation 
remains  the  same  except  for  a  time  shift  of  T  seconds.  This  indicates 

a  change  in  energy  of  the  approximation  from  one  subinterval  to  the 
next. 

(2)  As  the  location  of  the  point  associated  with  x^Ct)  in  the 
vector  space  changes,  a  quantitative  indication  is  given  of  the  changing 
influence  of  each  basis  vector  4^(0  upon  the  total  approximation. 

This  influence  can  be  measured  by  the  absolute  value  of  the  cosine  of 
the  angle  0^  between  the  ^(t)  vector  and  a  particular  basis  vector  as 
illustrated  in  Figure  3-11.  When  the  approximation  vector  k^(t)  points 
in  nearly  the  same  direction  as  the  angle  0^  will  be  very 

small  and  the  point  associated  with  ^(t)  will  lie  very  near  to  the 
line  along  which  e(t)  lies.  Consequently,  the  influence  of  <|i,  .  (t)  is 


expected  to  be  relatively  large  and  this  fact  is  borne  out  by  a  measured 
influence  value  close  to  one.  This  is  also  true  when  ^(t)  points 
nearly  in  the  opposite  direction  to  as  well.  However,  when 

j^(t)  points  almost  perpendicularly  away  from  a  particular  ^(t)  basis 
vector,  the  influence  of  that  basis  vector  will  be  measured  to  be  very 
small  as  expected  since  x^Ct)  will  now  be  pointing  almost  in  the  same 

direction  as  some  other  basis  vector.  This  other  basis  vector  will 

) 

now  carry  the  greatest  influence  with 


FIG-.  3-11  INFLUENCE  OF  y^  (t)  BASIS  VECTOR 


Using  the  energy  values  previously  discussed,  the  influence  of 

a  particular  basis  vector  ij /  (t)  can  be  expressed  as 

pm 


[pm  '  •C0S<VI 


(3-1-96) 


which  is  just  the  familiar  length  ratio  of  the  side  adjacent  to  the 
hypotenuse  in  an  M  dimensional  vector  space.  However,  upon  closer 


Inspection,  it  is  recognized  that  I  is  the  magnitude  of  the  complex 

pm 


th 


correlation  coefficient  between  the  p  segment  approximation  x  (t)  and 
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the  basis  vector  1 p  (t).  This  is  seen  in  the  definition  of  the  complex 
pm 


correlation  coefficient  p  which  is 
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Thus  the  definition  of  Influence  in  (3-1-96)  has  a  valid  interpretation 
as  a  measure  of  how  strongly  related  the  approximation  vector  x(t)  is 
to  a  particular  basis  vector. 

So  far  a  vector  space  representation  has  been  described  which 
establishes  the  justification  for  interpreting  a  set  of  squared  mag¬ 
nitude  components  as  a  time- frequency  energy  distribution.  The  next 
topic  concerns  the  question  of  representation  error  for  various  window 


functions. 


3.2  Completeness  of  the  Vector  Space  Representation 


By  completeness  for  the  representation  in  shaped  cisoidal 
functions  is  meant  that  the  approximation  error  goes  to  zero  as  the 
number  of  basis  vectors  is  increased  without  ?  >und.  Using  the 
expression  for  approximation  error  in  (3-1-71)  and  noting  that  the 
integer  M  controls  the  number  of  basis  functions  in  the  representation, 
it  follows  that  the  approximation  error  in  the  limit  is 
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The  behavior  of  the  second  term  on  the  right  side  of  (3-2-2)  is  now 
Investigated.  Using  the  definition  of  the  coefficient  c^  in  (3-1-69), 
there  results 
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Now  the  fact  that 
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can  be  used  to  determine  the  limiting  form  of  the  right  side  of  (3-2-2). 
Assuming  that  the  summation  and  the  integration  operations  can  be  inter¬ 
changed,  it  is  found  that 
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Noting  that  w(t-t^)  is  timelimited  to  the  interval  t.  <  t  <  t, 
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and  that  w(t  -  tk — — )  is  timelimited  to  the  interval 
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tk  +  |  (•“  -  i>  <.  t  <  tk  +  |  +  1)  , 
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it  follows  that  the  product  w(t-tk)v(t  -  t^  - ~)  is  zero  whenever  | A |  >  Q. 
Therefore  the  infinite  summation  in  (3-2-10)  becomes  a  finite  summation. 
Specifically, 
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This  expression  can  be  reduced  further  by  separating  the  £®0  term  from 
the  other  terms  in  the  summation.  Specifically,  the  right  side  of 
(3-2-13)  is 
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irfiere  the  first  term  results  from  1-0,  the  first  braketed  term  results 
when  £*±1,  and  finally  the  last  bracketed  term  results  when  £»±(Q-1) . 
The  general  term  above  is 
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It  Is  now  shown  that  this  general  term  consists  of  just  a  quantity  plus 
Its  complex  conjugate.  It  Is  noted  Immediately  that  the  first  quantity 
within  the  bracket  of  (3-2-15)  has  the  product  w(t-tk)w(t-tk--^)  under 
the  square  root.  The  two  factors  In  this  product  are  sketched  In 
Figure  3-12  for  Jt  a  positive  Integer. 
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It  is  seen  from  the  sketch  that  the  region  of  overlap  is  given  by 
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so  that  the  first  quantity  within  the  bracket  can  have  its  integration 
limits  changed.  Specifically,  the  first  quantity  is  more  properly 
expressed  as 
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Similarly,  the  second  quantity  within  the  bracket  of  (3-2-15)  can  be 

expressed  as 
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l  (  T  x(t)x*(t  +  ~)  i/w(t-t.  )w(t-t. +~)  dt,  (3-2-18) 


where  Figure  3-13  provides  a  sketch  of  the  factors  in  the  product 
AT 

w(t-t^)w(t-t^+-Q-)  and  the  corresponding  region  of  overlap  which  justi¬ 
fies  the  change  of  integration  limits. 
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Nov  the  change  of  variable  x  -  t  +  can  be  made  in  (3-2-18)  with  the 
result  that  the  second  quantity  in  (3-2-15)  is  just 
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which  is  recognized  as  the  complex  conjugate  of  the  first  quantity 
within  the  bracket  as  expressed  in  (3-2-17).  Therefore  the  general 
term  in  (3-2-15)  is  written  as 
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and  the  result  Incorporated  into  (3-2-14)  and  (3-2-13)  to  give 
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Substituting  this  back  into  the  limiting  form  of  the  approximation  error 
in  (3-2-2),  it  follows  that 


In  general  this  approximation  error  reduces  to  zero  only  for  the 
Rectangular  window  given  in  (3-1-38)  for  which  Q  was  shown  to  be  unity. 
Indeed,  when  Q«l,  the  second  term  on  the  right  side  of  (3-2-22)  is  non¬ 
existent  while  the  first  term  evaluates  to  zero  for  the  Rectangular 
window.  Thus  the  representation  error  goes  to  zero  and  the  set  of 
basis  functions  is  complete  for  the  Rectangular  window.  This  holds  for 
arbitrary,  finite  energy  signals  x(t).  This  important  result  guarantees 
that  whatever  representation  accuracy  is  desired  may  be  achieved  by 
using  the  Rectangular-shaped  cl sol dal  functions.  This  result  is  intu¬ 
itively  expected  once  the  Rectangular  window  case  essentially  reduces 
to  the  ordinary  complex  Fourier  Series  expansion  applied  separately  to 
each  tlmellmlted  segment  of  the  signal.  If,  for  a  finite  term  approxi¬ 
mation,  greater  accuracy  of  representation  is  desired,  all  that  needs  to 
be  done  is  to  Increase  the  number  of  basis  functions  and  calculate  the 
corresponding  new  coefficients  using  the  Fourier  coefficient  expression 
In  (3-1-69).  The  only  coefficients  that  need  to  be  calculated  are 
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those  associated  with  the  newly  added  basis  functions.  Since  the  set 
of  basis  functions  is  complete ,  this  procedure  guarantees  a  smaller 
error  of  approximation  each  time  it  is  applied.  Thus  the  rectangular 
window  expansion  is  an  excellent  tool  for  signal  representation. 

It  is  interesting  to  investigate  the  limiting  form  of  the 
approximation  error  in  (3-2-22)  for  non-rectangular  windows  w(t).  First, 
however,  a  convenient  form  for  the  signal  x(t)  will  be  written.  Intro¬ 
ducing  the  Rectangular  window  representation  for  the  signal  x(t)  as 
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where  the  Rectangular-windowed  coefficients  r^£  are  given  by 
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it  is  possible  to  substitute  this  form  for  the  signal  x(t)  into  any 
expression  involving  the  signal  since  this  representation  is  complete 
as  discussed  above.  Accordingly,  for  the  purposes  of  this  discussion, 
(3-2-23)  will  be  formally  substituted  into  (3-2-22)  merely  to  facili¬ 
tate  a  useful  form  for  the  non-rectangular  windowed  approximation  error. 
This  does  not  imply  a  contradiction  since  the  non-rectangular  windowed 
approximation  of  x(t)  is  just  being  analyzed  with  the  particular  x(t) 
given  by  (3-2-23).  As  will  be  seen,  this  enables  a  tractable  expression 
for  the  approximation  error  in  the  non-rectangular  case  to  be  developed. 
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Substituting  (3-2-23)  into  (3-2-22),  the  following  simplification 
can  be  made  immediately.  Since  the  integrals  in  (3-2-22)  are  concerned 
with  values  of  t  in  the  Interval 
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and  since,  from  (3-2-23), 
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in  this  Interval,  (3-2-26)  can  be  used  in  (3-2-22) .  There  results 
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The  first  term  on  the  right  side  above  is  more  conveniently  expressed 

as 
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vfatrc  Interchange  of  summations  and  Integration  operations  was  assumed 
possible.  Now  the  integral  in  the  first  triple  suamation  of  (3-2-28) 
la  just 

tk+2 

T  a  dt  -  T  sine  (i-i)  -  T6±  £  ,  (3-2-29) 

St  "2 


so  that  (3-2-28)  becomes 
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(3-2-30) 


Making  the  change  of  variable  x  -  t  -  t^,  dx  -  dt,  and  noting  that  w(x) 
has  even  syametry  about  x  ■  0  by  virtue  of  its  composition  in  (3-1-24) 
by  convolution  of  two  even  functions,  (3-2-30)  yields 
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i  rkirti  * 


#  »(t)  co.j^Y  U-»)t1  di.  (3-2-31) 


The  first  ten  on  the  right  side  of  (3-2-31)  is  the  total  energy  in 
the  signal  as  expected  by  careful  review  of  the  first  ten  in  (3-2-22) . 
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The  second  term  offsets  the  first  term  by  subtracting  a  sum  of  weighted 

* 

cross  products  of  the  Rectangular-windowed  coefficients  The 

second  term  on  the  right  side  of  (3-2-27)  can  be  simplified  by  making 
the  'variable  change 


*  -  A  - 1  ’  dt 


(3-2-32) 


which  gives 
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(3-2-33) 


where  Interchange  of  summation  and  integration  was  again  assumed  possible 

IT  IT 

Now  the  product  w(t  +  :^)w(t  -  2q)  can  be  shown  to  have  even  symmetry 
about  t -  0.  In  particular,  since  w(t)  -  w(-t)  for  any  t,  then 


v(t+t)  *  w(-r  -T> 


(3-2-34) 
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w(-T  +  T  ) 


(3-2-35) 


w(t  +  J  )  w(t  -  |)  -  w(— T  -  J  )  w(-t  +  |) 


■  w(-T  +  J  )  v(— T  -  J  )  , 


(3-2-36) 


(3-2-37) 


where  Che  right  side  of  (3-2-37)  is  the  sane  as  the  right  side  of 
(3-2-36)  except  that  the  factors  are  re-ordered.  But  the  right  side  of 


(3-2-37)  can  be  expressed  as 


v(-t  +  -  )  v(-r  -  J  )  -  w(t  +  J  )v(t 


-1)! 
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J  t“-T 


,  (3-2-38) 


thus  showing  that 


w(x  +  j  )  w(x 


w(t  +  J  )  w(t 


-  }  )  (3-2-3 
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which  proves  the  even  symmetry  about  t  ■  0.  This  permits  the  right 
side  of  (3-2-33)  to  be  expressed  as 

^f  OO  os  2ir 

-  4 1  i  i  i  «*  iv 


w(T  +|~)  w(T--g)  cos|j(p-n)Tjdt] 

(3-2-40) 


This  epuation  represents  the  second  tern  on  the  right  side  of  (3-2-27) 
end  it  can  be  combined  with  (3-2-31)  which  represents  the  first  term  on 
the  right  side  of  (3-2-27)  to  yield  a  tractable  form  for  the  approxi¬ 
mation  error.  Before  this  is  done,  the  indices  ±  and  in  the  triple 
summation  in  (3-2-31)  are  re-labeled  as  £  and  n,  respectively,  while 
the  Indices  in  (3-2-40)  remain  unchanged.  In  addition,  since  (3-2-31) 
is  real  by  virtue  of  its  derivation  from  the  first  term  on  the  right 
side  of  (3-2-27)  it  can  he  brought  within  the  Real  Operator,  RE{.}. 

This  gives  the  approximation  error  as 
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w(t  +|^)w(t  --~)cos  J^(p-n)xJ  dt 


(3-2-41) 


At  this  point  it  is  noticed  that  the  triple  summation  over  k,  n,  and  £ 
above  can  be  separated  into  two  kinds  of  summations;  one  for  which  n  -  p 
and  one  for  which  n  i  p.  Performing  this  separation  gives  the  final 


form  for  the  approximation  error  as 
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(3-2-42)“ 

As  an  illustrative  example,  the  sinusoidal  "burst"  signal  is  considered. 


Specifically , 


x(t)  -< 


A  cos  t  +  d>  ] ;  |t-tR|  <  | 


(3-2-43) 


0  ;  ,t-tK'  > 2  * 


where  D  is  a  positive  integer,  A  and  4>  are  arbitrary  finite  amplitude 
and  phase,  respectively,  and  is  a  specific  point  in  time.  A  sketch 
is  given  in  Figure  3-14. 


F\&.  3-14  SINUSOIDAL  *  BURST*  CENTERED  AT  tK 

Assuming  that  tR  is  the  centerpoint  of  a  particular  subinterval  in  the 
Rectangular  windowed  expansion  of  x(t) ,  the  coefficients  r.  in  the 
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Rectangular  window  expansion  are  determined  from  (3-2-24)  to  be  r  -  0 


for  k  +  K  and 
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(3-2-47) 


tdiere  T  •  t  -  tR,  dx  -  dt.  Integration  yields 
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It  is  noted  that  “  0  for  k  i4  K  since  the  signal  is  confined  to  the 
K**1  subinterval.  In  terms  of  the  vector  space  concept,  the  signal  is 
represented  by  a  vector  in  a  two  dimensional  space  with  equal  length 
projections  along  the  two  basis  vectors 


jfNt 


(3-2-50) 


The  energy  in  each  component  of  the  signal  vector  is 


rk,±N 
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(3-2-51) 


and 


lrk,J2  "  0  for  k  *  K’  1  *  ±N* 


(3-2-52) 


Substituting  (3-2-49)  into  (3-2-42) ,  the  approximation  error  for  the 
window  w(t)  with  its  associated  frequency  spacing  factor  Q  is  found  to 
be 


w(t)  cos(ytfT)  dr 


Ycofl<TNvq 
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(3-2-53) 


This  is  an  interesting  result  because  it  gives  a  tractable  expression 
for  the  approximation  error  for  an  arbitrary  window  provided  that  the 
signal  is  the  sinusoidal  "burst"  of  (3-2-43).  If  the  frequency  of  this 
sinusoid. 


f  -  f  ,  (3-2-54) 

were  not  an  integer  multiple  of  ^  ,  then  the  approximation  error  would 
be  a  complicated  function  of  many  coefficients  c^  but  (3-2-42)  could 
still  be  used  to  calculate  it. 

When  the  Bartlett  window,  given  in  (3-1-7),  is  used  in  (3-2-53) 
along  with  the  associated  Q  value  of  2,  the  approximation  error  is 
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(3-2-58) 


The  second  bracketed  tern  in  (3-2-55)  will  in  general  be  a  function  of 
the  value  of  the  integer  N  and  the  centerpoint  of  the  signal,  tg. 
However  the  first  Integral  in  this  second  bracketed  tern  is  evaluated 
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The  second  Integral  within  the  second  bracketed  tern  of  (3-2-55)  is 
plotted  in  Figure  3-15  "he  specific  integral  being  calculated  for 
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Fl€r.  5-\5  VALUES  OF  INTEGRAL  m  (Vl-4^  FOR 

integer  N 


It  la  noted  that  the  dashed  curve  la  given  for  comparison  purposes  only 
In  particular,  for  N  -  10,  the  absolute  value  of  the  Integral  Is  less 

_3 

than  1.3  x  10  T.  Thus  the  second  bracketed  term  In  (3-2-55)  makes  a 
relatively  small  offset  to  the  first  bracketed  term  In  (3-2-55) .  As 
shown  In  Figure  3-15,  the  maximum  value  of  the  Integral,  36.9  x  10  T, 
occurs  when  M  ■  1,  Combining  this  result  with  the  value  of  the  first 
bracketed  term  In  (3-2-55)  as  expressed  In  (3-2-58),  there  follows,  for 


it 


Jllae- Y  [*  +  f~  0-0738  cos(y<k)J 
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^integer' 
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l£  the  percentage  approximation  error,  Pg  ,  is  defined  as  the  actual 

approximation  error  in  (3-2-55)  divided  by  the  total  signal  energy  which 
2 

is  from  (3-2-51),  then  for  N  -  1, 


81. 9%  <  V  <  96. 7% 

e 


(3-2-65) 


Similarly,  when  N  -  2,  combining  the  graphical  value  from  Figure  3-15 
of  -13.8  x  10~3X  with  the  first  bracketed  term  from  (3-2-58)  gives  the 
bounds  on  the  percentage  approximation  error  as 


7.9%  <  <  13.5%. 

e 


(3-2-66) 


Table  3-1  presents  the  bounds  on  the  percentage  approximation  error  for 
odd  and  even  values  of  the  integer  N  up  to  14.  It  is  clear  from  the 
tabular  data  as  well  as  the  decaying  value  of  the  Integral  in  Figure 
3-14  that  P£  approaches  the  limiting  form  in  (3-2-58)  of 


Pe  -  ft  -  (-D"  f)  x  100% 


10.7%,  N  even 


89.3%,  N  odd 


(3-2-67) 
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*8  N  gets  larger  and  larger. 

TABLE,  3-1 

BOUNDS  ON  PERCENTAGE  ERROR,  P6>WITH  BARTLETT 

WINDOW 
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Since  (3-2-58)  is  the  dominant  term  which  gives  this  limiting  form, 
the  term  in  (3-2-55)  from  which  it  was  derived  will  be  called  the 
primary  sinusoidal  error  term.  Specifically,  the  primary  sinusoidal 
error  term,  rip»  is  defined  from  the  left  side  of  (3-2-55)  as 

T/2  T/4 _ 

rip  -  (1  -  j  w(T)dT  -  2(-l)N  |  ^(T  +  |)w(T-  j)dx}.  (3-2-61 

0  0 


A  / 

-1)N  j  ^(T  +  |-)w(T-|)dT}.  (3-2-68) 


Similarly,  the  secondary  sinusoidal  error  term,  n  >  is  defined  as  the 

8 

second  bracketed  term  on  the  left  side  of  (3-2-55),  that  is. 
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w(t)  cos(^-Nt)  dT 


(3-2-69) 


+  2  j  /w(t  +  ±)w(t-±)  C0«(^T)dT). 

0 

It  is  noted  that  these  error  terms  pertain  only  to  the  representation  of 
a  sinusoidal  burst  by  non-rectangular  windows  with  associated  frequency 
spacing  factors  of  Q  ■  2. 

It  is  interesting  to  consider  why  the  percentage  approximation 
error  in  (3-2-67)  should  depend  so  critically  upon  the  evenness  or 
oddness  of  N.  It  is  recalled  from  (3-1-6)  that  the  radian  frequencies 
in  a  representation  with  Q  -  2  are 
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(3-2-70) 


But  the  radian  frequency  of  the  sinusoid  under  analysis  here  is,  from 

2ffN 

(3-2-43) ,  just  -y  •  Thus  when  N  is  an  even  Integer,  the  frequency  of 
the  sinusoidal  burst  is  contained  in  the  set  of  frequencies  while  for 
odd  Integer  N  it  is  not.  It  should  be  no  surprise  then  that  the  per¬ 
centage  approximation  error  is  so  high  in  (3-2-67)  for  N  odd.  This 
effect  will  always  occur  when  the  actual  signal  contains  sinusoidal 
components  whose  frequencies  are  not  contained  in  the  set  of  frequencies 
of  the  clsoidal  basis  functions. 

As  another  example  of  non-rectangular  window  approximation 
error,  the  Hanning  window  given  in  (3-1-49)  is  considered.  Again,  Q  ■  2 
and  the  primary  sinusoidal  error  term  is  found  by  substituting  (3-1-49) 
into  (3-2-68)  which  gives 
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The  quantity  under  the  square  root  contains  the  two  terms 
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-  2  coaC^T)  cos(|)  -  0  ,  (3-2-73) 

where  use  was  made  of  the  fact  that,  for  any  9  and  ♦  , 

cos  (0  +  $)  +  cos  (0  -♦)  ■  2  cos  0  cos  (3-2-74) 

In  addition,  the  product  term  under  the  square  root  Is 

CO.  [f  <t+£>]co.[£ 

.  h  [cos(^)  +  cos(iT)j  ,  (3-2-75) 

where  use  was  made  again  of  the  trigonometric  Identity  In  (3-2-74) . 
Substituting  these  results  into  (3-2-72) ,  it  follows  that 
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Now,  for  any  6, 
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(3-2-77) 


so  that  the  integral  which  involves  the  square  root  in  (3-2-76)  is 
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Since  the  left  side  of  (3-2-55)  can  be  written  in  terms  of  and  ng 
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the  Hanning  windowed  approximation  error  is  found  by  substituting 
(3-2-82)  and  (3-2-93)  into  (3-2-94).  This  gives 
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(3-2-95) 


(3-2-96) 

The  contribution  from  n  g,  which  is  the  factor  multiplying  the  cos(^y  tR) 

-3 

term,  is  smaller  than  1.3  x  10  for  N  >  8.  Dividing  (3-2-96)  by  the 

a2t 

total  signal  energy,  -y  ,  and  multiplying  by  1002,  the  percentage 
approximation  error,  P£  ,  is  found  for  N  >  8  to  be 


18.22;  N  even 
81. 8Z;  N  odd 


(3-2-97) 


where  use  was  made  of  (3-2-83). 

Once  again  it  is  seen  that  the  percentage  approximation  error  is 
quite  large  for  odd  N,  as  expected.  Comparing  (3-2-97)  with  (3-2-67),  a 
lower  error  is  obtained  for  the  odd  N  case  while  a  higher  error  is 
obtained  for  the  even  N  case  by  the  Hanning  window  over  the  Bartlett 
window.  This  suggests  that  a  trade-off  may  be  made  by  using  certain 


3-58 


windows;  one  may  be  willing  to  suffer  a  higher  approximation  error  for 
representing  sinusoids  with  frequencies  in  the  set  of  represented  fre¬ 
quencies  (N  even)  so  that  a  lower  representation  error  may  be  achieved 
for  sinusoids  with  frequencies  outside  the  represented  set  (N  odd) .  Or 
one  may  find  that  the  converse  of  this  trade-off  is  desirable.  In  any 
case,  the  evaluation  of  n_  and  q  provides  a  useful  figure-of-merit  by 

I*  3 

which  representation  error  may  be  compared  for  various  non-rectangular 
windows . 

At  this  point  it  may  be  asked,  "Why  should  a  non-rectangular 
windowed  expansion  be  used  in  the  first  place  if  a  non-zero  signal 
representation  error  will  always  result?"  The  answer  is  that  although 
it  is  true  that  a  non-rectangular  windowed  basis  set  is  not  complete  in 
general  and  will  not  generally  give  zero  signal  representation  error, 
energy  spectral  analysis  is  not  primarily  concerned  with  signal  repre¬ 
sentation  but  with  spectral  representation.  For  example,  a  10%  or  even 
an  89%  signal  representation  error  may  be  tolerable  in  order  to  get  a 
more  sharply  peaked  and  more  steeply  graded  appearance  to  the  distribution 
of  energy  values  over  frequency.  These  properties  greatly  facilitiate 
the  process  of  estimating  the  true  frequency  of  a  sinusoidal  burst,  which 
is  of  prime  Importance  in  spectral  analysis.  As  will  be  seen  in  the 
following  examples,  the  windowed  expansions  enjoy  the  benefits  of  sharp 
spectral  peaks  and  rapid  decay  from  those  peaks.  These  advantages 
counterbalance  the  undesirable  aspects  of  poor  signal  representation  error . 

The  first  example  is  the  plot  in  Figure  3-16  of  the  Bartlett 
windowed  expansion  coefficients  for  the  case  of  the  sinusoidal  burst 
given  by 


I 


K  2 


(3-2-98) 


Specifically,  the  plot  shows  the  normalized  values  of  the  energies 


rK  2  r 


T 

t_,  +  —  /  '  |  i  ■  ,  4tt£ 

^  /TISI  cos[HI60i  t)  e‘3  T  dt|t> 


(3-2-99) 


which  is  just  the  squared  magnitude  of  the  expression  for  the  coefficients 
in  (3-1-69)  with  k  -  K  and  the  Bartlett  "windowed"  basis  functions 


Af. 


2  t-t. 


|t-tKi  <f 


(3-2-100) 


0  5  lt“tKl>2 


The  values  of  |c,r.|2  are  normalized  so  that  the  maximum  value  is 
assigned  a  zero  decibel  level.  Each  value  is  scaled  by  the  rule  10 
login  |c  |2  .  It  is  noted  that  the  maximum  energy  value  occurs  as 
desired  at  the  frequency  index  of 


2£max  "  60 


which  corresponds  to  the  basis  set  radian  frequency  of 


(3-2-101) 


MAX  .  2tt(60) 


(3-2-102) 


Thus  there  are  60  complete  cycles  of  the  sinusoidal  burst  In  the  K 
segment  and  the  frequency  of  the  sinusoid  Is  represented  In  the  set  of 


basis  functions  (N  Is  even) .  In  addition  the  centerpoint  of  the  K 


th 


segment,  t  ,  was  chosen  so  that  the  sinusoidal  burst  had  either  a 

Iv 


maximum  or  a  minimum  at  that  point.  It  is  worthwhile  to  note  that  the 
second  largest  energy  values  In  the  plot  are  smaller  than  the  maximum 
by  about  28  dB  and  that  the  decay  rate  is  about  7.8  dB  between  frequen¬ 
cies  10  and  20  Hz  away  from  the  center  frequency  of  60  Hz.  The  influence 


values  1^  are  also  plotted  for  each  frequency  in  the  representation  in 


Figure  3-17.  The  scale  for  this  plot  is  the  logarithmic  rule  10  log 


10 


I^£  and  therefore  the  influence  values  have  the  same  shape  in  Figure 


3-17  as  the  energy  values  in  Figure  3-16.  This  is  seen  in  the  fact 
that 


10  logio  hi  "  10  108 


10 


(3-2-103) 


M 


5  *og  lckJ2  "  5  £°8  <  l  lckp ! 2> 

p«-M  v 


(3-2-104) 


which  is  the  same  as  the  plotted  energy  values,  10  log1Q  }  ckJlI 2 »  except 
for  the  different  multiplication  factor  of  5  and  the  constant  offset 
term,  neither  of  which  affect  the  relative  shape  of  the  plot.  Therefore 
the  shapes  of  the  energy  value  plot  and  the  influence  value  plot  will 
always  be  the  same.  It  is  noted  that  the  smooth  curves  are  drawn 
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Fig.  3-17  Bartlett  Influence  Values  for  60  Cycle  Sinusoid 


between  the  points  only  for  ease  of  comparison  and  are  not  to  imply  a 
continuum  of  values  for  arbitrary  frequencies. 


A  second  example  is  given  in  Figure  3-18  in  which  the  frequency 
of  a  sinusoidal  burst  is  not  contained  in  the  set  of  represented  fre¬ 
quencies  of  the  Bartlett  windowed  expansion.  Specifically,  the  signal 
is  given  by 


x(t)  »< 


■—  1^1  ;  It-g  <  f 


«•  I'-'kI>7 


(3-2-105) 


and  the  plotted  values  are  almost  the  same  as  those  from  (3-2-99) 
except  that  the  radian  frequency  of  the  cosine  is  Instead  of 

deg,,  from  the  plot  that  there  are  two  maxima,  each  of 
which  is  approximately  15.5  dB  above  the  next  lower  values  and  that  the 
decay  rate  is  about  9.2  dB  between  frequencies  of  10  and  20  Hz  away 
from  the  most  significant  frequencies  of  60  and  62  Hz.  Essentially  the 
Bartlett  windowed  expansion  represents  this  61 -cycle  sinusoidal  burst 


with  two  equally  strong  shaped  sinusoids  having  60  and  62  cycles  within 
the  T  second  segment.  Here  is  a  case  where  the  "Bartlett  window"  per¬ 
centage  approximation  error  in  (3-2-67)  is  very  near  to  89. 3%  and  yet 
the  spectral  representation  is  still  quite  good.  A  more  precise  des¬ 
cription  of  spectral  representation  is  given  by  the  influence  value 
plot  of  Figure  3-19.  There  it  is  seen  that  the  basis  functions  with 


m 


60  and  62  cycles  within  the  a  segment  have  the  equal  influence  of 
0.490792  upon  the  approximation  signal  (taken  from  the  maximum  value 
data  in  the  block  at  bottom  of  the  figure) .  But  more  significant  is 
the  fact  that  the  next  lower  values  of  influence  are  more  than  7  dB  down 


from  these 


.  This  means  that  the  signal  is  being  approximated  by 


essentially  two  equally  weighted  basis  vectors.  It  is  also  important 
to  note  the  high  degree  of  synmetry  about  the  "correct"  frequency  of 
61  even  though  it  is  not  represented  in  the  set  of  plotted  points. 

This  gives  the  distinct  impression  that  a  sinusoid  with  61  cycles  in 
the  segment  is  the  cause  of  these  effects.  The  Influence  values  plotted 
for  the  60  cycle  case  in  Figure  3-17  also  indicate  a  sharp,  central 
maximum  value  with  next  lower  values  down  from  the  maximum  by  about 
13  dB.  In  each  case,  whether  the  frequency  is  within  the  set  or  not, 
the  influence  plots  show  central  maxima  which  are  clearly  dominant  over 
all  other  values  and  in  addition  have  sharp  symmetry  properties  which 
strongly  imply  the  correct  frequency  of  the  signal.  These  are  certainly 
desirable  spectral  characteristics. 

Now  it  should  be  pointed  out  that  the  rectangular  windowed 
expansion  is  superior  in  its  spectral  representation  to  all  non- 
rect angular  windowed  expansions  when  the  frequency  of  the  sinusoidal 
burst  under  consideration  is  in  the  set  of  represented  frequencies. 


This  is  evident  from  (3-2-44)  through  (3-2-52)  which  shows  that ,  for 
Integer  frequency  index  N  (that  is,  N  cycles  in  the  segment),  the  set 
of  rectangular  windowed  coefficients  consists  of  Kronecker  Delta  func 
tlons  located  at  i  ■  ±N,  the  exact  frequencies  of  the  sinusoid.  Thus 


Che  indication  of  the  true  frequency  of  the  signal  is  unmistakable  for 
this  case.  However,  when  the  frequency  of  the  sinusoid  is  not  repre¬ 
sented  in  the  set  of  basis  functions,  the  spectral  representation  by 
the  rectangular  windowed  expansion  is  not  as  "sharply  defined"  as  a 
non-rectangular  expansion  such  as  the  Bartlett.  By  "sharply  defined" 
is  meant  the  two  properties  of  having  a  narrow  width  about  the  mum 
value  and  also  having  a  rapid  decay  rate  or  strong  slope  away  from  the 
peak.  As  an  Illustration,  Figures  3-20  and  3-21  show  the  rectangular 
windowed  expansion  coefficients  and  associated  Influence  values  for  a 
60.5-cycle  sinusoidal  burst  in  the  K*”*1  segment.  Specifically,  the 
signal  is 

'  cos  (60.5)t]  ;  It-tJ  S  \ 

x(t)  -  (3-2-106) 

„  °  J  I *“*1^  >  2 

It  is  interesting  that  the  energy  values  as  well  as  the  influence 
values  for  these  rectangular  windowed  coefficients  are  nicely  centered 
in  magnitude  about  the  correct  frequency  index  of  60.5.  It  should  be 
seen  here  that  the  index  ranges  over  all  integers  as  expected  from  the 
rectangular  windowed  basis  functions.  The  second  largest  values  of 
energy  each  have  magnitudes  of  about  9.0  dB  smaller  than  the  maxima  at 
l  m  60  and  i  -  61,  and  the  energy  decay  rate  between  frequencies  10  and 
20  Hz  away  from  the  central  "point"  of  60.5  Hz  on  the  horizontal  axis  is 
about  5.7  dB.  Similarly,  the  influence  value  plot  in  Figure  3-21  has  its 
imnHma  at  1*60  and  £*61,  about  4.4  dB  above  the  next  smaller  values.  The 


decay  rate  is  about  2.7  dB  In  the  same  region  of  frequencies  measured 
before.  This  is  not  as  sharply  defined  as  the  Bartlett  windowed  spec¬ 
tral  representation  for  the  same  signal.  This  representation  is  given 
in  Figures  3-22  and  3-23.  The  energy  values  in  Figure  3-22  show  a  single 
maximum  at  &»60  although  there  is  a  slight  skewness  to  the  plot  indi¬ 
cated  by  the  single  second  largest  values  at  62  Hz.  This  seems  to  imply 
that  there  may  be  a  strong  contribution  to  the  signal's  frequency  at 
some  idex  value  to  the  right  of  60  Hz.  In  addition,  the  maximum  value 
is  more  than  20  dB  larger  than  the  second  largest  value,  in  contrast  to 
the  9.0  dB  value  found  for  the  rectangular  windowed  coefficients  in 
Figure  3-20.  The  decay  rate  for  the  Bartlett  windowed  expansion  is  also 
greater;  about  9.3  dB  of  decay  between  frequencies  10  and  20  Hz  away 
from  the  60.5  point  as  compared  to  the  rectangular  windowed  value  of  5.7 
dB.  The  most  significant  result  here  is  that  the  influence  values  for 
the  Bartlett  windowed  expansion  show  a  single  dominant  basis  vector  of 
60  Hz  frequency  with  Influence  of  about  10  db  above  the  next  significant 
basis  vector  with  62  Hz  frequency.  It  is  the  magnitude  of  influence 
which  is  desirable  here.  On  the  other  hand,  the  rectangular  expansion 
influence  values  have  their  maxima  only  4.4  dB  above  the  next  significant 
influence  values.  This  shows  that  the  Bartlett  windowed  expansion  is 
more  sharply  defined  in  its  spectral  representation  and,  although  it 
suffers  from  a  larger  spacing  between  its  frequencies,  its  strong  sym¬ 
metry  and  skewness  can  be  used  to  infer  the  presence  of  other  frequencies 
outside  the  represented  set.  Thus  the  spectral  representations  of 
windowed  expansions  like  the  Bartlett  are  still  quite  good  even  though 
the  associated  approximation  error  may  be  large. 


■MM  MUK  ■  .1«1M  MUM  NUK  • 

m mum ms  » mmm  bnv  *.  j an 
■nica  mw  a  a  it  t.um 


•  W.S 
4.42NK-43 


j  •  xo 


Fig.  3-20  Energy  Values  for  Rectangular  Windowed  Expansion 
for  60.5  Cycle  Sinusoid 
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Fig.  3-22  Energy  Values  for  Bartlett  Expansion  of  60.5 

Cycle  Sinusoid 
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Bartlett  Expansion 
Sinusoid 


It  is  worthwhile  to  note  that  the  spectral  characteristics  of 
windowed  expansions  stay  be  adjusted  by  changing  the  shape  of  the  func¬ 
tion  a(t)  in  (3-1-24)  and  (3-1-25) .  It  is  recalled  that  this  function 
must  be  symmetric  about  t*0  and  have  duration  B  and  area  where  6 
is  given  by  (3-1-33).  The  specific  shape  of  a(t)  is  arbitrary  within 
these  limits.  As  a  rule,  the  broader  the  "shoulders"  of  a(t),  the 
narrower  the  appearance  of  the  maximum  value  of  its  Fourier  Trasnform, 

A  (to) .  It  is  this  frequency  function,  A(to) ,  that  plays  an  important 
part  in  the  formation  of  W(to) ,  the  Fourier  Transform  of  the  window  w(t) . 
This  is  evident  in  Figure  3-5  which  shows  that  the  product  of  the 
adjustable  function  A(to)  and  a  fixed  function,  2  sin  [^-  (T-B)]/to, 
produces  W(u>) .  From  the  figure,  the  adjustable  function  A(to)  is  a 
slowly  varying  "envelope"  for  the  more  rapidly  varying  fixed  function. 
The  adjustable  function  can  be  constructed  to  have  a  specified  decay 
rate  and  "main  lobe"  width  so  that  W(oj)  will  follow  these  envelope 
characteristics.  Since  the  windowed  expansions  discussed  above  use  the 
square  root  of  the  window,  it  is  the  Fourier  Transform  of  /w(t)  that 
produces  the  spectral  representation  properties.  In  any  case,  the 
spectral  characteristics  desired  for  a  windowed  expansion  can  be 
achieved  by  judicious  choice  of  the  shape  of  the  symmetric  function  of 
time,  a(t) .  For  most  of  the  commonly  used  windows,  the  analytical 
expression  for  the  Fourier  Transform  of  A/(t)  does  not  exist  so  that 
a(t)  in  (3-1-25)  or  A(u)  in  (3-1-27)  must  be  adjusted  experimentally 
by  computer  analysis  to  accomplish  the  narrow  width  and  rapid  decay 
specifications  desired. 


Discrete  Fourier  Transform  Techniques 


This  section  concerns  the  five  questions  posed  in  Section  2.4. 
Although  the  answers  are  implicitly  connected  with  the  preceding  parts 
of  this  chapter,  a  detailed  treatment  of  them  is  given  here.  Each 
question  will  be  stated  follows  by  its  answer. 

(1)  "What  is  the  relationship  between  the  energy  spectrum  of 
the  entire  signal  and  the  energy  spectra  of  individual  subintervals?" 
The  Fourier  energy  spectrum  of  a  signal  x(t)  defined  over  an  entire 
observation  interval  [t^,  tfi]  is  the  squared  magnitude  |x(u>)|2  where 


x(t)  e  dt 


(3-3-1) 


(3-3-2) 


Using  the  proposed  segmentation  scheme  of  Figure  3-1,  there  results 
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(3-3-3) 


“  j  U)  ( t — T  ) 


(3-3-4) 


This  shows  that  the  sum  of  the  Individual  segment  Fourier  energy  spectra 
in  the  first  term  on  the  right  side  above  is  not  in  general  equal  to  the 
energy  spectrum  of  the  entire  signal.  On  the  other  hand,  using  the  pro¬ 
posed  vector  space  approach,  it  was  shown  in  (3-1-94)  that  the  energy  in 
the  total  approximation,  x(t),  is  the  sum  of  the  energies  of  the  indi¬ 
vidual  basis  functions  in  the  expansion.  This  Important  result  allows 
the  comparison  of  contributions  from  each  segment  to  the  overall  energy 
distribution.  Moreover  the  vector  space  approach  approximates  the  sig¬ 
nal  in  a  given  segment  by  a  vector  whose  length  squared  is  the  energy  of 
the  approximation  for  that  segment.  As  time  advances  in  the  signal, 
another  vector  is  formed  to  represent  the  signal  in  the  next  segment. 
Since  both  vectors  are  members  of  the  same  vector  space,  their  relative 
orientation  Indicates  th<*  changing  influence  of  different  sinusoidal 
frequencies.  On  this  bauis,  the  frequency  distribution  of  energy  in  the 
signal  can  be  said  to  be  changing  with  time.  Thus,  the  vector  space 
approach  provides  a  clear,  useful  concept  for  time-variant  spectral 
analysis.  As  a  specific  example,  the  rectangular  windowed  expansion 
reduces  to  a  sampled  version  which  is  identical  to  the  unwlndowed 
Discrete  Fourier  Transform  approach  of  Section  2.4. 

(2)  "What  interpretation  should  be  given  to  the  values  of 
|  Xjjft 2  at  arbitrar?  frequencies  f  ?"  As  mentioned  previously 
in  this  chapter,  the  proposed  signal  expansion  is  valid  only  for  the 
discrete  radian  frequencies  in  (3-1-6) .  A  physical  significance  for 
energy  distributed  at  other  frequencies  is  not  possible  with  this 
vector  space  approach.  This  is  consistent  with  the  Heisenberg  Uncer¬ 
tainty  Principle  mentioned  in  Section  2.3  in  that  an  arbitrarily  small 


frequency  spacing  cannot  be  made  for  the  representation  of  arbitrarily 

short  time  signals.  Thus,  for  the  unwindowed  case,  I  ^dfT^A’Sc^ 2 

£ 

(2-4-6)  has  the  significance  of  energy  only  for  f ^  .  The  windowed 

case  will  be  treated  under  Question  (5)  below. 

(3)  "At  what  specific  frequencies,  f^,  should  I  Xjjpj.(f 1 2 
be  used  to  provide  a  physical  Interpretation  of  a  'time-variant  energy 
spectral  density'?"  Although  an  energy  density  function  is  a  convenient 
device  for  visualizing  the  concentration  of  energy  in  frequency,  it  has 
the  disadvantage  of  requiring  an  integration  operation  to  determine 
actual  energy  in  some  band  of  frequencies.  This  is  a  disadvantage 
because  the  energy  density  function  must  be  specified  for  a  continuum 
of  frequencies,  which  is  not  possible  for  an  arbitrary  length  signal 
segment  by  the  Heisenberg  Uncertainty  Principle.  In  addition,  the 
density  function  approach  involves  two  steps  in  the  calculation  of 
energy;  the  first  step  is  to  calculate  the  density  function  itself  and 
the  second  is  to  calculate  the  area  under  its  curve  for  some  given  set 
of  frequency  limits.  In  contrast,  the  vector  space  approach  requires 
the  calculation  of  expansion  coefficients  for  a  discrete  set  of  fre¬ 
quencies — a  one  step  operation.  In  addition,  the  valuable  influence 
data  is  produced  for  the  small  extra  effort  of  division  by  a  scale 
factor  and  a  square  root  operation  as  Indicated  in  (3-1-96).  Therefore, 
the  vector  space  approach  with  its  discrete  frequencies  provides  a 
more  economical  way  than  the  density  function  approach  for  physically 
interpreting  the  change  of  energy  with  time. 


(4)  "What  are  the  effects  of  averaging  several  sets  of 

I  XDFT^f£;tk^2  from  dlfferent  tlffle  sub intervals  tfc?  That  is,  what  sig¬ 
nificance  can  be  given  to  the  expression 


I  ?" 

k»K^ 


(3-3-5) 


This  expression  represents  a  function  of  frequency  f  associated  with 

i  ~ 


a  number  of  contiguous  signal  segments  centered  at  ,  t  . . t_  . 

Based  only  upon  its  definition,  there  is  no  justification  for  its 


interpretation  as  the  energy  at  an  arbitrary  frequency  f.  due  to  a 


sampled  signal  in  the  observation  interval  [  t  -  -r  ,  t_  +  -r-] .  However , 
there  is  a  sound  justification  for  the  proposed  vector  space  view  of  a 
series  of  approximation  vectors  x^Xt),  \  +1(t),...,  ^  (t) .  Each 
approximation  vector  has  a  component  along  a  specific  basis  vector 
ii>kA(t)  which  is  associated  with  the  discrete  frequency  f^  and  the  kth 
time  segment.  The  squared  value  of  each  component  ^s  the  energy  of  the 


approximation  for  the  frequency  f  .  These  squared  values  may  then  be 

X> 


averaged  from  several  approximation  vectors  to  determine  an  average 


energy  associated  with  frequency  f^.  As  mentioned  under  Question  (2) 


above,  only  certain  equally  spaced  frequencies  qualify  for  this  inter¬ 
pretation.  Therefore,  in  light  of  the  vector  space  concept  proposed 
here,  the  expression  in  (3-3-5)  is  recognized  as  the  sum  of  the  energies 

(t) 


at  f^  of  the  approximation  vectors  x^Ct),  x^  +^(t),...,xR 
provided  that  X^^(f t^)  is  Interpreted  as  the  sampled  version  of  the 


expansion  coefficient 
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(3-3-6) 


(5)  "What  Interpretation  should  be  given  to  the  family  of 
energy  spectra,  l2»  w^ien  the  si?1131  is  multiplied  by  a 

'weighting  window'  before  it  is  transformed?"  Outside  of  the  family 
of  windowed  expansions  proposed  here,  there  is  no  sound  reason  why  the 
signal  x(t)  should  be  multiplied  by  an  arbitrary  window  before  a 
Fourier  Transformation  and  then  treated  as  if  this  process  revealed  the 
true  frequency -distribution  of  energy.  It  would  seem  that  this  practice 
would  introduce  error  in  general.  The  typical  explanation  offered  is 
that  the  "windowing"  process  sharpens  the  spectral  representation  and 
facilitates  estimation  of  sinusoidal  frequencies  present  in  the  signal. 
Although  the  end  result  is  good,  the  approach  taken  is  unjustified. 

On  the  other  hand,  the  proposed  windowed  expansion  of  a  signal  in 
(3-1-4)  and  (3-1-5)  offers  a  sound  vector  space  approach  to  the  analysis 
of  energy  distributed  in  shaped  cisoidal  basis  functions.  The  radian 
frequencies  of  these  basis  functions  are  limited  by  the  rule  in  (3-1-7) 
which  is 
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This  shows  that  for  larger  values  of  Q,  a  larger  frequency  spacing 
between  basis  functions  is  required.  In  addition,  these  non-rectangular 
windowed  expansions  produce  sharp  frequency  representations  which  are 


adjustable  through  the  symmetric  function  a(t).  Thus  only  certain 
frequencies  and  certain  windows  are  allowed  in  these  basis  function 
expansions.  Finally,  it  should  be  remembered  that  the  non-rect angular 
windowed  expansions  do  not  in  general  possess  the  completeness  property 
but  their  sharp  spectral  representations  may  make  their  high  approxi¬ 
mation  error  tolerable.  The  problems  of  energy  spectral  analysis, 
which  involve  many  sinusoidal  signals  of  various  durations  and  even 
frequency  modulated  sinusoids,  can  often  be  made  more  manageable  by 
taking  advantage  of  the  sharp  spectral  representation  properties  of 
non-rectangular  windowed  expansions.  The  application  of  these  expan¬ 
sions  to  several  realistic  time-variant  spectral  analysis  problems  is 
the  subject  of  Chapter  4. 
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Chapter  4 

APPLICATION  OF  PROPOSED  APPROACH  TO  TIME-VARIANT 
SPECTRAL  ANALYSIS 
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In  order  to  apply  the  vector  space  concepts  of  Chapter  3, 
rationale  must  be  established  for  choosing  the  values  of  the  various 
parameters  in  the  signal  expansion,  such  as  the  signal  subinterval 
width,  T.  In  addition,  means  of  interpreting  the  appearance  of  two- 
dimensional  plots  of  the  energy  values  Ic^l2  must  be  established.  The 
purpose  of  this  chapter  is  to  develop  rationale  for  choosing  the  para¬ 
meters  of  the  expansion  and  rules-of-thumb  for  interpreting  the  plot 
of  | | 2 .  Finally  a  realistic  signal  example  is  given  as  an  illus¬ 
tration  of  time-variant  energy  spectral  analysis.  The  computer  programs 
used  are  given  in  the  Appendix. 

4.1  Types  of  Signal  Observations 

In  order  to  efficiently  use  the  proposed  windowed  expansions  of 
Chapter  3,  the  type  of  signal  observation  must  be  clearly  recognized. 
Four  types  of  observations  will  be  defined.  For  each  type,  a  specific 
relationship  exists  between  the  observed  signal  and  leading  and  trail¬ 
ing  "edges"  within  the  observation.  These  leading  and  trailing  "edges" 
will  be  shown  to  have  a  significant  effect  upon  the  interpretation  of 
the  plotted  energy  values  in  the  time-frequency  plane. 

As  a  first  step,  the  observed  signal  will  be  denoted  as  x(t) 
but  the  total  signal,  which  is  only  partially  observed,  will  be  denoted 

4-1 


as  y(t).  In  effect,  the  observer  knows  x(t)  only.  In  reality,  however 
there  is  a  signal  y(t)  that  accounts  for  both  observed  and  unobserved 
data.  Figure  4-1  clarifies  these  definitions. 


■i(t) 


FIG.  4-- \  THE  OBSERVED  S\G-NAL,  OC(-t') , 

AND  THE  TOTAL  S\GWAL 

The  total  signal  y(t)  is  assumed  to  be  square  integrable  and  to  have 

the  same  smoothness  properties  of  x(t)  mentioned  in  Chapter  1.  Again 

it  is  emphasized  that  x(t)  is  defined  only  during  the  observation 

interval  [t.,t-];  the  signal  x(t)  is  all  that  is  available  for  analysis. 
A  D 

As  in  Figure  3-1  of  Chapter  3,  the  observation  interval  is  partitioned 
into  N  subintervals  each  of  width  T. 

In  order  to  define  the  four  types  of  observations,  the  following 
terms  are  needed. 

(1)  Leading  Edge  point.  Consider  a  time  Instant  t^.  If 


there  exist  a  time  t,  and  a  time  t„  where 


such  that 


(4-1-6) 


(b)  x(t)  -  0  for  all  te  , 

then  Is  a  trailing  edge  point  of  x(t).  Figure  4-3  illustrates  the 
idea.  It  is  noted  that  x(t_)  -  0  which  holds  for  all  trailing  edge 

AX* 

points. 
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FIG.  4-3  EXAMPLE  OF  A  TRAILING  EDGE  POINT 


(3)  Edge  pair.  Consider  a  leading  edge  point  tT_,  and  a 
trailing  edge  point  t^.  If  there  are  no  leading  edge  points  or 
trailing  edge  points  in  the  open  interval  (tT_,t__),  and  x(t)  i  0  for 
some  t  in  this  open  Interval  then  t__  and  t_  are  said  to  form  an  edge 

LE  XE 

pair.  An  example  is  given  in  Figure  4-4  where  tT_  and  t__  con- 

LE^  TEj_ 

stltute  one  edge  pair  and  t__  and  t_  constitute  another  edge  pair. 
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FIG.  4-4  EXAMPLE  OF  2  EDGE  PAIRS  AND 
LEADING  EDGE  DOWINANT  OBSERVATION 


With  this  background,  the  four  types  of  observations  are  defined. 

(1)  Leading  Edge  Dominant  Observation.  This  is  defined  by 
the  fact  that  x(tA)  ■  0  and  that  there  is  one  leading  edge  point  not 
part  of  an  edge  pair.  An  example  Is  given  in  Figure  4-4  above. 

(2)  Trailing  Edge  Dominant  Observation.  This  is  defined  by 
the  fact  that  x(tg)  -  0  and  that  there  is  one  trailing  edge  point  not 
part  of  an  edge  pair.  An  example  is  given  in  Figure  4-5. 


F\G.  4-5  EXAMPLE  OP  TRAILING 
LDGE  DO  Ml  WANT  OBSERVATION 


(3)  Paired  Edge  Observation.  This  situation  is  defined  by 
the  fact  that  x(t^)  “  0  «  x(tg)  and  that  all  leading  and  trailing  edge 
points  lie  in  edge  pairs.  An  example  is  shown  in  Figure  4-6. 
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f\&.  4-6  EXAMPLE  OF  A  PAIRED  ED6E 
OBSERVATION 


1. 


(4)  Unpaired  Edge  Observation.  This  is  defined  by  the  fact 

that  x(t.)  +  0  and  x(t_)  i*  0  and  that  leading  and  trailing  edge  points 
A  B 

if  they  exist,  do  not  all  lie  in  edge  pairs.  An  example  is  given  in 
Figure  4-7. 


FIG.  4-7  EXAMPLE  OF  An  UNPAIRED 
ED6E  OBSERVATION 


Every  physical,  finite-length  signal  observation  falls  into  one  of  these 
four  categories.  It  is  important  that  the  type  of  category  is  recognized 
before  time-variant  spectral  analysis  is  begun  because  the  interpre¬ 
tation  of  the  energy  value  plot  of  | c^ I  2  will  depend  upon  the  presence 
of  leading  and  trailing  edge  points. 


4.2  Characteristics  of  the  Plotted  Values  of  Energy  in  the  Time-Variant 


Frequency  Plane 


As  an  illustration  of  the  preceding  discussion,  a  total  signal 
y(t)  is  considered,  where 


A  :  ‘a  < 

y(t)  ■  <  A/2;  t  «  tfl  and  t  -  t& 

0  :  t  <  t  and  t  >  t.  . 


(4-2-1) 


This  total  signal  Is  sketched  In  Figure  4-8 


F\&.  4-&  A  TOTAL  SIGNAL  UNDER  CONSIDERATION 

IN  (4-2-1) 

For  a  Leading  Edge  Dominant  Observation,  the  observed  signal  would 
likely  appear  as  shown  in  Figure  4-9,  where 
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FIG.  4-9  A  LEADING  EDGE  DOMINANT  OBSERVATION 

It  Is  noted  that  the  time  Instant  t  •  t  -  e,  with  £  an  infinitesimally 
small  positive  constant,  is  a  leading  edge  point.  The  constant  e  is 
needed  to  ensure  that  x(ta  -e)  ■  0  as  required  by  the  definition  of  a 
leading  edge  point.  Using  the  procedure  indicated  in  Figure  3-1,  the 
observed  signal  is  partitioned  into  N  equl-length  segments  each  of 
width  T.  It  is  required  that  T  be  such  that  N  be  an  integer. 
Specifically, 


must  be  an  Integer.  I£  T  is  chosen  such  that  a  segment  boundary 


happens  to  fall  exactly  at  the  leading  edge  point  t  *  t^  -e  ,  then  the 

subintervals  of  the  observed  signal  appear  as  in  Figure  4-10,  where 

t_  -  t.  has  been  defined  as  4T  for  convenience. 
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FIG-.  4-10  SEGMENTlON  OF  THE  OBSERVED  SIGNAL 


Applying  the  Rectangular  Windowed  Expansion,  the  coefficients  rkJ  are 
given  from  (3-2-24)  as 
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rVl 
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x(t) 


-iflt 

e  dt 


(4-2-4) 


Their  magnitude  squared  is  easily  seen  to  be 


Ir.  .I2 


ro  ;  k  - 1, 
J  t.  +? 


,  2lT  n*. 


(4-2-5) 


More  precisely, 


f  0;  k  -  1,  2 


A*T«a>0  !  *  -  3,  4 


(4-2-6) 


This  is  plotted  in  Figure  4-11  where  the  appearance  of  the  energy  values 
nicely  lends  itself  to  the  interpretation  as  a  "DC"  pulse. 


EXPANSION  ENERGY  VALUES  FOR  THE  OBSERVED 

SIGNAL  IN  FIG.  4-10 


Applying  the  Bartlett  Windowed  Expansion  from  (3-1-72)  with  w(t)  given 
by  (3-1-7)  and  Q  ■  2,  the  coefficients  are  found  to  be 

£t+7  / — 2TE=T|  t 

T  x(t)/l - j-*-  e  dt.  (4-2-7) 
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It  is  clear  again  that  c^A  ■  0  for  k  ■  1,  2.  For  k  *  3  or  4,  however, 
there  results 
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4-10 
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,4tt£ 
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K  e  1  dt.  (4-2-8) 


2  2 

Letting  r  -  — (t-t^),  dr  -  —  dt,  there  results 


Ck£  "  V  2 
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-ctft  e  1  2'  k  d t 


(4-2-9) 
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,4fr£,  . 

-  A  /I  e  T  k  j  e~i2*lT  dr  (4-2-10) 
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I  /1-T  cos(2tt£t)  dr 
J0 


(4-2-11) 


This  expression  is  reducible  to  an  interesting  formby  letting  z~*(1-t) 
so  that  2zdz  ■  dt  .  Performing  the  substitution,  it  is  found  that 
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z 2  cos  [2ir£(l-z^) ]  dz. 


(4-2-12) 


At  this  point  it  is  useful  to  consider  separately  the  case  for  £  -  0. 
Setting  £  -  0,  there  results 


c,  ■  AvW 
k,o 


2  . 
z  dz 


(4-2-13) 
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(4-2-14) 


Therefore,  the  energy  in  the  component  of  the  Bartlett  vindoved 
expansion  for  £  ■  0  is 
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(4-2-15) 


It  Is  noted  that  the  basis  function  component  for  l  *  0  is  just  the 
square  root  of  a  triangle,  as  sketched  in  Figure  4-12. 
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FIS.  +-12  BARTLETT  WINDOWED  EXPANSION  BASIS 

FUNCTION  FOR  A-O 

When  1^0,  (4-2-12)  can  be  simplified  and  then  integrated  by  parts. 
Specifically,  the  integral  is 


2 

z  cos 


X 

[2iri(l-z2)dz]  ■  [  z2  cos[2ir|  i|z2)dz  (4-2-16) 
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where  u  ■  z;  du  ■  dz 


(4-2-17) 


(4-2-18) 


"  sin(2irl*lz2); 
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dv  ■  z  cos  (2ir|  l\z  ) 


(4-2-19) 


This  process  gives 
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z  cos(2ir|i.|z  )dz  -  ^ j  z  sin(2ir|i|z  ) 
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Finally,  when  the  following  variable  change  is  made,  a 
expression  is  produced.  In  particular,  with 


y  -  2/JTfz, 
dy  *  2/jTjdz  , 


It  is  found  that 


I  1 

II  8|£|3/2 


s(2/17T), 


where  S(t)  is  the  Fresnel  Sine  Integral  defined  as 


S(t) 


I. 


sin(-|  y2)  dy. 


Since  the  Fresnel  Sine  Integral  is  tabulated  [17],  (4- 
evaluated  for  £  j*  0.  Inserting  this  result  back  into 
found  that,  for  i  +  0 


(4-2-20) 

(4-2-21) 

tractable 

(4-2-22) 

(4-2-23) 

2)dy  (4-2-24) 

(4-2-25) 

(4-2-26) 

-25)  can  be 
(4-2-12),  it  is 
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so  that 
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Kil2  m  *2T  w28j  £|  3  S2(2/[IT  )  .  (4-2-28) 


Combining  (4-2-15)  and  (4-2-28) ,  the  Bartlett  windowed  energy  values 


are  given  by 
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(|)  A2T  ;  A  -  0 


S2(2/TTP 


8*21113 


A2T  ;  l  +  0 


(4-2-29) 


when  k  -  3,4.  These  values  are  listed  in  Table  4-1  for  the  4  largest 
values  of  energy. 

TABLE  4-1 

Bartlett  windowed  energy  values  in  (4-e-*9) 


Ul 

l«Mi* 

o 

0.89 

OdB 

\ 

1.49  *IO'1 

-  87.  e  dS 

2 

-35'. 7  ib 

3 

7.  BE  *10’* 

.40.6  dB 

Note:  These  values  are  the  factors  multiplying  A2T  in  (4-2-29). 


It  is  interesting  that  the  values  listed  indicate  an  abrupt  drop  in 
energy  from  A  •  0  to  A*  ±1.  This  shows  that  the  energy  values  are 


essentially  insignificant  compared  to  the  energy  in  the  £  -  0  component. 
This  is  a  very  nice  spectral  representation,  although  it  is  not  as  nice 
as  the  Rectangular  windowed  expansion  result  in  (4-2-6).  It  is,  however 
tolerable  in  that  a  human  observer  would  tend  to  estimate  the  presence 
of  a  constant  signal  since  the  £  -  0  component  contains  the  maximum 
energy.  A  sketch  of  (4-2-29)  is  given  in  Figure  4-13. 


FIG.  4-13  BARTLETT  WINDOWED  EXPANSION 
ENER6T  VALUES  FOR  THE  OBSERVED  SIGNAL  IN  FlG.A-IO 
It  is  noted  that  the  index  £  in  this  plot  does  not  represent  the  same 

frequencies  as  the  index  £  in  Figure  4-11.  Here  the  index  values 

i  -1,2,3,  for  example,  represent  the  radian  frequencies  of  ~  , 

12ir 

— •  In  Figure  4-11,  these  index  values  represent  the  radian 

-  .  2ir  3ir  .  4ir 

frequencies  — ,  and  “ • 

To  Illustrate  the  effects  of  leading  and  trailing  edge  points, 
a  situation  is  now  considered  in  which  no  subinterval  boundary  falls 
on  the  leading  edge  points  in  Figure  4-9.  The  leading  edge  point 
will  now  lie  within  a  segment.  Specifically,  a  new  value  for  N  will 
be  chosen,  N  -  5,  and  the  subinterval  width  T  will  be  such  that 


(4-2-30) 


which  causes  the  leading  edge  at  t  •  t  -  e  to  lie  one  quarter  of  a 
subinterval  width  before  the  point  t  ■  tj  as  shown  in  Figure  4-14. 
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FIG:  4-14  ANEW  PARTITION  OF  THE  SIGNAL 

The  most  significant  result  of  this  change  will  be  seen  in  the  pattern 
that  develops  for  the  magnitude  of  the  energy  values  in  the  Rectangular 
windowed  expansion.  Speclfically(  the  squared  magnitude  of  the 
Rectangular  windowed  expansion  coefficients  from  (4-2-4)  are  found  to 


0  ;  k  -  1,2 


a2I  ;  k  ■  4.5 

1,0 


(4-2-31) 


where  the  case  k  ■  3  is  discussed  shortly.  These  values  in  (4-2-31)  are 
identical  the  corresponding  values  in  (4-2-6)  except  for  k  *  3.  When 


k  -  3, 
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-  3  -JTt(t3  +  |>  3 

a/F  -r  e  sine  (7 •£)  . 


(4-2-35) 


Therefore  the  squared  magnitude  Is 


lr3  J2  *  A2T  (~)  sine2  (-|i)  , 


(4-2-36) 


which  is  sketched  in  Figure  4-15. 
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F 16.  4-15  SKETCH  OF  ENERGY  VMUES  IN  (4-2-36) 


It  is  noted  that  energy  appears  to  be  distributed  all  across  the 


frequency  axis  with  magnitudes  which  follow  the  envelope  of  the  sine  (.) 


function  sampled  at  integer  multiples  of  y 


Applying  the  Bartlett  Windowed  Expansion  from  (4-2-29) ,  it  is 


seen  that  the  energy  values  for  k  *  4  and  5  are  given  by 
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(4-2-37) 
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and  that 
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(4-2-38) 


when  k  -  1  or  2.  However,  when  k  -  3,  it  is  found  from  (4-2-7)  that 
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Making  the  variable  change 


t  -  T  (t-t3)  -  4 


(4-2-41) 
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there  results 
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so  that 


3/4 

|cj  |*  *  T1  1  f  e-j2’*TdT|*  .  (4-2-45) 
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The  Integral  Involved  here  was  evaluated  for  various  Integer  values  of 
1  by  numerical  Integration.  The  results  are  shown  In  Figure  4-16  along 
with  the  logarithmically  scaled  values  of  the  Rectangular  expansion 
coefficients  in  (4-2-36) .  The  Rectangular  expansion  energy  values  are 
indicated  by  the  triangles  while  the  Bartlett  expansion  energies  are 
given  by  the  asterisks.  It  is  quite  interesting  that  the  Bartlett 
expansion  energies  of  (4-2-45)  actually  fall  below  the  Rectangular 
expansion  energies  at  so  many  frequencies.  The  reason  why  there  are 
twice  as  many  triangles  as  asterisks  in  the  plot  is  because  the  Rect¬ 
angular  expansion  represents  frequencies  which  are  integer  multiples  of 

whereas  the  Bartlett  represents  frequencies  which  are  integer  multi- 

4  it 

pies  of  —  .  The  smooth  curve  connecting  the  asterisks  is  only  for 
comparison  purposes  and  is  not  to  imply  a  continuum  of  represented 
points.  It  is  noted  that  the  Bartlett  windowed  expansion 
exceeds  the  fall-off  rate  of  the  Rectangular  expansion  for  this  plot. 
Thus  the  Bartlett  spectral  representation  is  quite  pleasing. 

In  order  to  interpret  this  plot  it  is  very  important  to  under¬ 
stand  why  neither  spectral  representation  produced  the  desirable  result 
of  a  single  nonzero  energy  component  at  1  *  0  with  the  rest  of  the  com¬ 
ponents  zero.  This  is  termed  the  desirable  result  because  it  would 
instantly  suggest  to  the  observer  a  square-wave  pulse  encountered  in 
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Fig.  4-16  Comparison  of  Bartlett  and  Rectangular  Windowed 

Expansions 
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the  third  time  subinterval.  Instead,  a  "splashing-out"  of  energy 
occurs  across  the  entire  frequency  axis.  The  reason  is  simply  that  no 
basis  function  in  either  Rectangular  or  Bartlett  expansions  contains  a 
leading  edge  point  within  the  subinterval  of  its  definition.  Therefore 
neither  expansion  will  represent  the  signal  in  the  third  subinterval 
with  a  single  basis  function.  Furthermore,  both  expansions  seem  to 
need  significant  energy  contributions  from  several  low  frequency  com¬ 
ponents,  the  Rectangular  expansion  especially,  in  order  to  approximate 
the  signal.  Specifically,  the  Rectangular  expansion  needs  the  component 
at  frequency  2/T  to  have  an  energy  of  about  13  dB  below  the  maximum 
energy  value  at  zero  frequency.  At  the  frequency  of  3/T,  it  needs  a 
component  with  an  energy  of  about  20  dB  below  the  maximum  energy  value. 
On  the  other  hand,  the  Bartlett  expansion  never  needs  more  than  about  a 
-17  dB  energy  value  in  any  component,  relative  to  its  maximum  at  zero 
frequency.  Thus  the  Bartlett  windowed  expansion  is  shown  here  to  have 
a  more  desirable  spectral  representation  because  of  its  more  rapid 
decay  rate.  Although  these  spectral  representations  may  both  be  con¬ 
sidered  troublesome,  at  least  they  can  be  anticipated  when  leading  edge 
points  are  noticed  within  a  signal  segment.  Trailing  edge  points  can 
be  shown  to  have  identical  spectral  representation  effects  as  leading 
edge  points.  The  lesson  learned  here  is  that  the  signal  must  be  care¬ 
fully  scanned  in  the  time  domain  to  detect  where  leading  and  trailing 
edges  may  fall.  Once  this  is  done,  then  the  basic  subinterval  width, 

T,  may  be  chosen  and  any  leading  or  trailing  edge  points  lying  on  sub- 


interval  boundaries  must  be  specifically  accounted  for  in  the  time- 
frequency  plot  interpretation. 


Another  interesting  investigation  concerns  the  relationship 
between  the  subinterval  width,  T,  and  the  frequency  of  a  sinusoidal 
signal.  The  question  may  be  asked,  "What  is  the  minimum  number  of 
cycles  of  a  sinusoidal  signal  that  needs  to  be  enclosed  by  the  basic 
subinterval  in  order  to  achieve  a  nicely  'peaked'  spectral  representa¬ 
tion?"  The  answer  depends  upon  the  choice  of  window  function.  For 
simplicity  the  signal  is  given  as 
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x(t) 
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(4-2-46) 


where  the  frequency  parameter,  v,  is  not  necessarily  an  integer.  If  v 
is  an  Integer,  the  Rectangular  windowed  expansion  gives  the  pleasing 
spectral  representation  of  the  single  nonzero  energy  value  at  the 
correct  frequency  index,  as  discussed  in  Chapter  3  after  (3-2-49) . 
This  is  due  to  the  fact  that,  for  integer  v,  (4-2-46)  is  identical  to 
the  real  part  of  the  basis  function  <i>v  (t)  in  the  Rectangular 
expansion.  For  noninteger  v,  the  Rectangular  expansion  coefficients 
are  found  by  substituting  (4-2-46)  into  (3-2-24) .  There  results 
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The  last  expression  can  be  simplified  by  the  change  of  variable 
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Therefore  the  Rectangular  windowed  energy  values  are 
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It  Is  recalled  that  these  values  were  plotted  In  Figure  3-18  for 
v  -  60.5.  Here  it  is  of  interest  to  set  v  equal  to  small  values  and 
Interpret  the  resulting  plot  of  energies.  Before  reviewing  the  plots, 
however,  it  is  important  to  note  the  presence  in  (4-2-54)  of  the  term 
involving  the  product  of  two  sine  functions.  This  "cross-product" 
reveals  that  there  is  a  certain  contribution  to  the  energy  value  at  a 
specific  frequency  index  l  due  to  the  "Interaction"  of  the  tails  of 
two  separated  sine  functions.  Although  this  term  may  have  negligible 
effects  for  wide  separations  of  the  two  sine  functions,  that  is,  large 
values  of  v,  the  effects  for  small  values  of  v  are  quite  significant. 

In  addition,  the  cosine  factor  which  multiplies  this  "cross-product" 
term  in  (4-2-54)  may  have  some  noticeable  effect  upon  the  spectral 
representation.  An  example  is  given  in  Figures  4-17  and  4-18  of  the 
Rectangular  expansion  energy  values  for  v  «  10.25,  Figure  4-17  shows 
(4-2-54)  for  the  maximum  contribution  from  the  cross  product  term  while 
Figure  4-18  shows  (4-2-54)  for  the  minimum  contribution  from  the  cross 
product  term.  More  precisely.  Figure  4-17  corresponds  to  the  case 
when 


4ttv 
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while  Figure  4-18  corresponds  to 


tR  «  (2n+l)ir 
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(4-2-56) 


The  two  plots  appear  very  roughly  the  same  although  there  are  slight 
differences.  In  Figure  4-18,  the  minimum  cross  product  case  of  (4-2-56) 
the  energies  for  £  =  10  and  £  ■  11  are  practically  the  same  within  the 
vertical  resolution  stated  in  the  caption.  In  Figure  4-17,  the  maxi¬ 
mum  cross  product  case,  the  energies  at  £  ■  10  and  £=11  are  shown  at 
different  decibel  leves  but  the  vertical  resolution  given  in  the  caption 
is  smaller  than  that  of  Figure  4-17  which  may  account  for  the  difference 
In  fact,  the  only  significant  difference  between  the  plots  occurs  at 
£■0.  In  this  case,  (4-2-54)  becomes 
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When  (4-2-55)  holds,  this  reduces  to 
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lrK,0l2  "  A2T  sinc2(v)  (4-2-59) 

which  is  shown  in  Figure  4-16  to  be  more  than  16  dB  down  from  the  maxi¬ 
mum.  On  the  other  hand,  when  (4-2-56)  holds, 

lrK>0|2  -  0  (4-2-60) 


as  shown  in  Figure  4-18.  Thus  the  overall  effect  of  the  value  of  the 
cross  product  term  is  very  slight,  at  least  at  this  value  for  v. 
Smaller  values  for  v  will  be  discussed  shortly,  but  first  the  Bartlett 


windowed  expansion  coefficients  are  plotted  for  the  same  value  of  v. 
Using  the  expression  for  the  Bartlett  windowed  expansion  In  (4-2-7)  for 
the  signal  given  In  (4-2-46) ,  the  Bartlett  coefficients  are  found  to  be 
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For  k«R,  there  results 
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Making  the  variable  change 
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it  is  found  that 
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Using  Che  notation 
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It  follows  that 
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Again  a  "crosa-product"  term  is  seen  here  involving  the  product  I+I_ 
with  the  same  cosine  factor  as  in  (4-2-54) .  Figures  4-19  and  4-20  give 
the  logarithmically  scaled  values  of  (4-2-71)  with  v  ■  10.25  for  the 
two  cases  of  maximum  and  minimum  cross  product  contribution,  respectively. 
For  the  maximum  "cross  product"  contribution  in  (4-2-71) ,  (4-2-55)  is 
imposed  and  there  results 


which  is  plotted  In  Figure  4-19.  Vhen  the  minimum  "cross  product"  con 
trlbution  condition  in  (4-2-56)  holds,  (4-2-71)  becomes 


which  is  plotted  in  Figure  4-20.  When  1*0  in  (4-2-73),  inspection  of 
(4-2-69)  shows  that  I_*I+  for  all  v  so  that 


Ic^qI2  -  0  (4-2-74) 

in  Figure  4-20.  It  is  seen  that  Figure  4-19  and  Figure  4-20  are  prac¬ 
tically  identical  within  the  vertical  resolutions  indicated.  Thus  the 
cross  product  term  is  not  very  significant  for  the  Bartlett  windowed 
expansion  at  this  frequency.  It  is  surprising,  however,  to  compare 
the  Bartlett  windowed  expansion  to  the  Rectangular  windowed  expansion. 
Specifically,  comparing  Figure  4-17  with  Figure  4-19,  it  is  surprising 
to  find  a  tremendous  increase  in  spectral  sharpness  for  the  Bartlett 
windowed  expansion.  For  example,  although  both  plots  have  their  maxi¬ 
mum  at  1*10  (which  is  desirable  since  v*10.25).  the  Rectangular 
expansion  energy  at  1*12  is  only  4.6  dB  down  from  maximum  compared  to 
about  21  dB  down  for  the  Bartlett  expansion  energy  at  1*12.  On  the 
lower  frequency  side  of  the  peak,  the  Rectangular  expansion  shows  an 
energy  of  about  8  dB  down  at  1*8  compared  to  the  Bartlett  expansion 
energy  of  about  36  dB  down  from  maximum  at  1*8.  This  indicates  a 


significant  advantage  in  using  the  Bartlett  windowed  expansion  for  this 
frequency  of  the  sinusoidal  signal. 

The  spectral  sharpness  advantage  of  the  Bartlett  windowed 
expansion  is  also  evident  at  other  frequencies.  Figures  4-21  and  4-22 
present  the  Rectangular  windowed  energies  for  v=2.5  in  the  two  cases  of 
maximum  and  minimum  cross  product  term  contribution.  More  precisely, 
Figure  4-21  presents  the  energy  values  of  the  Rectangular  expansion  in 
(4-2-54)  for  the  maximum  cross  product  contribution  condition  of 
(4-2-55)  while  Figure  4-22  presents  the  Rectangular  expansion  energies 
for  the  minimum  cross  product  contribution  in  (4-2-56) ,  both  plots 
representing  a  sinusoidal  signal  with  v-2.5.  The  plots  are  somewhat 
similar,  but  there  is  a  marked  difference  in  the  decay  patterns  on  the 
high  frequency  side  of  the  maxima.  Figure  4-21  shows  subsidiary  peaks 
at  £”7,11  and  14  with  magnitudes  of  approximately  13  dB,  17  dB,  and 
19  dB  down  from  the  maxima,  respectively.  On  the  other  hand.  Figure 
4-22,  the  minimum  cross  product  contribution  case,  indicates  subsidiary 
peaks  at  1*7,  10,  and  13  with  magnitudes  of  about  11  dB,  15.5  dB,  and 
21.2  dB  down,  respectively.  Thus  the  cross  term  has  a  noticeable  effect 
on  magnitude  and  location  of  "sidelobes."  The  reason  for  this  is  that  for 
low  frequencies,  the  tail  of  the  sine  function  at  the  negative  frequency 
of  V--2.5  in  (4-2-54)  starts  to  merge  with  the  sidelobes  of  the  sine 
function  centered  at  the  positive  frequency  v-2.5.  This  causes  the 
cross  product  term  to  increase  in  magnitude  and  so  its  effects  will  be 
seen  to  a  greater  extent.  It  is  important  to  note  that  the  presence 
of  significantly  large  energies  at  1-1  and  1-4  in  both  plots  makes  the 
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Rectangular  expansion's  spectral  representation  somewhat  undesirable. 

In  fact,  once  the  signal's  frequency  becomes  smaller  than  10,  the 
Rectangular  expansion  begins  to  produce  a  troublesome  spectral  repre¬ 
sentation.  It  is  interesting  at  this  point  to  review  the  Bartlett 
expansion  plots  for  the  same  signal  frequency  of  v-2.5.  Figure  4-23 
presents  the  Bartlett  expansion  energies  for  the  maximum  cross  product 
term  while  Figure  4-24  presents  these  energies  for  the  minimum  cross 
product  term  contribution.  The  difference  between  these  Bartlett 
spectral  representations  and  the  Rectangular  spectral  representations 
is  quite  pronounced  now.  For  example,  comparing  the  maximum  cross 
product  term  plots  in  Figure  4-21  and  Figure  4-23 ,  it  is  readily 
apparent  that  the  Bartlett  windowed  expansion  is  sharper  and  has  a 
more  rapid  decay  rate  than  the  Rectangular  windowed  expansion.  In  par¬ 
ticular,  the  Bartlett  expansion  energies  at  £-0  and  £-4  are  down  about 
14  dB  and  18  dB,  respectively,  from  the  maximum  whereas  the  Rectangular 
expansion  energies  are  down  only  7.7  dB  and  4.3  dB  from  their  maximum 
value  at  the  same  frequencies.  In  addition  the  energy  values  have 
decreased  by  about  41  dB  from  £-2  to  £-18  In  the  Bartlett  expansion  of 
Figure  4-23  compared  to  only  about  24.8  dB  over  the  same  frequency 
band  in  the  Rectangular  expansion  of  Figure  4-21.  The  difference  is 
even  more  pronounced  for  the  minimum  cross  product  term  cases  of 
Figure  4-22  and  Figure  4-24.  At  £-4,  the  Bartlett  expansion  of 
Figure  4-24  is  down  about  24.5  dB  from  maximum  whereas  the  Rectangular 
is  down  about  only  1.4  dB  from  maximum  at  £-4!  Furthermore  the  decay 
rate  of  the  Bartlett  expansion  in  Figure  4-24  is  about  62  dB  from  £-2 
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to  1*18  compared  to  about  only  31  dB  for  the  Rectangular  expansion  over 
the  same  frequency  band.  (The  deep  null  at  1* 6  in  Figure  4-22  cannot 
be  interpreted  as  the  indication  of  decay  rate  from  1-2  to  1-18.)  These 
facts  show  that  the  Bartlett  expansion  is  far  superior  in  spectral 
representation  to  the  Rectangular  expansion  when  the  basic  subinterval 
contains  only  2.5  cycles  of  the  sinusoidal  signal.  Even  though  the 
signal  is  not  represented  in  the  set  of  basis  functions  of  either  expan¬ 
sion,  the  Bartlett  expansion  still  gives  a  tolerable  spectral  repre¬ 
sentation.  Of  course  when  there  are  an  integer  number  of  cycles  of  the 
sinusoid  in  the  basic  subinterval,  the  Rectangular  expansion  produces 
an  ideal  spectral  representation,  as  mentioned  before.  However,  the 
Bartlett  expansion  is  still  useful  even  for  these  cases  of  integer  v. 

For  example,  when  v-1,  the  smallest  frequency  for  which  the  Rectangular 
expansion  gives  an  ideal  spectral  representation,  the  Bartlett  expansion 
is  still  somewhat  tolerable,  as  shown  in  Figures  4-25  and  4-26.  In 
Figure  4-25,  the  maximum  cross  term  case  is  plotted  while  in  Figure  4-26 
the  minimum  cross  term  case  is  given.  In  Figure  4-25  the  maximum  value 
occurs  at  1-0  and  the  next  significant  value  is  at  1*2,  4  dB  down  from 
this  maximum.  The  decay  rate  is  quite  rapid,  having  decreased  about 
38  dB  from  1-0  to  1*18.  Thus  it  is  likely  that  an  observer  would 
estimate  the  presence  of  a  sinusoid  with  frequency  between  zero  and  2/T 
in  the  signal.  On  the  other  hand,  the  minimum  cross  product  term  case 
for  v*l  plotted  in  Figure  4-26  has  a  pronounced  peak  at  1*2  because 
the  cross  product  term  causes  the  Bartlett  energy  value  to  vanish  at 
1*0  in  accordance  with  (4-2-74) .  Again  the  decay  rate  is  quite 


pleasing,  being  about  51  dB  from  2-2  to  £>18.  It  is  likely  that  here 
the  observer  would  estimate  the  sinusoidal  signal's  frequency  to  be  2/T, 
off  from  the  true  frequency  by  1/T  Hz.  This  may  be  troublesome  but  at 
least  the  spectral  representation  is  fairly  sharp. 

As  a  final  example  in  this  series,  the  Rectangular  expansions 
are  given  for  v-0.5,  the  maximum  cross  product  term  being  used  in 
Figure  4-27  and  the  minimum  cross  product  term  being  used  in  Figure  4-28 . 
In  both  cases  the  spectral  representations  are  not  very  desirable, 
especially  in  the  minimum  cross  product  term  case  of  Figure  4-28 .  In 
that  plot  the  energy  peaks  at  1*2  and  the  decay  rate  is  about  29  dB 
from  the  2-2  to  2.-19  values.  Similarly,  Figure  4-27,  the  maximum  cross 
product  term  case,  shows  only  a  little  better  representation  with  the 
maximum  value  at  2,-0,  off  by  only  0.5/T  Hz  from  the  correct  frequency. 

The  value  of  energy  at  2.-1  is  only  1  dB  down  from  this  peak  which  is 
acceptable.  Additionally,  the  decay  rate  is  about  30.6  dB  from  2.-0 
to  2.-19  which  is  about  the  same  as  that  for  Figure  4-28.  Thus  the 
Rectangular  expansion  for  v-0.5  is  not  very  safe  to  use  in  general, 
mainly  because  of  the  dominant  effect  of  the  cross  product  term.  It 
is  interesting  to  compare  the  Bartlett  expansion  energies  for  v-0.5, 
plotted  in  Figures  4-29  and  4-30.  Figure  4-29  gives  the  maximum  cross 
product  term  contribution  while  Figure  4-30  gives  the  minimum  cross 
product  term  contribution.  The  result  in  Figure  4-30  is  not  very 
desirable  because  the  maximum  occurs  for  2-2,  off  by  1.5/T  Hz  from 
the  correct  frequency.  The  result  in  Figure  4-29  is  much  more 
acceptable,  with  a  maximum  value  at  zero  frequency  and  a  rapid  de^ay 
rate  of  about  47  dB  from  2=0  to  2-18,  Despite  these  drawbacks,  the 


Bartlett  windowed  expansion  has  Important  advantages  to  offer  over  the 
Rectangular  windowed  expansion.  First,  the  sharpness  of  spectral 
representation  Is  quite  evident  and  Is  certainly  noteworthy  for  non- 
integer  values  of  v  down  to  about  v=2.5.  For  lower  non- integer  values 
of  v,  both  Bartlett  and  Rectangular  spectral  representations  become 
unacceptably  "contaminated"  by  the  cross  product  term.  Secondly,  the 
Increased  frequency  spacing  of  the  Bartlett  energy  values  over  the 
Rectangular  energy  values  seems  to  enable  a  smoother,  more  rapid  decay 
of  the  spectral  plot  away  from  the  maximum.  Thirdly,  the  cross  product 
term  itself  appears  to  be  more  attenuated  and  confined  to  a  narrower 
frequency  band  for  the  Bartlett  spectral  representation.  This  last 
point  is  evident  from  the  comparison  of  plots  in  the  v-0.5  case.  In 
conclusion,  it  should  be  evident  that  the  desirable  spectral  character¬ 
istics  of  the  Bartlett  windowed  expansion  for  an  integer  number  of 
sinusoidal  cycles  as  well  as  a  non-integer  number  of  cycles  in  the 
basic  subinterval  make  it  a  useful  windowed  expansion.  This  provides 
greater  flexibility  for  the  choice  of  the  basic  subinterval  width,  T. 
Using  only  spectral  sharpness  as  the  criterion  for  choosing  T,  the 
Bartlett  windowed  expansion  is  acceptable  with  as  few  cycles  as  2.5  in 
the  subinterval.  On  the  other  hand,  if  T  is  adjustable  so  that  an 
integer  number  of  cycles  can  be  contained  within  the  basic  subinterval, 
the  Rectangular  expansion  is  superior.  In  the  final  analysis,  the 
choice  of  T  depends  upon  knowledge  of  local  time  characteristics  in 
the  signal,  that  is,  the  number  of  cycles  per  second  in  a  given  time 
interval.  This  can  be  estimated  roughly  by  the  number  of  zero  crossings 
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Fig.  4-29  Bartlett  Windowed  Expansion  for  v*0.5,  Maxiimun 

Cross  Product  Contribution 
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Fig.  4-30  Bartlett  Windowed  Expansion  for  v-0.5,  Minimum 

Cross  Product  Contribution 


Secondly,  It  should  be  remembered  that  the  smaller  the  value  of  T 
chosen,  the  larger  will  be  the  frequency  spacing  in  Hertz  in  the  energy 
distribution  plot.  Finally  it  should  be  realized  that  the  windowed 
signal  expansions  express  the  signal  in  the  given  subinterval  by  a 
linear  combination  of  shaped  cisoldal  functions  so  that  interpretations 
of  energy  versus  frequency  plots  must  Include  a  reference  to  the  shape 
or  amplitude  modulation  of  the  basic  building  block  of  the  expansion. 

4.3  Example  of  Time-Variant  Energy  Spectral  Analysis 

This  section  discusses  the  analysis  of  an  FM-tone  modulated 
signal.  Specifically, 


x(t) 


A  cos  (2itv  t  +  Ssin  (2ttv  t)  ] ;  t.  <  t  <  t_ 
c  m  A  B 


unknown  ;  t  <  t.  and  t  >  t_  , 
A  B 


(4-3-1) 


where  6  is  known  as  the  modulation  index  defined  by 


6  - 


A  v . 
m  d 


m 


(4-3-2) 


In  this  context,  the  modulating  signal  is 


m(t)  ■  A  cos  (2irv  t)  ,  (4-3-3) 

m  m 

The  carrier  frequency  in  Hertz  is  vc,  and  the  signal  x(t)  is  another 
form  of  the  more  general  expression 


PIG.  4-31  INSTANTANEOUS  FREQUENCY  FUNCTION  OF  (4-3-5) 

A  sketch  of  the  signal  corresponding  to  this  instantaneous  frequency 

function  is  given  in  Figure  4-31. 


FIG.  4-32  AN  FM  TONE.  -MODULATED  SUSNAL 


The  Bartlett  windowed  expansion  will  be  determined  first.  The  Bartlett 
windowed  coefficients  from  (4-2-7)  are 


c 


k£ 


T 

rtk+2 


't  -- 
k  2 


(4-3-6) 


Making  the  change  of  variable 


T  -  J  (t-tk) 


dx  -  -  dt. 


(4-3-7) 

(4-3-8) 


there  results 


"k  l 


•A\  *<■¥ +  v  ^ 

-l 


x  e 


dx 


(4-3-9) 


For  simplicity,  let 
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NL 


(4-3-11) 
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and  N  la  a  positive  integer.  This  assumption  means  that  the  basic 
subinterval  is  exactly  as  long  as  N  cycles  of  the  carrier  frequency, 
follows  that 

*  (r  +  V  "  x(2T-r  +  V  (4-3-12) 

c 

■  A  cos  f2irvc  (^pr  +  t^)  +  8  sin  +  (4-3-13) 

c  c 

V 

«  A  cos  {wNt  +  2vvct^  +  8  sin  [»(~)Nr  +  2™^^]}.  (4-3-14) 

c 

Substituting  this  for  x(yx  +  t^)  in  (4-3-9),  there  results 

Cj^  ■  A^  e  * *  f  cos  {wNt  +  2irv{,tk  +  8  sin  [w(p)Nx 


+  2™^]}  /1-|t|  e"J2ir£t  dx  . 


(4-3-15) 


Notation  is  simplified  by  the  following  expressions.  Letting 


♦  -  2wv  t.  ,  (4-3-16) 

nc 


♦_  ■  2irv_t, 


(4-3-17) 


and 


■  ?  *'  ’T-  f-v  9jr*j>  wv*  v^  jv 


R"“ * 
c 


(4-3-18) 


where  R  is  not  necessarily  an  integer*  it  follows  that 

j4ir^t  1 

‘Tel  "  vf  m  (  cos  (»Nt  +  ♦  +0  sinfwRNx  +  4  ]} 


.  /iqRf  «‘3J’lT  4t 


(4-3-19) 


and  so  the  Bartlett  windowed  energy  values  are  found  to  be 


Icj^l2  -  1 1  cos{wNt  +  +0  sin  [irRNr  +  ♦  ]) 


•  e’J2,rAT  drj 2  . 


(4-3-20) 


Without  loss  of  generality*  the  observation  Interval  [tA»tg]  will  be 


defined  as 


‘A’0 


(4-3-21) 


tB  "  12T 


(4-3-22) 


so  that,  using  (4-3-10),  the  centerpoints  of  the  subintervals  are 


•  *  ■*. 


'  vi 


t  -I-JL 

C1  2  2v 


t  .  21 . 3N_ 
c2  2  2v 
.  •  c 


23T  _  23N 


(4-3-23) 


Substituting  (4-3-23)  Into  (4-3-16),  It  follows  that 


♦  ■  2irv_t,  -  irN 

ci  c  1 


4  ■  2irv  t.  “  3irN 

c2  ,  c  2  . 


vSl 


«Vi! 


♦  -  2irv  t_,  -  23irN. 

c^2  c 


(4-3-24) 


Similarly,  substituting  (4-3-23)  Into  (4-3-17)  and  using  (4-3-18),  there 
results 


4  ■  2wv  t,  •  wRN 

"l  "  1 


4  ■  2irv  t,  ■  3*RN 

■2  ■  2 


4  •  2irv  t,  ■  23irRN 

*12  "1J 


(*-3-25) 


As  an  exaaple,  the  Bartlett  energy  values  for  the  first  subinterval 
T 

centered  at  t^  ■  j  are  given  by 


IC1  J2  -  ~  IJ  CO*  {itNt  +  wN  +  0  sin  [trRNx  +  irRN]} 

-1 

.S=R  .‘J2,rtI  4t|*  (4-3-26) 

and  the  Bartlett  energy  value*  for  the  second  subinterval  are  given  by 

1 

lc2  J2  -  ^  |J  cos  {itNt  +  3wN  +  0  sin  [irRNt  +  3irRNj} 

-1 

/l- f t |  e"J2lTlT  dt|2  .  (4-3-27) 

When 

N  -  14 
0-48 

and 

R  “  (12H14)  » 

the  Instantaneous  frequency  function  In  (4-3-5)  and  In  Figure  3-30  has 
the  appearance  shown  In  Figure  3-32. 
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*  (4-3-28) 

a 


Mi 


their  corresponding  instantaneous  frequencies  are  given  in  Table  4-2. 


Table  4-2 


Centerpoints  of  Subintervals  and  Corresponding 
Instantaneous  Frequencies 


Centerpoint 

t. 


u10 


’ll 


'12 


Instantaneous  Frequency 
17. 864 /T 
16.828/T 
15.035/T 
12.965/T 
11.172/T 
10.136/T 
10.136/T 
11.172/T 
12.965/T 
15.035/T 
16.828/T 
17. 864 /T 


The  graphs  of  the  Bartlett  expansion  energy  values  for  the  first  three 
subintervals  are  given  in  Figures  4-33,  4-34,  and  4-35.  A  condensed 
plot  of  these  energy  values  for  the  first  eight  subintervals  is  given  in 
Figure  4-36.  In  Figure  4-36,  the  dashed  line  represents  the  Instantaneous 
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frequency  end  the  vertical  axes  Indicate  values  of  energy  greater  than 
-10  dB  down  from  the  anxious  of  all  subinterval  energy  values.  Thus 
a  "global”  scaling  is  used  in  Figure  4-36.  It  is  interesting  to  note 
that  the  significant  values  of  energy  follow  the  instantaneous  frequency 
quite  well.  An  Important  feature  is  that  the  magnitude  of  the  energy 
values  decreases  when  the  rate  of  change  of  instantaneous  frequency 
increases.  For  example,  from  the  first  subinterval  centered  at  t«t^  to 
the  next  subinterval  at  t^tj*  the  maximum  energy  value  decreases  by 
about  6  dB.  The  maximum  energy  value  continues  to  stay  near  this  lower 
level  for  the  third  through  the  fifth  subintervals  as  the  Instantaneous 
frequency  changes  from  17/T  to  11/T.  Then,  in  the  sixth  and  seventh 
subintervals  the  maximum  value  becomes  large  again,  essentially  equal  in 
magnitude  to  the  maximum  of  the  first  subinterval.  Thus  when  the  instan¬ 
taneous  frequency  changes  to  about  (1.5/T)  Hz  or  more  during  a  particular 
sub Interval,  the  peak  value  decreases  by  about  6  dB  from  the  maximum 
value  of  the  "steady"  tone  case.  Conversely,  when  the  instantaneous 
frequency  is  relatively  "steady,"  the  frequency  distribution  of  energy 
shows  a  strong  peak  at  or  near  the  instantaneous  frequency.  In  terms  of 
Figure  4-36,  a  rule  of  thumb  for  definition  of  "steady"  is  the  condition 
in  which  the  instantaneous  frequency  changes  by  less  than  l/T  Hz  over 
the  subinterval  being  considered.  Subintervals  1,  6,  and  7  would  qualify 
for  "steady"  tone  cases.  The  basic  fact  being  presented  here  is  that 
each  windowed  expansion  coefficient  Involves  the  average  of  the  product 
of  the  signal  and  a  shaped  cisoldal  function.  Thus  Instantaneous  infor¬ 
mation  is  being  "smoothed"  out  over  the  basic  subinterval.  As  a  final 
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Fig.  4-33  Bartlett  Wlndowad  Expansion  for  First  Sublntsrval 
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comment  on  Figure  4-36,  It  should  be  noted  that  a  basic  pattern  is  shown 
in  the  first  6  subintervals;  the  pattern  repeats  itself  in  a  mirror  image 
fashion  in  subintervals  7  through  12. 

The  Rectangular  windowed  expansion  is  now  presented  for  the 
sinusoidally  modulated  FM  signal.  From  (4-2-4),  the  expansion  coeffi¬ 


cients  are 


T 

vl 


_I  *  2 


-#lt 

x(t)  e  dt  . 


(4-3-29) 


Making  the  change  of  variable 


T  -  f  (t-t^) 


(4-3-30) 


dT  -  |  dt  , 


(4-3-31) 


there  results 


■  2v..  .1 


-f  1  ttk  {  x(fx  +  tk) 


(4-3-32) 


I**!1 


-?  1 1  x(fr  +  V  .-*tTdT|* 


(4-3-33) 


'.M'l’iyu  iff 


pi^FTmjn^rrr  - ‘  IJ’  ,"1  "T' 
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Substituting  x(-^r  +  t^)  from  (4-3-9)  and  using  the  simplifications  of 
(4-3-16),  (4-3-17),  and  (4-3-18),  it  follows  that 


Llc!' 


¥lj 


-1 

•  e-3"lTdrl2. 


cos  {irNr  +  <fi  +8  sin  [ttRNt  +  <p  ]} 
Ck  °k 


(4-3-34) 


These  energy  values  are  plotted  in  Figures  4-37,  4-38,  and  4-39  for 
k-1,  2,  and  3,  respectively.  The  same  parameter  values  of  the  Bartlett 
windowed  expansion  case  In  (4-3-28)  were  used  here.  Thus  the  first 
three  subintervals  produce  energy  values  which  are  shown  in  the  three 
figures  mentioned.  In  addition,  a  condensed  plot  of  the  first  eight 
subintervals  is  given  in  Figure  4-40.  In  this  condensed  plot  it  is 
seen  that  the  peak  values  of  the  energy  "track"  the  instantaneous 
frequency  quite  well.  In  fact,  because  of  the  smaller  frequency  spacing 
of  the  Rectangular  windowed  expansion,  the  condensed  plot  of  Figure  4-40 
seems  to  give  a  better  "track"  than  the  condensed  plot  of  Figure  4-36. 
Upon  close  examination,  several  interesting  differences  between  the 
Rectangular  and  Bartlett  expansions  emerge.  First,  when  the  energy 
plots  for  the  first  subinterval  (k-1)  are  compared,  i.e..  Figure  4-37  of 
the  Rectangular  expansion  with  Figure  4-33  of  the  Bartlett  expansion,  it 
is  seen  that  the  Bartlett  expansion  has  the  sharper  spectral  representa¬ 
tion.  This  is  shown  by  comparing  energy  values  at  similar  frequencies. 
Specifically,  the  Bartlett  energy  value  at  4-16  is  20.1  dB  below  the 
maximum  while  the  corresponding  Rectangular  energy  value  is  18.4  dB  below 
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1 1 


Its  Mil—.  A  similar  trend  la  seen  st  other  frequency  Indexes  also. 
It  la  noted  that  the  Instantaneous  frequency  la  fairly  constant  during 
this  first  tins  subinterval,  beginning  at  t-0  with  the  value  of  (18/T) 


Hz  and  ending  with  the  value  of  (17.46/T)  Hz  at  t*T.  Now  for  the  second 
subinterval  (lc-2) ,  the  Rectangular  expansion  of  Figure  4-38  Is  compared 
to  the  Bartlett  expansion  of  Figure  4-34.  The  Instantaneous  frequency 
at  the  centerpolnt  t-t2  Is  (16.828/T)  Hz  from  Table  4-2  and  the  Rect¬ 
angular  expansion  appears  to  have  the  better  spectral  representation 
with  Its  peak  at  (17/T)  Hz.  But  more  Importantly,  the  Rectangular 
expansion  plot  has  a  slightly  more  desirable  decay  rate.  At  the 
frequency  (14 /T)  Hz,  for  Instance,  the  Rectangular  windowed  energy  is 
down  21.3  dB  from  its  maximum  while  the  corresponding  Bartlett  windowed 
energy  is  only  down  16  dB  from  its  maximum.  On  the  high  frequency  side 
of  the  peak,  however,  both  expansions  decay  quite  rapidly.  In  particu¬ 
lar,  at  20/T  Hz,  the  Bartlett  energy  is  about  24.6  dB  below  Its  maximum 
and  the  Rectangular  energy  is  about  21.3  dB  below  maximum.  The  instan¬ 
taneous  frequency  decreases  from  Its  value  of  17.46/T  Hz  at  the  beginning 
of  the  second  subinterval  to  a  value  of  16 /T  Hz  at  the  end  of  this  sub- 
interval,  a  change  of  almost  1.5/T  Hz.  Finally,  in  the  third  subinterval, 
the  Rectangular  expansion  of  Figure  4-35  is  compared  to  the  Bartlett 
expansion  of  Figure  4-39.  The  Rectangular  energy  plot  is  clearly  more 
desirable  here  with  its  sharp  spectral  peak  at  15/T  Hz.  From  Table  4-2 
the  inatantaneoua  frequency  at  the  third  aubinterval  centerpolnt  (t-t^) 
la  15.035/T  Hz.  The  Bartlett  energy  plot  is  tolerable  with  maxima  at 


14/T  Hz  and  16/T  Hz,  although  ita  energy  at  12/T  Hz  la  only  18.2  dB 
down  from  maxi wna  compared  to  about  25.6  dB  down  from  maximum  for  the 
Rectangular  expansion  at  the  same  frequency.  Thus  the  Bartlett  expan¬ 
sion  has  a  slower  decay  characteristic  than  the  Rectangular  expansion 
for  this  case.  It  should  be  noted  that,  for  the  third  subinterval,  the 
instantaneous  frequency  changes  by  2/T  Hz,  starting  with  the  value  of 
16/T  Hz  at  the  beginning  and  ending  with  the  value  of  14/T  Hz  which  is 
the  carrier  frequency.  Figure  3-32  is  helpful  in  considering  these 
values  of  Instantaneous  frequency. 

The  analysis  performed  here  indicates  that  there  are  advantages 
and  disadvantages  for  using  any  windowed  expansion.  The  Rectangular 
expansion  is  more  sensitive  to  changes  of  instantaneous  frequency  than 
the  Bartlett  expansion  and  provides  more  precision  for  estimating 
instantaneous  frequency  with  its  smaller  frequency  spacing.  Also  the 
Rectangular  expansion  is  complete  for  arbitrary  finite-energy  signals 
while  other  windowed  expansions  are  in  general  not  complete.  On  the 
other  hand,  the  Bartlett  windowed  expansion  provides  a  sharper  spectral 
representation  than  the  Rectangular  expansion  when  the  Instantaneous 
frequency  is  fairly  constant  within  the  subinterval  being  considered. 
This  8 harp  spectral  representation  entails  a  rapid  numerical  decay  from 
maximum  as  well  as  a  monotonlcally  smooth  fall -off  on  both  sides  of  the 
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4.4  Example  of  the  Resolution  of  Two  Sinusoidal  Signals 


Since  non-rectangular  windowed  expansions  like  the  Bartlett  are 
expected  to  offer  sharper  spectral  representations  when  the  signal 
instantaneous  frequency  is  fairly  constant,  it  is  natural  to  expect 
these  properties  to  provide  good  resolution  for  two  or  more  sinusoidal 
signals.  In  this  context,  resolution  means  the  process  of  estimating 
the  true  frequency  and  energy  of  a  sinusoidal  signal  in  the  presence 
of  other  sinusoidal  signals.  In  this  section,  two  sinusoidal  signals 
are  assumed  to  occupy  the  same  subinterval  of  width  T.  The  energy 
values  for  both  the  Bartlett  and  Rectangular  windowed  expansions  are 
calculated  and  an  interesting  comparison  results. 

The  signal  under  consideration  is  given  by 

f  A  {cos(-~~  vt)  +  £  cos[^“  (v-A)t] } ;  It-tgif  j 


(4-4-1) 


where 


o;  It-t^f 


0  <  £  <  1 


0  <  4  <  v  . 


The  parameter  £  represents  the  ratio  of  the  amplitude  of  the  weaker 
sinusoid  which  has  frequency  (v  -  A)/T  Hz  to  that  of  the  stronger 
sinusoid  which  has  frequency  v/T  Hz.  The  parameter  A  is  the  frequency 
separation  factor  between  the  two  sinusoids.  Using  (4-3-9) ,  the 


\*3 


Bartlett  windowed  expansion  coefficients  are  given  by 


-/f  e_j^T  £tk  J^cfr  +  t^)  A  -  |t|  e"j2irtT  dx  . 

~l  (4-4-2) 

Substituting  the  signal  froa  (4-4-1)  into  this  expression,  there 
results 


\l  “  °i  k  *  K 


(4-4-3) 


a 

c„  .  A  /f.-l4?  *«fc  ♦  «,.>]  /ZR  €“J2fftTAT 


cos  [~(v-A)  (jr  +t^)  ] /l-| 


(4-4-4) 


TT  e-J1,r£TdT 


/sr  a4it1 

|  e“J  T  [  Ix  +  €  IA  1  (4-4-5) 

where 

Ij  ■  |  cos[wvt  +  ~v  t^]  A-|t|  e  ^2ir£r  dx  (4-4-6) 

-1 
and 

XA  "  j  co*fir(v"A)‘r  +  ^(v-AJt^]  A-|x|  e“^27r£T  dx  . 

-1  (4-4-7) 

Therefore, 

|cM|2  mi¥  1  *1  +  S  JA  I4  •  <*-«-« 

Choosing  t^*0  for  slaplidty,  the  values  of 
in  Figures  4-41  and  4-42  for  the  cases  when  £  ■  1/10  and 
C  •  1/16,  respectively.  In  both  figures,  the  sinusoidal  frequency 


IcjJ2  are  plotted 
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IMTIiTT  NINBONES  EXPANSION  FOR  TOO  SIMIS01H 

MX1IMI  VALUE  «  .220941  MINIUM  VALUE  >  7.30011E-00 

EACH  SHAKO  HORIZONTAL  SAM  HAS  MOTH  IN  N  Of  4.40059  LONEST  FHEQUENCY  •  * 

VERTICAL  RESOLUTION  IN  H  IS  1.49393  MINEST  FREQUENCY  •  24 

SINUSOIDAL  FREQUENCIES  NEK  18  /T  AM  14  /T 

TK  RATIO  OF  ANPL1TUK  OF  THE  14  /T  HI  SINE  TO  THE  IS  /T  HI  SINE  IS  .1 


figure  4-41  Bert let t  Windowed  Expension  for  Two  Sinusoids 

with  C  -  1/10 


parameters  were  chosen  as 


v  -  18 


and 


(4-4-9) 


A  -  4.  (4-4-10) 

This  choice  of  parameters  Indicates  that  a  dominant  18/T  Hz 
sinusoidal  burst  and  a  weaker  14/T  Hz  sinusoidal  burst  are  present 
in  the  same  subinterval  and  that  the  lower  frequency  sinusoid  has 
an  amplitude  which  is  smaller  by  the  given  factor  £.  It  is  Interesting 
that  the  weaker  sinusoidal  burst  Is  barely  detectable  when  £  -  1/16 
■  0.0625  In  Figure  4-42.  In  terms  of  the  energy  of  each  sinusoidal 
burst  alone,  the  weaker  sinusoid  has  an  energy  of  A2£2T/2  while 
the  stronger  sinusoid  has  an  energy  of  A2T/2.  Thus  the  ratio  of  their 
energies  Is  £2  and  the  rule  10iog^g(£2)  gives  -20.0  dB  for  £  »  1/10 
in  Figure  4-41  and  -24.1  dB  for  £  -  1/16  in  Figure  4-42.  It  is  noted 
that  even  though  the  frequencies  involved  here  are  contained  in  the 
set  of  represented  frequencies  of  the  Bartlett  windowed  expansion, 
the  spectral  properties  exhibited  are  not  ideal.  For  example,  the 
graphs  indicate  energy  at  other  frequencies  than  the  those  in  the 
signal  and  the  weaker  sinusoid  is  almost  hidden  from  view  when  £  ■  1/16 
The  reason  for  this  is  that  the  Bartlett  windowed  basis  functions  are 
shaped  sinusoidal  bursts  with  the  general  shape  sketched  in  Figure  3-2. 
In  contrast,  the  Rectangular  windowed  expansion  gives  the  ideal 
spectral  representation  for  this  case.  To  show  this  fact,  the 
expression  for  the  signal  in  (4-4-1)  is  substituted  into  the  expression 
for  the  Rectangular  windowed  energies  in  (4-3-33)  to  give 


-  0;  k  f  K 


(4-4-11) 


i-kJ2 


jl  j1  (*  co.1^1  v(f  +  tjj) 

-1 

+  AC  cos[^(v-A)(fx  +  t^)]}  e"jir*T  dx|2  . 


(4-4-12) 


Setting  t^*0  for  simplicity,  it  follows  that 

kKii2  ■  1 1  cos(irvx)  e~^lT  dx 

-1 

+  C  |  cos(ir(v-A)x]  dx|2 

-1 

■  ^4^  {  cos[ir(l-v)x]  dx 

0 

+  |  cos(ir(l+v)x]  dx  +  C  |  cos[ir(i-\>  +t 


(4-4-13) 


+  J  cos(ir(l+v)x]  dx  +  C  I  cos[tt(1-\>  +A  )x  ]  dx 
0  0 

f1  ”2 

+  Cl  cos[ir(l  +v  -A)x]  dx  ( 

0 


(4-4-14) 


— [  sinc(i-v)  +  sinc(i  +v  )  +C  slnc(l  -v  +  A) 


+  C  sinc(i  +v  -  A)]2  . 


When  v  end  4  are  both  Integers,  there  results 


(4-4-15) 


"  ¥  f  \v  +  \-v  +  «5l,(v  -A  )  +  *«£  ,-(v  -  A)  J2  * 


(4-4-16) 


This  Is  expected  since  the  two  sinusoidal  bursts  are  identical. 


BARTLETT  HINDONED  EXPANSION  FOR  TW  SINUSOIDS 

RAJtIHUH  VALUE  *  .123046  HINIHUH  VALUE  *  9.22401E-06 

EACH  SHADED  HORIZONTAL  BAR  HAS  N1DTH  IN  DB  OF  4.12S15  LONEST  FREQUENCY  *  6 

VERTICAL  RESOLUTION  IN  DB  IS  1.37505  HI6HEST  FREBUENCY  «  24 

SINUSOIDAL  FREQUENCIES  HERE  18.5  /T  AND  14.5  /T 

THE  RATIO  OF  AMPLITUDE  OF  TIE  14.5  /T  HZ  SINE  TO  THE  18.5  /T  HZ  SINE  IS  .0625 


Figure  4-44  Bartlett  Windowed  Expansion  for  18.5/T  Hz  and  14.5/T  Hz 

Sinusoids  with  £  ■  1/16 
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rectangular  hummed  expansion  for  tho  sinusoids 

IMXlHUH  VALUE  «  .101514  NINIHUN  VALUE  >  3.5843E-04 

EACH  SHADED  HORIZONTAL  BAR  HAS  HIDTH  IN  EB  OF  2.45212  LONEST  FREQUENCY  >  6 

VERTICAL  RESOLUTION  IN  OB  IS  .817374  HIGHEST  FRE8UENCY  >  25 

SINUSOIDAL  FREQUENCIES  HERE  18.5  /T  AND  14.5  /T 

THE  RATIO  OF  AHRUTUDE  OF  THE  14.5  /T  HZ  SINE  TO  TK  18.5  /T  HZ  SINE  IS  .1 


Figure  4-45  Rectangular  Windowed  Expansion  for  18.5/T  Hz  and  14.5/T  Hz 

Sinusoids  with  £  *  1/10 
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RECTAN6ULAR  NINDONED  EXPANSION  FOR  TNO  SINUSOIDS 

MAXIMUM  VALUE  «  .102486  MINIMUM  VALUE  >  4.36141E-04 

EACH  SHADED  HORIZONTAL  BAR  HAS  NIDTH  IN  EB  OF  2.37104  LONEST  FREOUENCY  «  6 

VERTICAL  RESOLUTION  IN  DB  IS  .700346  HIGHEST  FREOUENCY  «  25 

SINUSOIDAL  FREQUENCIES  HERE  1B.5  ,'T  AND  14.5  /T 

THE  RATIO  OF  AMPLITUDE  OF  THE  14.5  /T  HZ  SINE  TO  THE  18.5  /T  HZ  SINE  IS  .0625 


Figure  4-46  Rectangular  Windowed  Expansion  for  18.5/T  Hz  and  14.5/T  Hz 

Sinusoids  with  E  ■  1/16 


I  »  14  auch  more  prominent.  The  Bartlett  windowed  energy  plot  does 
not  present  such  a  deep  "valley"  and  is  therefore  less  desirable 
for  the  resolution  problem  in  this  case.  It  should  be  noted  that  the 
"valley"  effect  is  made  possible  by  the  interference  of  the  "tails" 
of  two  neighboring  sine  functions  in  (4-4-15) .  As  a  final  example, 
the  Bartlett  and  Rectangular  windowed  energies  are  compared  for  the 
case  when 

v  -  18.5  (4-4-20) 


and 


4*5 


(4-4-21) 


so  that 

v  -  A  -  13.5  .  (4-4-22) 

Figure  4-47  and  Figure  4-48  present  the  Bartlett  windowed  energies 
for  5  ■  1/10  and  £  »  1/16,  respectively,  while  Figure  4-49  and 
Figure  4-50  present  the  Rectangular  windowed  energies  for  £■  1/10 
and  5  ■  1/16,  respectively.  In  the  latter  two  plots,  the  deep 
"valleys"  occur  on  the  lower  frequency  side  of  the  weaker  sinusoid's 
frequency  of  13.5/T  Hz  and  not  between  the  two  correct  frequencies 
as  in  Figures  4-45  and  4-46  for  which  A  ■  4.  This  makes  the  task  of 
resolution  more  difficult  for  the  Rectangular  windowed  expansion  in 
this  case.  Again,  this  is  caused  by  the  degree  of  overlap  of  the  tails 
of  two  neighboring  sine  functions  in  (4-4-15).  It  is  surprising. 


however,  that  the  Bartlett  windowed  energy  values  plotted  in  Figures 
4-47  and  4-48  with  A  -  5  exhibit  practically  the  same  characteristics 
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DARTLETT  MNDONED  EXPANSION  FOR  TOO  SINUSOIDS 

NAX1NUH  VALUE  *  .124112  HINIRUH  VALUE  «  5.47451E-06 

EACH  SHADED  HORIZONTAL  BAR  HAS  NIOTH  IN  DB  Of  4.33547  LOWEST  FREQUENCY  *  4 

VERTICAL  RESOLUTION  IN  DB  IS  1.43102  HI6HE5T  FREQUENCY  •  24 

SINUSOIDAL  FREQUENCIES  HERE  18.3  /T  AND  13.5  /T 

THE  RATIO  OF  AMPLITUDE  OF  THE  13.3  /T  HZ  SINE  TO  THE  IB.5  /T  HZ  SINE  IS  .1 


Figure  4-47  Bartlett  Windowed  Expansion  for  18.5/T  Hz  and  13.5/T  Hz 

Sinusoids  with  £  »  1/10 
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ENEMY  -  SQM  NAGNITUK  OF  T  SIB  K,L 


RECTANGULAR  HINDONEO  EXPANSION  FOR  TWO  SINUSOIDS 

HAXIHUH  VALUE  «  .106082  NININUH  VALUE  »  2.42807E-06 

EACH  SHADED  HORIZONTAL  BAR  HAS  NIDTH  IN  EB  OF  4.64038  LOWEST  FREQUENCY  * 

VERTICAL  RESOLUTION  IN  DB  IS  1.5467V  HIGHEST  FREQUENCY  >  25 

SINUSOIDAL  FREQUENCIES  MERE  18.5  /T  AND  13.5  /T 

THE  RATIO  OF  AMPLITUDE  OF  THE  13.5  /T  HZ  SINE  TO  THE  18.5  /T  HI  SINE  IS  .1 


Figure  4-49  Rectangular  Windowed  Expansion  for  18.5/T  Hr  and  13.5/T  Hz 

Sinusoids  with  C  ■  1/10 


«s  the  Bartlett  windowed  energies  of  Figure  4-43  and  Figure  4-44  for 
which  A  -  4.  Thus  an  Increase  of  the  frequency  separation  factor 
has  produced  a  noticeable  difference  In  appearance  for  the 
Rectangular  energy  plot  but  not  for  the  Bartlett  energy  plot.  This 
Indicates  that  the  tails  of  the  sine  functions  in  the  Rectangular 
windowed  expansion  are  ouch  more  significant  in  value  compared  to 
the  corresponding  properties  of  the  Bartlett  windowed  expansion. 

The  Bartlett  properties  referred  to  are  the  rapid  decay  effects 
noted  previously  in  Chapter  3  and  in  Section  4.2. 

In  conclusion,  it  must  be  stated  that  neither  Rectangular  nor 
Bartlett  windowed  expansion  is  clearly  superior  for  the  resolution 
problem  in  all  of  its  aspects.  But  it  is  important  to  understand  the 
spectral  characteristics  of  various  windowed  expansions  so  that 
maximum  benefit  can  be  derived  from  their  different  advantages.  In 
the  resolution  problem,  the  closer  frequency  spacing  of  the 
Rectangular  windowed  expansion  must  be  carefully  weighed  against  the 
benefits  of  a  more  rapid  decay  from  maximum  offered  by  the  Bartlett 
windowed  expansion.  In  addition,  maximum  use  must  be  made  of 
knowledge  gained  by  a  careful  Inspection  of  the  signal  by  eye.  Such 
knowledge  consists  of  the  dominant  shapes  of  envelopes  within  the 
signal  and  characteristic  durations  of  these  shapes.  These  observations 
Influence  the  choice  of  window  and  the  choice  of  subinterval  length, T, 
respectively. 
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CHAPTER  5 

SUMMARY  AMD  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 

This  research  presents  a  sound  approach  to  time-variant  energy 
spectral  analysis  using  vector  space  concepts.  Chapter  1  briefly  des¬ 
cribes  the  motivation  for  the  time-frequency  description  of  signals 
along  with  some  basic  definitions.  Chapter  2  contains  a  survey  of 
previous  concepts  proposed  for  joint  time-frequency  signal  analysis. 
Several  important  problems  with  these  concepts  Involve  their  interpre¬ 
tation  and  implementation.  The  conventional  Finite  Fourier  Transform 
approach  is  presented  and  a  number  of  new  questions  are  raised  which 
have  not  been  answered  in  the  open  literature.  Chapter  3  introduces 
the  vector  space  approach  to  time-variant  energy  spectral  analysis. 

The  essential  points  covered  involve  the  segmentation  of  the  observed 
signal  Into  contiguous  signal  "pieces",  the  representation  of  each 
signal  "piece"  by  a  vector  in  a  specific  vector  space,  and  the  inter¬ 
pretation  of  the  components  of  this  vector  as  time  progresses.  A 
useful  family  of  windowed  signal  expansions  is  presented  which  provides 
theoretical  justification  for  using  weighting  windows  with  the  Finite 
Fourier  Transform.  The  questions  posed  in  Chapter  2  are  fully  answered 
providing  a  foundation  for  use  of  the  Discrete  Fourier  Transform  for 
"piece-wise"  energy  spectral  analysis.  Chapter  4  establishes  some 
ground  rules  for  applying  the  proposed  vector  space  approach.  A  real¬ 
istic  example  is  also  given  of  the  joint  time-frequency  energy  spectral 
analysis  of  an  FM  tone-modulated  signal.  The  appendix  lists  the 
computer  programs  used  in  Chapters  3  and  4. 
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A  number  of  Interesting  aspects  of  time-variant  energy  spectral 


analysis  remain  for  further  study. 

(1)  The  spectral  representation  effects  of  windowed  signal 
expansions  on  noisy  signals  is  an  important  area  of  study.  Strengths  and 
weaknesses  of  various  windows  may  be  revealed  when  specific  noise-like 
signals  are  added  to  the  signal  whose  time-variant  energy  is  being 
analyzed. 

(2)  The  windowed  expansions  may  be  optimized  for  specific  pur¬ 
poses  through  adjustment  of  the  function  a(t)  discussed  in  (3-1-24) . 

The  desired  placement  of  deep  nulls  in  the  spectral  representation  of 
sinusoids  may  be  conveniently  achieved  by  adjustment  of  the  shape  of 
a(t). 

(3)  The  problems  of  signal  detection  and  signal  parameter 
estimation  in  noise  might  be  attacked  more  fruitfully  using  the  proposed 
vector  space  concept  of  signal  representation.  In  general  this  may 
involve  windowed  expansions  of  random  signals  as  well  as  noise. 

(4)  The  general  area  of  bi-orthonormal  signal  expansions  such 
as  Bastlaans '  [2]  approach  should  be  more  thorougly  investigated  to 
determine  usefulness  for  representing  special  large  time-bandwidth 
product  signals. 


Appendix 


All  Integrals  in  the  text  were  evaluaced  by  the  numerical  Inte¬ 
gration  method  known  as  Filon's  Method  [18].  This  method  Involves 
Interpolation  of  the  integrand  at  a  number  of  equally  spaced  points. 

The  number  of  interpolation  points  is  chosen  based  upon  the  number  of 
cycles  of  the  sinusoidal  part  of  the  integrand  contained  within  the 
Integration  interval.  Specifically,  two  types  of  integrals  are  numeri¬ 
cally  Integrated,  I  and  I  ,  defined  by 

c  s 

ic  -  l  g(t>  cos  at)  dt  (A-i) 

and 

I,  -  J  g(t)  «ln  at)  dt,  (A-2) 


€  being  a  constant  for  the  Integration  process.  In  this  notation,  g(t) 
is  assumed  to  be  non-slnusoldal .  Tranter  [18]  shows  that  Fllon's 
method  is  just  a  modified  form  of  Simpson's  Rule.  Programs  were  written 
in  BASIC  and  used  on  a  Radio  Shack  TRS-80  computer.  Graphs  were  plotted 
on  an  Epson  MR-80  printer  with  132  character  capability.  The  computer 
had  48K  of  core  memory  and  was  used  in  conjunction  with  an  Aerocomp 
disk  drive.  Computer  programs  are  now  listed  by  corresponding  figure 
number. 


MRTliTT  EIPMSIOR  OF  SINGLE  SINUSOIDAL  TONE 
FILON  INTEGRATION  METHOD  USED 


40  BIN  N(123)|C(123) 

*0  PI»3. 141592*5358979324 
10  P2*2IPI 


100  INPUT  ‘NON  MNY  POINTS  TO  BE  PLOTTED'; NO 
120  1H-1.0E21 
140  NRM.0E21 


1*0  F2“0. 4342945 
ISO  SH*0 

190  INPUT  *NHAT  IS  THE  SINUSOIDAL  FREBUENCY  OF  THE  SIGNAL* :F9 


2*0  FOR  J*0  TO  123 
200  C(J)«0 


FlGOtMS  1-14  TWA.ufli 
3-H  ?-4». 


340  LPRINT  'FREOUENCY  ENERGY  VALUE 

BE*. *8*335218* 


300  GA* 1.312884582 
420  FOR  8*0  TO  <NQ-1> 
425  S0*<2»Q)-F9 


430  Sl*(2t8)+F9 
440  FOR  F1*0  TO  1 

4*0  IF_S0*0  AND  F1*0  THEN  480:ELSE  505  _ 

480  Al«0.***********7 
500  GOTO  1300 

505 JF  F1*^THEN  A0*8IABS(S0)tEL5E  A0*8IABS(S1) 
512  IF  AOCB  THEN  514:ELSE  520 

514  H*0.12S 

515  A0*8 


520  H*l/AO 
540  NI*lNT(A0/2) 


*20  FOR  N*1  TO  N1 
*40  AR*AR*H 

6*0  IF  F1*0  THEN  GOSUB  U20tELSE 
*80  R1*R1+Z 
700  AR«AR*H 

720  IF  F1*0  THEN  GOSUB  1120tELSE 
740  RO-RO *1 
7*0  NEXT  N 


820  IF  F1*0  THEN  GOSUB  !120iElSE 

B40  RO*ROMZ/2) _ 

940  IF  Fl*l  THEN  1000 
9*0  Al*<lR0tBE)+lRlt6A))tH 


1000  Bi*KR0tBE)*(8It6A>)tH 

1020  Cm*<AltAI)+<BltBm(2IASS<AltBI» 

1100  SOTO  1300 


1120  BOSUI  1180 
1140  2-21tCOS(SOtPItAR) 

1160  RETURN _ 

1180  IF  AR>1  THEM  MM 
1200  21-SQR<1-<A0S<AR))) 
1220  RETURN 


1240  60SUB  1180 
1260  2-2ltCOS(SltPItAA) 
1280  RETURN 


1320  SR-SIHC1I) 
1340  I-I+l 
1300  NEXT  8 


1370  IF  C101-0  THEN  C(0)-NN 

1380  FOR  1-0  TO  (NO-1) 

1400  A(l)-10tF2tL06(C<I)) 

1420  NEXT  I 

1440  HX-10IF2IL06<XHI 

1400  NN-10«F2tL06(NH) 

1480  REN  BBiINNlNP  OF  THE  HORIZONTAL  GRAPHING  NODULE 
1500  SP-INT1 124/NO) 

1320  RC-(HX-HN)/30 

1340  LPRINT  CHR0113) 

1560  LPRINT  *0  DO  *} 

1380  FOR  N-0  TO  29 

1600  QU-HX-(RCtN) 

1020  OL-OU-RC 

1040  IF  N-29  THEN  0L-200 

1600  PL-0 
1080  10-8 
1700  OP— 1 


1720  IF  IM  THEN  1400 

1740  IF  »3KINT(N/3))  THEN  1760I&SE  1840 

1760  DI-H)llNT«<RCtN)+.05)tl0)/10 _ 

1780  LPRINT  88) *88*) 

1800  IF  08)1-10)  THEN  LPRINT  CHW132)} 

1820  IF  8B-INT1B8)  THEN  LB-2>aSE  1900 


1840  FOR  L-l  TO  L8 
I860  LPRINT  CHRf(32)| 

1880  NEXT  L  _ 


2060  NEXT  J 


2140  LPRINT 
2160  KIT  K 

2190  LPRINT  'BARTLETT  jjINDONED  |XPAHSjON _ 

2200  LPOtlNT  'HAX1HUH  VALUE  •  '}XN}'  HINIHUH  VALUE  »  a|NH 
2220  LPRINT  'EACH  SHADED  HORIZONTAL  BAH  HAS  WDTH  IN  DB  OF  ';RC!3;* 

2240JLPRINT  'VERTICAL  RESOLUTION  IN  DBJS  ';RC;' _ 

2260  LPRI NT  "TOTAL  ENERGY  IN  THE  COHPONENTS  PLOTTED  NAS  'jSH 


_ S I 6NAL  FREQUENCY  «_'iF9j'/T' 

LOWEST  FREQUENCY  «  O' 
JUGHEST  FREQUENCY  JL 1  [2*  (NO-1 ) 


5  REN  RECTANGULAR  EXPANSION  OF  SINGLE  SINUSOID 
7  REN  FILON'S  INTEGRATION  HETHODJS  NOT  USED 
10  DIN  A(123),C(123) ~ 

20  PI«3. 141592 
30  Nl«1.0E21 


40  XI-1.0E21 

45  INPUT  *HON  HANY  POINTS  TO  BE  PLOTTED* ;NH 

47  F1«0. 434294402 _  __ 

50  INPUT  *NHAT  IS  THE  FREQUENCY  INDEX  NUHftER*;TT 

52  SH-0 

53  LB-1  _  _ 

59  LPRINT  *  L  VALUE  ENERGY  VALUE 

60  l«0 

70  FOR  3«0  TjMNMJ _ _ 

00  TL«ABS(J*TT) 

90  AR*TL 

100  60SUB  250 _  _ 

110  QH-9 

120  IF  J«0  THEN  QP«QN:60T0  170 

125  TIMBS(J-TT) _ _ 

127  IF  TIM  THEN  QP*1:80T0  170 
130  AR*TU 


160  QP-0 

170  A4I)«4  4QPSQP)  ♦  (ANION) ♦  (2IA9S (QPtQH)  M  /4 

190  CtI)«FH10tL06(A(l ) )_ _ 

200  IF  A(I)MI  THEN  XI-A(i) 

210  IF  A(IX«NI  THEN  Nl-A(I) 

215  SN»SIHA9SIA(I)) 


ITERATION 


234  RX«FII10I100<XI> 
236  HHFlllOILOG(Nl) 
240  SOTO  290 


275  REM  HORIZONTAL  GRAPHING  NODULE 
280  PRINT  'MAXIMUM  OF  GRAPH  MINIMUM  OF  GRAPH  >*,111 

200  REM  D£GINNIN6_0F  HORIZONTAL  8RAPHZN6  HDjjULE _ 

300  REN  MX,HN  ARE  THE  ASSUMED  MAI I MUM  AND  MINIMUM,  RESPECTIVELY 
310  REM  NH  IS  THE  NUMBER  OF  POINTS  PLOTTED 
320  SP-INTU24/NH) 


330  RC*(MX-MN)/30 
340  LPRINT  CHRI(IS) 
350  LPRINT  *0  DB 
360 


370  OU*HX-(RCtN) 

380  BL«flU-RC 
390  IF  N»29  THEN  OL 
400  PL«0 


420  OP-1 

430  IF  N*0  THEN  450 

431  IF  M«3l(lMT(N/3>)  THEN  432:  aSE  442 

432  D8>M)tINT((<RCtN>+.05><10)/10 

433  LPRINT  DB;*DBa; 

434  IF  IB) (-10)  THEN  LPRINT  CHR4(32)i 


435  IF  DD-INT(DB)  THEN  Lfl>2iELSE  450 

442  FOR  L»1  TO  LB 

443  LPRINT  CHRK32)  j _ 

444  NEXT  L 

450  FOR  J*0  TO  (NH-1) 

470  IF  C(JX«OU  AND  C(J))BL  THEN  PL-SPIJsaSE  540 


510  NEXT  K 
520  LPRINT  •»•} 
530  OWL 


540  NEXT  J 
550  LPRINT 
560  NEXT  N 


570  FOR  N*1  TO  5 
SBO  LPRINT 
590  NEXT  N 


600  LPRINT  'RECTANGULAR  NINBONED  EXPANSION 

610  LPRINT  ■MAXIMUM  VALUE  >  *|XI|*  MINIMUM  VALUE  • 

620  LPRINT  "EACH  COLORED  BAND  ISJMANGE  IN  DB  OF  *iRC»3 


630  LPRINT  'VERTICAL  RESOLUTION  IN  DB  IS  a|RC 

640  IF  L6»l  THEN  END 

650  REM  INFLUENCE  MEASUREMENT  MODULE 

660  XI— 1.0E21 

670  N1M.0E2! 

680  LPRINT  '  FREQUENCY 


SIGNAL  FREBUENCY  >  'jTT 


FREQUENCY 


MTIWW+TI/7) - 

MO  00-0.73 

MO  IF  F1*0  THEN  M89  1240:ELSE  USUI  1120 

— mr>i - 

920  00*0.73 

940  IF  FI-0  THEN  SOTO  900 

— WOUSJJTflO - 

935  OO-OAMY-Z) 

930  MTO  980 


970  OO-MKZ-Yl 

973  C(I)*l(R0IIE)«<RIt6A)+M>IH 

47«  00 ro  iJoO 

980  C(l)*((((R0tK)«<RltM)4A0)tH)(2>+<<C(im2) 
1100  60T0  1300 
1120  80SIII  1180 
1140  Z*ZltCOS<THtAR) 


1180  IF  00)0.73  THEM  00*0.73 
1200  Zl*SflR(l-(«BS(AR*0.25))> 

1220  RETURN 
1240DSUB  TIM 
1260  Z*ZltSIN(THIAR> 

1280  RETURN 

- 1300  TEITF1 - 

1302  C(I)*C(I)/2 

_ 1304  IF  C(I> >— IN  THEN  XIK1I) 

1306  TF  niX-NH  THEN  NH*CM) - 

1308  LPRINT  210;*  *jC(I) 

1320  SWIHC1I) 

- I340T«T*1 - 

1360  NEXT  0 

1380  FOR  1-0  TO  118-1) 

iwnrn)*ioiF2»LdG(C(T)) - 

1420  NEXT  I 

1440  NX*10IF2tL06(XN) 

1460  18I5I0IF2IL061W)  ~ - 

1480  REN  RE6INNIN6  OF  THE  HORIZONTAL  GRAPHING  NODLNJE 
1483  LPRINT  CH04U3) 

- lOTTPHNT 

1488  LPRINT 

1490  LPRINT*  ENEMY  -  3000  NA0NITU9E 

I3WSMNT190/N81 - 

1320  RC*(NX-NN)/30 

1360  LPRINT  *  0  81  *} 

13B0FDR  TWTir29 - 

1390  IF  N>0  THEN  LPRINT  *  *1 

1600  0U-HX-1RCIN) 


ItN  U.*W-nC 

1640  IF  N-29  TIEN  OL— 200 

1660  PL-0 


3060  NEXT  3 
2000  LPRINT 
2100  NEXT  N 


2107  FOR  J-l  TO  NO-1 
2100  FOR  l»l  TO  SP-1 
2111  LPRINT  CMRt(32); 

2113  KIT  K 
2115  LPRINT  'I'j 
2117  KXT  J 
2120  FOR  K*1  TO  5 
2140  LPRMT 
YlbO  next  K 

2100  LPRINT  *  MRTLETT  NINOONO  EXPANSION  FOR  EOWTION  (4-2-45)' 

2200  LPRINT  '  MXINUH  VOLK  ■  '|XHj'  NINIHUH  VALUE  ■  *|NR 

iPRiNr*  ^i«A0ErR)RT7WTCTjri^  TnroF'^fti3iT—EWEsr  frequency  .  o 

2240  LPRINT  '  VERTICAL  RESOLUTION  IN  N  IS  *|RC|'  NIBHEST  FREQUENCY  «  'J2NN0-1) 

22*0  LPRINT  '  TOTAL  ENERGY  IN  THE  CONPOKNTS  PLQTTEB  NAS  *(SH 


20  RED  TNO  SINUSOIK  IN  THE  SME  SUO INTERVAL  -  RECTANGULAR  HUMMED  EXPANSION 
40  IIN  A(V0)|C(90) 

- A4-PK.U1S92A5350070324 - 

00  P>2tPI 

100  INPUT  'HON  MANY  POINTS  TO  K  PLOTTED'|NO 
- lUMB*lMT4NQ/2) - 


* » j  • •  i  • *  » V " •  « *•*»*»  ' < ■ » '  • 


120  IR-1.0E21 
140  M-1.0E21 

- 160-F3*4342945  - 

100  SIM 

100  INPUT 'WAT  IS  TNE  FREQUENCY  OF  THE  1ST  SINE'lFO 

- 300-INPUT  !NHAT  IS -THE -FREQUENCY  OFFSET  OF -THE  3M84ME'HJ0 - 

220  INPUT'RATIO  OF  2MB  SINE  AHFIITUOE  TO  1ST  SIC  ANPL!TU8E';K0 
240  01BINT<(F0-(U0/2)H 

- 350-SOMfQ-UW - 


280  C4J)«0 

TOO  NEXT  3  . 

320  1-0 

340  LPRINT  'FREQUENCY  ENERSY  VALUE 

-420  FOR  8“<BI-N2}  TO  (0I*H2-1I 

■ 

440  AR-O+FO 

500  MNttS(O-FO) 
510  80SUB  1040 
— 520-R>R*Z  - 


560  60SUB  1040 
HO-MKUOUl- 
600  AR-AOSIO-SO) 


640  Mt*(KOU>  - 
660  C(I)B<R«R)/4 
740  IF  CUIMH 
260  IF-C(IK*NN 

THEN  IIP 
THEM -MR* 

780  LPRINT  0;* 

785  I>1«1 

•00  SOTO  1380 

1040  IF  ARXO  THEN  SOTO  1100 

- 1060  l-l - 

1080  SOTO  1200 
1100  AR-ARIPI 

- U304*4SIN(MUVNR - 

1200  RETURN 

1380  FOR  1-0  TO  INB-1) 

- i4oo~«  m*iotf2nmziw — 

1420  REIT  I 

1440  NXa10iF2HJB(IH) 


1480  REN  RE6INNIN6  OF  TIC  HORIZONTAL  0RAPMM6  NOSUi 
1490  LPRINT  *  ENERSY  -  SOM  MON  I  TUBE  OF  C 

1W  SP-INTIVO/I*) 

1520  RC*<M-M)/30 
1540  LPRINT  CNROUS) 

1560  LPRINT  *  0  18  ') 

1580  FOR  MO  TO  29 

1590  IF  N>0  THEM  LPRINT  1  *t 


A- 10 


1600  0IHUH8CIN) - 

1620  0L*8U-RC 
1640  IF  *29  THEN  BL— 200 

-466Q-PW - 

1680  18*7 
1700  OP— 1 

- W20-4F-IM  4WU900- - 

1740  IF  N*3NINT<N/3)1  THEN  1760:  &SE  1730 
1750  LPRINT 

- 6755  SOTO  1840 - —  ... 

1760  DD«<-niINT<<IRCIN>+.05)tl0)/10 
1780  LPftINT  08; *88'; 

- 480QTF-08X-WI-THEN  LPR4NT4HRHJ2); - 

1820  IF  8B*INT (08)  THEN  LQ*2:ELSE  1900 
1840  FOR  L*1  TO  LO 

- 1860  -LPRINT  OHR44J2); - 

1880  NEXT  L 

1900  FOR  J«0  TO  (NB-1 ) 

- 4920  IF  AUK*QU  AND  A4J)>QL  THEN  Pl*SP$J»ELSE  2060 - 

1940  IF  (PL-0P)*1  THEN  2020 
1960  FOR  K«0  TO  (PL-OP-2) 

- 1980  LPRINT  43NU32I; - 

2000  NEXT  K 
2020  LPRINT  ***; 

- 2040  OP*PL - - 

2060  NEXT  J 
2080  LPRINT 

- 2100  NEXT-N - - .....  - 

2105  LPRINT  *  I  !•; 

2107  FOR  J*1  TO  NO-1 

- 2109  TOR  K*1  TO  SP-1 - - - 

2111  LPRINT  CHRK32); 

2113  NEXT  K 

- 21 15  -LPRINT  ■I*; - - 

2117  NEXT  j 
2120  FOR  K-l  TO  5 

- 2440  LPRINT -  - 

2160  NEXT  K 

2180  LPRINT  *  RECTANGULAR  WINDOWED  EXPANSION  FOR  TNO  SINUSOIDS* 

- 2200  LPRINT  * -NAXIRUH -VALUE  ^  -H«»! - WNINUN  -VALUE  -  -HNN - 

2220  LPRINT  *  EACH  SHADED  HORIZONTAL  OAR  HAS  NIDTH  IN  EO  OF  *}RCI3j*  LOWEST  FREQUENCY  *  *;DI-N2 
2240  LPRINT  *  VERTICAL  RESOLUTION  IN  DO  IS  ■jRC;*  MISHEST  FREQUENCY  *  *;0I*N2-1 

- 2260  LPRINT  *  4INUS04BAL  -FREQUENCIES  HERE  ^ffO;VT - AND  - 

2270  LPRINT  *  THE  RATIO  OF  AMPLITUDE  OF  TIE  *j  IFO-UO>;VT  HZ  SINE  TO  THE  *;FO;*/T  HZ  SI*  IS  *;K0 
2275  LPRINT  '  TOTAL  ENERGY  IN  THE  COMPONENTS  PLOTTED  NAS  *;SN 
- 2389  {NO - — - - 


20  REH  TNO  SINUSOIDS  IN  SANE  DATA  HINDOO 

JO_OIH  A(90),C<90) _ _ 

60  PI-3. 14159265350979324 
80  P2-2IPI 


^Xft*etAT«r»jwrTV _ 

*4-94, 


100  INPUT  'HON  HMY  POINTS  TO  *  PLOTTED'lNO _ 

110  N2MNTlN0/2> 

120  IH--1.0C21 

J40  NHM.0f21 _ 

100  F 2-0. 434m5 
180  SH-0 

190  INPUT-HOI JS_THE  FREQUENCY JF  THE  1ST  SINE’jFO _ 

200  INPUT'NNAT  IS  TIC  FREQUENCY  OFFSET  OF  THE  2ND  SINE'fUO 
220  INPUT-RATIO  OF  2ND  SINE  AMPLITUDE  TO  1ST  SINE  AHPLITUK'jKO 

MO  DI-INT (IFO- (80/2)1/2) _ _ 

250  SO-(FO-UO) 

200  FON  J-0  TO  90 

280  C(J)»0  _ 

300  NEXT  J 
320  1*0 

340  LPftINT  'FREQUENCY _ ENER6Y  VALUER 

300  DE-O. 0803352180 
380  60*1.312884562 

420  FOR  O-IDHCMO  <01*2-11 _ 

424  T1«PIIU2»)«F0> 

420  T2«P1»U2»Q)*S0> 

428  T3«PIJ((2*8)-F01 _ 

430*  T4-PII(|2I8)-S0> 

440  A0-BHA8S((2I0)+F0)) 

445  H-I/AO 
450  N1-A0/2 
400  FOR  Fl-4  TO  1 
580  RR-0 
000  RI-0 
020  RO-O 

J40  .FDR  N-UO  JU _ 

000  AR-AR+H 

080  IF  Fl-0  THEN  OOSUD  1040iELSE  60SU8  1100 

_700_RJ«R!t! _ 

720  AR-AR* 

740  IF  F1*0  THEN  OOSUD  1040iEl.SE  60SUS  1100 
700  RO*RO+Z 


780  REIT  N 
800  R0-R0-(Z/2) 

120  ARM)  _ _  _ 

840  IF  Fl-0  THEN  60SUD  1040: ELSE  OOSUD  1100 
000  R0-R0MZ/2) 

1000  C(1>MUR0IDE)»(RH6A))1H)*C11) 

1020  OOTO  1220 
1040  OOSUD  1100 

JOOO  I«I1I< (C081T1MR)  1*1X01  (C0S(T2fAR) ) ) ) 
1080  RETURN 
1100  IF  AR>1  THEN  ARM 
112#  Zl-fQR(l-(ADS<AR)H 


1300  LPRINT  2Wj! 

1320  SH-SIHC(I> 

1340  I«I*1 

_ 1340  KXT  B _ 

1380  FOR  I«0  TO  IM-ll 
1400  A<I)«10IF2IL06<C<1») 

_ 142:  KIT  J _ 

1440  HX«I0«F2»L08<XR) 

1440  IM*10IF2I10G<M) 

_ 1 480 JEN  BEGINNING  0F_TK  HORIZONTAL  GRAPHING  NODULE _ 

1490  LPRINT  *  ENERGY  -  SflRD  IM6NITUK  OF  C  SUB  K,L' 

1500  SPMNT(90/NQ> 

_ 1520  RC-JIfXHW) /30 _ 

1540  LPRINT  CHRIUS) 

1540  LPRINT  '  0  DB  *| 

_ 1580  FOR  N*0  TO  29 _ 

1590  IF  N>0  THEN  LPRINT  *  *| 

1400  BU*HX-(RCtN) 

_  1420  0L*OO-RC _ 

1440  IF  N>29  THEN  01-200 
1440  PL>0 

_L480  LO-7 _ 

1700  OP-l 

1720  IF  N-0  THEN  1900 

_1740  JF  N*3t  (JNT(|/3)J JMEN  1740;ELSE  1750 _ 

1750  LPRINT  *I*j 
1755  60T0  1840 

1740  D8»(-1)HNT(((RC>N)*.05) OLQX/JLSl _ 

1780  LPRINT  DB}*BB*| 

1800  IF  BBH-10)  THEN  LPRINT  CMRt(32){ 

1820  IF  BB»INT(BB)  THEN  LQ»2»ELK  1900 _ 

1840  FOR  L*1  TO  LB 
1840  LPRINT  CHRt(32){ 

1880  KXT  L 

1900  FOR  J-0  TO  (NB-1) 

1920  IF  A(JX«OU  AND  A(J)>8L  THEN  Pl«SPUiELSE  2040 
1940  IF  (PL-OP)*i  THEN  2020 
1940  FOR  K«0  TO  (PL-OP-2) 

1980  LPRINT  CHRII32); 

ifiOOEXLK _ 

2020  LPRINT  *t*j 
2040  OP»PL 

2040  KXT  J  _ 

2080  LPRINT 

2100  KXT  N 

2105  LPRINT  * _ 1  I*» 

2107  FOR  J»1  TO  NO-1 
2109  FOR  K>1  TO  8P-1 
2111  LPRINT  0114(32)} 


4-13 


2113  KIT  K 
2113  MINT  * 
2117  KIT  i 


2121  FOR  1-1 
2140  LMINT 
2100  KIT  K 


2110  LMINT  1 
2200  LMINT  1 
020  LMINT  1 
2240  LMINT  1 
2200  LMINT  * 
H70  LMINT  1 
2773  LMINT  ’ 


MRTLETT  NINMNED  E1PONSION  FOR  TMO  SINBSOISS* 

M1IINJH  VOLK  >  *|IN}*  NININUN  VOLK  >  *|M 

EOCN  SHOOED  HORIZONTAL _M0  NOS  I11TH  IN  M  K  a>RCI3:a  _  LOKST  FKOKNCY  -  •J2I(DIHC) 
VERTICAL  RESOLUTION  IN  M  IS  a|ACja  HI6HEST  FKOKNCY  >  "j2t(II«CMl 

SINUSOIDAL  FKOKNCIES  NEK  ’jFOjVT  OK  a|(FO-W);*/Ta 
TNE  RATIO  OF  ONFLITUK  OF  THE  *»(FO-UO>tVT  N2  SIK  TO  JNE  'jFOj'/T  H2  81K  IS  ,|K0 
TOTAL  EKRBY  IN  TK  COHMKNTS  PLOTTER  NOS  a|SN 


537  SOTO  539 

538  NlaINT (41  <AK(  (2tfl)*F7) ) ) 

539  PRMT  *N1  »  *}N1 _ _ 

540  FOR  FI-0  TO  1 

550  RBI  Ft  CONTROLS  THE  REAL  (-01  OR  HM6IHARY  CD  INTEGRALS 

551  REN  J  CONTROLS  THE  PLUS  <«0>  OR  WNI1S  («1)  ARGUMENTS _ 

SM  FOR  M  TO  1 

580  Ti>P2t<fl«U(-l)H)t(N9/2>» 

HO  AR-l _ 

080  R1*0 
900  R0»0 

920  FOR  N»1  TO  N1 _  . _ _ 

940  ARWM+H 

940  IF  FI»0  THEN  G0SU8  1320J&SE  60SUB  1460 

900  RHU»Z  _ 

1000  MT-AR«M 

1020  IF  Fl-0  THEN  GOSUB  1S20:ELSE  60SUB  1460 

J040JNNR0+Z _ 

1060  NETT  N 
1080  INMNMZ/2) 

1100  AR»-1  _  _ 

1120  IF  F1«0  THEN  60SUB  1320: ELSE  GOSUB  1460 
1140  R0-R0+(Z/2) 

1160  IFJF1-1  THEN  1200  _ 

1180  C(I)>( ( IROIBE) ♦IRHGO) )  IH)«C(I) 

1190  GOTO  1540 

J200  IF  J«1  THEN  1260 _ _ _ 

1220  T3»C<mC<l) 

1240  C(I)»0 

1260^I1)»HIR0IBE)*IR1IS0))»H)4CJI) _ 

1280  IF  J«1  THEN  C(I)«(C(I)IC(imf3 
1300  GOTO  1540 

^320  GOSUB  1400  _  _ 

1330  RBI  T4-PHI  SU8  C  AND  T2-PHI  SUB  H 
1340  Z2»  ( 1-1)  I  J)  I HB9ISINI  (S2tAR)+T2> I  +T4) 

1360  Z«Z1*CKIIT1MR)*Z2) _ 

1380'  RETURN . 

1400  IF  AR>1  THEN  AR*1 

1420  Z1«S0RU-(ABS<AR») _ 

1440  RETURN 
1460  GOSUB  1400 

1480  Z2»<(-1HJ)I( IB9I5IN ( (S2IARI+T2) I +T4) 

1500  Z«Z1ISINI(T1IAR)*Z2) 

1520  RETURN 
1540  NEXT  J 
1560  NEXT  FI 
1580  CllXdl/l 

1MUE  mnaa  THER-UfeClIl _ 

1620  IF  C(I)<*NH  THEN  HKI!) 

1640  LPRIIT  218)'  •}CII> 

1660  >WN*CI1) _ 

1610  M*1 


1729  FOE  W>  TO  (MO-1) _ 

1749  A(I)-10§F2tL06(C(I» 

17M  NEXT  1 

1780  NI»10IF2IL08(IN)  _ 

1900  NN-10fF2tU)6<NN> 

1120  ISX  K6IHNIN6  OF  THE  HORIZONTAL  GRAPHING  NODULE 
1840  5P-INT1 124/NO) 

IttO  RC-(HHNI)/30  ' 

1N0  LFR1NT  CHR»(I5> 

1900  LPRINT  ‘0  H 
1920  FOR  N-0  TO  29 
1940  OU-NX-<ftCtNI 
1960  OL-GJ-RC 


1980  IF  N-29  THEN  OL— 200 
2000  PL-0 
2020  LO-8 


2040  OP-1 

20*0  IF  N-0  THEN  2240 

J080  IF  IW»(!NT(N/3))  THEN  2 1 00 1 ELSE  2190  _ 

2100  H*(-l)fTNT<((RCIN)+.09)ll0)/10 
2120  LPRINT  Hi1 •»••* 

2140  IF  0IX-10)  THEN  LPRINT  CMH(32) 


2160  IF  OI-INKDI)  THEN  LB-2:ELSE  2240 
2180  FOR  L-l  TO  LO 

2200  LPRINT  CHRI(32)|  _ 

2220  NEXT  L 

2240  FOR  7-0  TO  (NO-1) 

2260  IF  AUK-OU  A»  AJJ)X)L  THEN  PL«SPIJ»£LS£  2400 
2280  iF  (PL-OP)-l  THEN  2360 
2300  FOR  K«0  TO  (PL-OP-2) 

2320  LPRINT  CHR«(32)i _ _ 


2360  LPRINT 
2300  OP-PL 


2400  NEXT  J 
2420  LPRINT 
2440  NEXT  N 


2460  FOR  K-l 
2480  LPRINT 
2500  NEXT  K 
2520  LPRINT 
2540  LPRINT 
2560  LPRINT 
2580  LPRINT 
2605  LPRINT 
2610  LPRINL 
2615  LPRINT 
2620  END 


*MRTLETT  NINOONED  EXPANSION  FOR  SINUSOIDAL  FH  (TONE  HODULATION) 
•HAXIHUH  VALUE  •  •jXH}'  HINIHUH  VALUE  >  '|NH 

■EACH  SHADED  HORIZONTAL  BAR  HAS  WIDTH  IN  IB  OF  ,;RCI3;' _ 

"VERTICAL  RESOLUTION  IN  H  IS  -|RC|" 

■HODULATION  INDEX  BETA  -  'jl9 

^CARRIER  FREQUENCY  »  ■iN9!,/TV _ HODULATIAB  FREQUENCY  » 

■TOTAL  ENERGY  IN  THE  COMPONENTS  PLOTTED  HAS  *;SH 


TIME  INDi 

LONEST  FRE8UENCY  -  ';2 
HIGHEST  FREOUENCY  >  'j2l 


A-16 


20  REH  SINUSOIDAL  FH  TONE  MODULATION  HITH  RECTANGULAR  H2ND0NED  EXPANSION 
40  REN  FILON  INTEGRATION  NETHOD  USED 
60  DIR  A(I23),C(123) 

GO  PI-3.14139265356979324 
100  PG-PI/G 
120  P2-2IPI 

140  INPUT'HON  RANT  POINTS  TO  BE  PLOTTED'jNO 
160  XR--1.0E21 

190  MH.0E21 _ 

200  F2-0. 4342945 
220  SR-0 

240  N2-IRT (NO/2) _ _ 

250  INPUT'NHAT  IS  THE  TIRE  INDEX'jY 

260  INPUT'NHAT  IS  THE  HODULATION  INDEX  BETA*;B9 

280  INPUT'RATIO  OF  HODULATIM  TO  CARRIER  FREQUENCY* ;R9 

300  IHPUT'NUHBER  OF  CARRIER  CYCLES  IN  THE  DATA  N1ND0N'|N9 

320  NG-INT (N9/2) 

340  FOR  3-0  TO  123 
360  C(J)-0 


ASSoci  At«-C> 

f  friMLOS  M-  3*7 
THAt  U6-H  *4  -«f<* 


380  NEXT  3 
400  1*0 

420  LPRINT  'FREQUENCY  ENERGY  VALUE* 

430  REN  F8-1MBETAMR 

440  F8«1*(B94R9) _ _ _ 

442  REH  FILM'S  THETA  IS  PI/4  HERE 
445  BE-O. 735220407 

_450  60*1.252878002  _  _  _  _ 

455  REH  F7*N9t<l+(BETA)tR) 

460  F7-N9IF8 

470  T2*U2IY)-l)tPI»R9tN9 _ 

480  T4«( (2IY)-i ) IPI IN9 
490  S2-PIIR9IN9 

520  F0R  Q*(N9-N2)  TO  (N9+N2H1 _ 

530  H«l4t(ABSIQ*F7)))[(-l) 

532  IF  H<«0.25  THEN  GOTO  538 

534  H-0.25  _ 

336  Nl-8  ' 

537  60T0  539 

538  Nl«INT14»JADS(Q+F7)2) _ _ 

539  PRINT  *N1  «  *;N1 

540  FOR  Fl-0  TO  I 

350  RER  FI  CONTROLS  THE  REAL  (-01J1R  I RAGINARY  4-11  INTEGRALS 

351  REH  3  CONTROLS  THE  PLUS  (-0)  OR  RINUS  <-f)  ARGUMENTS 
560  FOR  3*0  TO  1 

580  T1«PII(Q*(((-1)[3)IN9)) _ 

860  AR»-I 
880  R1*0 

900  R0*0  _ 

920  FOR  N*1  TO  HI 
940  AR*AR*H 


960  JF  F 1  »0  _THEN  60SUB  1320:ELSE  608UB  1460 
f90  R1«R1+Z 
1000  NMMM 

1020  IF  F1«0  THEN  SOSUB  I 320 i ELSE  SOSUB  1400 
l040  RO^RO+Z 
1060  NEXT  N 

1060  KMHHZ/2) _ 

1100  AR— 1 

1120  IF  F1«0  THEM  60SU6  1320: ELSE  60SUB  1460 

JJ40  RO»ROMZ/2> _ 

1160  IF  Fl«i  THEM  1200 

1160  C(IJ»(((R0tBE)MRU60))tH)+C(l) 

1100  60TJM540 _ 

1200  IF  J«I  THEM  1260 
1220  T3»C(I)IC(t) 

1240  C(1)»0 _ 

1260  C(I)*(((R0t8E)+(Rll60l)tH)*C(I) 

1280  IF  J»1  THEN  C(I)«(C(I)tC(imT3 
1300  60T0  1540 
1320  21«1 

1330  REH  T4-PHI  SUB  C  AND  T2*PHI  SUB  M 
1340  Z2»  (  (-1)1  J)  t  ( (B9ISIN  ( (52IAR)  *T2)  1 +T4 ) 
1360  ZaZltC0S( (T1IAR1+Z2) 

1380  RETURN 

_1460Jlal _ 

1480  Z2*(H)[J)t((80ISIN((S2tM)iT2))«T4) 
1500  ZaZltSIM((TIIAR)*Z2) 

1520RETWH _ 

7540  NEXT  J 
1560  NEXT  FI 

1580  C(1)»C ( I m 6 _ 

1600  IF  CtI)>«XH  THEN  XH-C1I) 

1620  IF  C(IK>NH  THEN  NK(I) 

1640  LPRINT  Bjl _ *;C(I) 

1660  9MH+CU) 

1680  I«I+1 

J700  NEXT  0 _ 

1720  FOR  I«0  TO  (NO-1) 

1740  A(I)a10tF2tl06(C(l)) 

1760  NEXT_J_ _ 

1780  HX*10IF2tL06(XH) 

1800  HNa10tF2tL06<NH> 

1820  REH  BE6INNIN6  OF_THE_HORnONTAL_SRAFHJ NS 
1840  SMUT  (124/NO) 

I860  Rr«(HX-HH)/30 

1J80L  i?  ;NT_CHR»(  15) _ 

1900  LPRINT  *0  DB  'j 
1920  FOR  IN)  TO  29 

1940  flU«HX-(RCIN) _ 

1960  Ql-dU-RC 

1980  IF  N»29  THEN  flL-200 


2020  ie*a 

2040  OP-1 

_2060_IF.N«0  THEN  2240 _ 

2080  IF  N-3UINT1N/3))  THEN  2100:&SE  2180 
2100  DB»l-l)tINT(URCIN)+.05)tl0)/10 

2120  LPRINT  BBl'M’j _ 

2140  IF  DBM-10)  THEN  LPRINT  CHRK32); 

2160  IF  BB*INT (08)  THEN  LB*2:EISE  2240 

2180  FOR  L*1_T0  LB _ 

2200  LPRINT  CHRK32); 

2220  NEXT  L 

2240  FOR  J*0  TO  1N0-1J _ 

2260  IF  A(JK«9U  AND  A(J))OL  THEN  PL*SPIJ:£LS£  2400 
2280  IF  (PL-0P)*1  THEN  2360 

2300  FOR  K*0  TO  (PL-OP-2) _ 

2320  LPRINT  CTRK32); 

2340  KXT  K 

2360  LPRINT  ‘ft _ 

2380  OP*PL 
2400  NEXT  3 

2420_ LPRINT  _ _ 

2440  NEXT  N 
2460  FOR  K«1  TO  5 

_2480_ LPRINT  _ 

2500  NEXT  K 

2520  LPRINT  'RECTAN6ULAR  WINDOWED  EXPANSION  FOR  SINUSOIDAL  FN  (TONE  MODULATION) 

_2540_  LPRINT  ‘HAXIHUH  VALUE  »  *:XH;‘ _ HINIHUN  VALUE _i_!iNR _ 

2560  LPRINT  'EACH  SHADED  HORIZONTAL  BAR  HAS  WIDTH  IN  DB  OF  ';RC>3;* 

2580  LPRINT  ‘VERTICAL  RESOLUTION  IN  DB  IS  a)RC;‘  HI 

2605 JLPRINIJHODULATIPN  INDEX  BETA  »  ‘;B9 _ 

2610  LPRINT  ‘CARRIER  FREQUENCY  »  •iN9;"/T‘;‘  H0DULATIN6  FREQUENCY  »  *; 

2615  LPRINT  'TOTAL  ENER6Y  IN  TK  COMPONENTS  PLOTTED  WAS'*; SR 
2620  END 


TINE  INDEX 


LOWEST  FREQUENCY  -  ‘}(N‘ 
HIGHEST  FREQUENCY  «  (N9- 


•;N9IR9; ‘/T‘ 
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5  MISSION  j 

i  0/  i 

Jj  Rome  Air  Development  Center  ^ 

n  RA VC  p-ta.ni>  and  executes  research,  de.veZopme.nt,  test  and  v 

•  4  elected  acquisition  programs  in  support  of  Command,  Control  Q 

8  Communications  and  Intelligence  (C3I)  activities .  Technical  i 
%  and  engineering  support  within  areas  o f  technical  competence  Q 
J>  li>  provided  to  ESP  Program  0 ibices  ( POs )  and  other  ESP  $ 

elements.  The  principal  technical  mission  areas  are  \ 

»  communications,  electromagnetic  guidance  and  control,  sur-  § 
*S  veillance  o f  ground  and  aerospace  objects,  intelligence  data  \ 

*  collection  and  handling,  information  system  technology,  v 

K  ionospheric  propagation,  solid  state  sciences,  microwave 
%  physics  and  electronic  reliability,  maintainability  and 
S  compatibility. 


